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Prologue 

What is Physical Chemistry? Simply spoken, it is a scientific branch located between Physics and 
Chemistry. By using the principles of physics and mathematics to obtain quantitative relations, physical 
chemistry deals with the structure and dynamics of matter. These relations are, in most cases, either 
concerned with phase and chemical equilibrium, or dynamic processes such as phase transitions, reaction 
kinetics, charge transport, and energy exchange between systems and surroundings. 

To describe a physical chemical system and its dynamics or evolution towards an equilibrium state, only 
a limited set of variables of state is needed: volume, temperature, pressure and amount of material. The 
equilibrium state itself is based on the simple principle of minimizing the free enthalpy of the system. 
Free enthalpy is thus one of the most important concepts in physical chemistry. The change of free 
enthalpy is based on the change of two fundamental quantities: enthalpy which provides a measure for 
the energy of the system, and entropy which, qualitatively, characterizes the state of order of our system. 
Before going into too much detail here, let me stress one point to be kept in mind during the study of 
this book: physical chemistry is based on a small number of fundamental quantities and general physical 
concepts, which once fully understood by a student, it should be quite easy to pass an examination in 
physical chemistry without the need to explicitly learn a lot of detailed equations. 

This book, which tries to provide a complete overview of physical chemistry on the level of a Bachelor 
of Science degree in Chemistry, is organized as following: 

I already mentioned that physical chemistry deals with the quantitative description of chemical 
phenomena. Consequently, it affords some knowledge in mathematics. In Chapter 1, I try to briefly 
introduce all the mathematical concepts needed to understand the formalisms as well as the example 
problems presented in subsequent other chapters. 

Chapter 2, the main part of this book, deals with phenomenological thermodynamics, a topic also most 
dominant in all lectures in Physical Chemistry on bachelor level. Thermodynamics mainly deals with 
equilibrium characterized on a macroscopic level by variables of state, such as temperature, pressure or 
composition of mixtures. Here, the most important aspect related to the preferred equilibrium state of 
a system is its exchange of energy and entropy with its surroundings. First, a single component system 
and its phase behavior are discussed to introduce the fundamental concepts. We progress to binary 
systems and their phase equilibrium, an aspect which is very important for instance in the purification 
of chemical mixtures by distillation or crystallization. The chapter concludes with the equilibrium of 
chemical reactions and energetic aspects of chemical reactions. 
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In Chapter 3, I will briefly discuss the kinetics of chemical reactions. To keep the mathematics simple, 
only fundamental single-step reactions are treated here. As an example of more complex chemical 
processes involving multiple reaction steps and chemical equilibrium, this chapter will address the 
Lindemann formalism and the Michaelis-Menten kinetics, the latter being a very important topic in 
Biology (enzymatic reactions). 

Chapter 4 introduces the two fundamental aspects of electrochemistry: ion mobility and its relation to 
electrical conductivity, and the electrochemical equilibrium as the basis of the conversion of chemical 
reaction energy into electrical energy (battery). 

Chapter 5, finally, is a brief introduction into quantum chemistry and spectroscopy. Whereas the 
preceding chapters mainly treat systems on the macroscopic level, ignoring the detailed structure of 
matter, here it will briefly be shown how atoms and molecules are described in modern physics on a 
single particle level. Finally, this chapter will also provide a very brief introduction into spectroscopic 
methods, which are important for the experimental determination of molecular parameters such as bond 
lengths, strength of a chemical bond, etc. 

For illustration of the Physical-Chemical concepts and their relevance, each chapter (2-5) will also 
contain some quantitative example problems, presented in a box. 
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1 Mathematical Basics 


Physical Chemistry is frequently regarded as mathematically very complicated. Generally, students start 
with various levels of mathematical knowledge. Therefore, I consider it appropriate to start my book 
on the basics of physical chemistry with a chapter to briefly introduce the mathematical basics. Besides 
some technical terms, such as the total differential, in this first chapter of the book the students are 
introduced, as a refresher on fundamentals, how to solve simple differential equations, use the logarithm in 
mathematical calculations, use powers and trigonometric functions, as well as how to linearize functions 
of one variable, and treat physical-chemical data in terms of units or error analysis. 

1.1 Differentials and simple differential equations 

For example, if z =f(x,y) is a function of the variables x, y, then the total differential of z is given as: 

dz = (0) y dx+ (|), < e 9-i.d 

So-called “quantities of state”, i.e. basic functions of two variables x, y, fulfill the theorem of Schwarz: 

((0)(0)J y = ((0)(0)J, (Eq.1.2) 

, i.e. the result is independent of the sequence of differentiation, or, alternatively, a change of z = f(x,y) 
(= change of the quantity of state!) is independent of the route but only depends on the change of the 
variables x, y! 

Examples of such quantities of state are, for example, energy or enthalpy, whereas quantities dependent 
on the process itself, and therefore not quantities of state, are work or heat. 


Since physical chemistry needs just a limited set of simple differentials and integrals, I would like to list 
the most important ones for basic functions y(x) of one variable x below: 


y = x n -> ^ = n- x n 1 -» f ydx = ■ x n+1 

(Eq.1.3) 

y = e *->^. = e *-»J ydx = e x (e = 2.7182818) 

(Eq.1.4) 

y _ e k-x -> = /f ■ e k x -> f ydx = - ■ e k ' x 

^ dx J k 

(Eq.1.5) 

1 dy 1 f , 1 

y ~ X n ^ dx~ 1 ^ J ydx- 

(Eq.1.6) 

y = l^% = ~h^^ ydx = 

(Eq.1.7) 

y = sin x -> ^ = cos x -» f ydx = — cos x 

(Eq.1.8) 
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To determine the differential or integral of more complex functions, the following basic rules of 
differentiation and integration may come handy: 


i. How to differentiate a product (multiplication rule of differentiation): 

y = u(x) ■ y(x) -> — = u(x) ■ + v(x) ■ 

ii. How to differentiate a ratio (division rule of differentiation): 

. >. dv(x) , r ^du O) 

_ uW dy _ -“W^WlT 

^ v(x) dx v(x ) 2 

iii. Partial Integration: 

y = u(x) ■ / ydx = u(x) ■ y(x) — / i;(x) ■ du ^ dx 


(Eq.1.9) 


(Eq.1.10.) 


(Eq.1.11) 


Next, let us solve a couple of simple differential equations, i.e. determine the function itself if the 
differential is given. Note that these two examples will be found again in the chapter on reaction kinetics: 


Example 1: 


dy , 

i;=- k 'y 

(Eq. 1.12) 

To solve this equation, we first have to separate the two variables x, y: 


—= —k■dx 

y 

(Eq. 1.13) 

Next, we have to determine the integrals, using a set of starting variables x Q .y Q 
x, y as respective boundaries: 

and arbitrary variables 

f y - dy = j x —k ■ dx 

J y o y J J *o 

(Eq.1.14) 

X 

II 

(Eq. 1.15) 

In y — In y 0 = — k ■ x + k ■ x 0 

(Eq.1.16) 

y = y 0 ■ exp (—fc • (x - x 0 )) 

(Eq.1.17) 


Note that exp(x) corresponds to e x ! This example will be used again when we discuss first order basic 
reaction kinetics. 
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Example 2: 


d y i 2 

~r = -k ■ v z 

dx J 

(Eq.1.18) 

To solve this equation, again we first have to separate the two variables x, 

r- 

* 

1 

II 

^ 1 

(Eq.1.19) 

Analogous to example 1, we then have to determine the integrals, using a 
and arbitrary variables x, y as respective boundaries: 

set of starting variables x g .y g 

\ y \dy = j X —k ■ dx 

J y 0 y 2 J J *o 

(Eq.1.20) 

7 [-^] = 4 [ - fe 

y o 7 

(Eq.1.21) 

l l 

-1 — = —k-x + k-Xo 

y y o 

(Eq.1.22) 


This example will be used again when we discuss second order basic reaction kinetics. 


Within the context of this book, we will not need to solve more complex differential equations as shown 
here. However, the interested reader should keep in mind that there exist a variety of strategies to address 
more complex differential equations as might be encountered in a lecture on advanced physical chemistry, 
e.g. in quantum mechanics or complex reaction kinetics. Just to name one of these mathematical methods: 
partial fraction analysis: 

f 7-7-7 —-dx = f —^—dx + f ri C2 - dx i 

J (a-x)-(b-x) J (a-x) J ( b—x ) 

We will briefly sketch this method later in the book when we discuss the exact kinetics of 2 nd order 
reactions of two different chemical components. 

1.2 Logarithms and trigonometric functions 

Besides simple differential equations, in physical chemical calculations and derivations of important 
formulae it is essential that the student is capable of the logarithm rules as well as that he/she can handle 
trigonometric functions. These are therefore summarized in this subchapter. 

e a+b _ e a . e b J n ( a ■ ft) = \ n a + In b (Eq.1.24) 

e a ~ b = ^ l n (^) = In a -In b (Eq.1.25) 
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0 e a ) b =e a ' b «-> \n(a b ) =b-\na (Eq.1.26) 

sin 2 a + cos 2 a = 1 (Eq.1.27) 

sin a = cos(90° — a ) (Eq.1.28) 

sin(a + /?) = sin a ■ cos (3 + cos a ■ sin/? (Eq.1.29) 

cos(a + /3) = cos a • cos (3 — sin a ■ sin (3 (Eq.1.30) 


1.3 Linearization of mathematical functions, and Taylor series expansion 

In some cases, a plot of the function y(x) can simply be linearized if the reverse function is used, 
for example: 

y = y 0 • exp(-/c ■ x) lny = lny 0 - k ■ x (Eq.1.31) 
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0.4 




Figure 1.1: Linearization of an exponential (top) by plotting the logarithm (bottom) (Plots were prepared using 
Microsoft Excel (MS Office 2007)) 

Any function can be linearized if you use the Taylor series expansion given as: 

y(x 0 + x) = yO 0 ) + ’ (x - x 0 ) + y (§) Xo ' (x - x 0 ) 2 +1 ■ (g)^ ■ (x - x 0 ) 3 + - (Eq.1.32) 

One important (and the only!) example for a Taylor series expansion we use in this textbook concerns 
the logarithm, i.e.: 

ln(l - x) = In 1 + (j^_ o • (at - 0) +1 • (^)^_ o ■ (x - 0) z + - * -at (Eq.1.33) 

For small x, the logarithm In (1 — x) is simply given by -x, a relation we will use for instance to obtain 
quite simple formulae for colligative phenomena in dilute solutions. 
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1.4 Treating experimental data - Sl-system and error calculation 

(A) Units - the Si-system: 

The standard system of physical units (system international (French), SI) contains only a limited number, 
namely kg (kilogram) for mass, m (meter) for length, s (second) for time, A (Ampere) for charge per 
time (= electric current), and K (Kelvin) for (absolute) temperature. All other units can be expressed in 
terms of these basic units, and any student of science should be capable of this, since it comes in handy 
in (i) deriving simple formulae, (ii) understanding physical quantities, and (iii) calculating quantities 
based on certain variables. As one important example, let us consider several ways of expressing energy: 

1) potential energy equals force times length, and force is acceleration times mass, therefore 
energy should be: kg m s' 2 m = kg m 2 s' 2 = J (Joule) 

2) alternatively, energy is power times time, therefore W s (or, not SI, but more common: kWh) 

3) in electrostatics, energy is voltage times charge, therefore A s V (or keV (kilo electron Volt, 
not SI, with the elementary charge 1 e = 1.6 10 19 A s)) 

The most common derived SI units are: N (Newton) for force, Pa (Pascal) for pressure, J (Joule) for 
energy, W (Watt) for power, C (Coulomb) for electric charge, V (Volt) for voltage, U (Ohm) for electric 
resistance. 


(B) Error calculation: 


Experimental data are not perfect, that is, even if you repeat a given experiment several times, you will 

get at least slightly different results. If x x> x 2 , x 3> ., x N is a set of these different results obtained from a 

total of N single measurements, these data are evaluated as following: 


First, one has to calculate the average of the experimental quantity x, i.e. 



N 

i =1 




(Eq.1.34) 


To determine the reliability of this average, we then have to calculate the standard deviation. The apparent 
error of a single experiment is given as Ax t = (x) — x h and the standard deviation then is obtained by 
averaging the absolute values of all apparent errors, i.e. 


5 = 



with N > 1 


(Eq. 1.35) 


Finally, our experimental result is more reliable the larger the number of individual measurements N, 
which is taken into account within the statistical error of the average quantity (x) calculated from the 
standard deviation as: 


Ax = 


\ N yj N-(N-l) 


(Eq.1.36) 
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Within a written report of the experiment, the result must then be given as x = (x) ± Ax, providing 
both the average and the reliability. 


(C) Error progression: 


In most experiments, a given physical quantity depending on several individually measured parameters 
has to be calculated. As one example, let us consider a quantity z depending on two parameters x, y, that 
is z — f (x, y). In this case, the calculated error of z depends on the statistical errors of the two physical 
parameters x, y as following (error progression): 


A z = (j) m Ax + (y-) -Ay (Eq.1.37) 

UxJ (x),(y) KdyJ {x),{y) 

Note that the derivatives could be positive or negative, and therefore partially compensate each other, 
which does not make much sense. Therefore, one has to calculate the absolute values, i.e. 


A z 2 


(j) 2 ■ Ax + (f-) 2 ■ Ay, 


(x),(y) 


(Eq.1.38) 


or 


A z = 



(x),(y) 


2 


■ Ax + 



(x),(y) 


2 


' Ay 


(Eq.1.39) 
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As we have seen before, often physical-chemical data are represented in a linear relation, and the 
intercept and the slope provide the quantities of interest. To determine the error of this linear analysis, 
formally linear regression has to be used. This requires plugging all of your individual data points into 
the computer, and trusting the numerical calculation of the errors of intercept and slope. Additionally 
simple linear regression might ignore the statistical errors of individual data points. In practice the 
following approach is much more convenient and reliable. 

First, plot all your data points considering the statistical errors of each individual point (= error bars), 
respectively. 2 nd , draw two extreme lines which contain all data points with error bars, as well as an 
intermediate line. From this simple picture, you then can determine average slope and intercept, as well 
as the error of these two. The procedure is illustrated in figure 2. 



Figure 2: Graphical linear regression of data points including error bars 
(Plot was prepared using Microsoft Excel (MS Office 2007)) 
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2 Thermodynamics 

Phenomenological thermodynamics is mainly concerned with equations and variables of state, energy 
exchange at physicochemical processes, and phase- and chemical equilibrium. The goal of this chapter 
is to give a complete overview of all important concepts used in phenomenological thermodynamics, 
and provide some illustrative examples when appropriate. 

We will start with some basic definitions concerning physico-chemical systems and processes, then 
introduce the ideal and real gas equations of state, and finally discuss the fundamental principles of 
thermodynamics. Here, physical models for the heat capacity of ideal gases and solids will allow us to 
calculate the energy, heat and work exchange for any given process, even cycles such as the Carnot process. 

Next, the phase equilibrium of both pure components and binary mixtures will be addressed. For the 
latter, in the case of high dilution, simple formula for so-called colligative phenomena such as lowering 
of the freezing temperature, increase of the evaporation temperature, or osmotic pressure will be derived 
and discussed. Finally, we will talk about the equilibrium of chemical reactions, as well as about energetic 
aspects of chemical reactions. 

2.1 Definitions 

In thermodynamics, systems are classified according to their ability to exchange matter and/or energy 
with their environment: 



open closed isolated 

Fig. 2.1 : definitions of physical-chemical systems 
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An example for an open system, which can exchange both matter and energy with the environment, 
is an open beaker filled with boiling water (heat is lost to the environment, and matter is lost due to 
evaporation of water), whereas an example for a closed system would be a pot of boiling water tightly 
closed with a lid (only heat is lost). An isolated system would be, for instance, a closed thermostat bottle 
filled with hot tea: neither energy nor matter is exchanged with the environment. 

On the other hand, physico-chemical processes are classified as following: 1. Isochor, i.e. the volume of 
the system is kept constant (dV = 0), 2. Isobar, i.e. the pressure is kept constant (dp = 0), 3. Isotherm, i.e. 
the temperature is kept constant (dT = 0), or 4. Adiabatic, i.e. no exchange of heat with the environment 

(Q = o). 

Finally, we have to consider the different aggregate states of matter: solid (defined volume and shape), 
liquid (defined volume, undefined shape), and gas (neither volume nor shape is defined) 
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2.2 Gas equations 

2.2.1 The ideal gas 

In 1679, Boyle and Mariotte found experimentally that for a gas undergoing an isothermal process (either 
compression or expansion) its volume times its pressure remain constant: 

p ■ V = const (Eq.2.1) 


, namely 22.12 1 bar for 1 mole gas at temperature t = 0 °C. 

In 1800, Gay-Lussac discovered that for isobaric or isochoric processes, the volume or the pressure of a 
gas, respectively, increases linearly with temperature: 

1 st law of Gay-Lussac: V = Vq ■ (1 + cr 0 ■ t) (Eq.2.2) 


Here, Vq is the volume of the gas at standard conditions (p = 1.013 bar, t = 0°C). Importantly, different 
gases all have the same normalized thermal expansion coefficient: 



l 

273.16 


°c 


(Eq.2.3) 


, and we obtain therefore: 


V = V 0 '(l+ — L ~) = V 0 ■ P 7316+t ) = V 0 -a 0 -T 

U V 273.16/ U V 273.16 / U U 


(Eq.2.4) 


Here, T is a new temperature scale with unit “Kelvin”. Both scales can simply be converted into each 
other according to T/[K] = t/[°C] + 273.16. Note that, for the description of the isochoric processes 
of the ideal gas the Kelvin scale is more convenient, since simple proportionality is obtained. Also, since 
no negative volume can exist, it is obvious that T = 0 K (or t = - 273.16 °C) is the smallest temperature 
possible. Therefore, the Kelvin scale is also called the scale of absolute temperature, whereas the Celsius 
scale is based on the properties of one specific chemical, water, which freezes under standard pressure 
at t = 0 °C, and boils at 100 °C. 


For the isochoric process, in analogy the 2 nd law of Gay-Lussac can be formulated: 


2 nd law of Gay-Lussac: 


p = p 0 ■ (l + /?o ■ 0 


(Eq.2.5) 


Here, Po is the pressure of the gas at standard conditions (p = 1.013 bar, t = 0°C). Importantly, again 
different gases all have the same normalized thermal compressibility coefficient: 


Po= — = °C 


Po 273.16 


(Eq.2.6) 
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and we obtain therefore: 


p = Po • (l + = Po • (Irir) = Po • ft • 7- (Eq.2.7) 

To obtain the general equation of state of the ideal gas irrespective of the process, we next combine the 
laws of Boyle, Mariotte and Gay-Lussac, using formally a two-step process to change from a given set 
of state variables (po, To) to a completely different set of variables (pi, V 1 , Tj_): 


dT= 0 dp= 0 

(Po. Vb, r 0 )-♦ (pl V x , T 0 ) -♦ (pl Vi, Ti) 


(Eq.2.8) 


For the first isotherm step we use the law of Boyle, Mariotte to express the unknown volume V x : 


Pi ' Vx = Po ■ Vo 


V X =V 0 -^ 


Pi 


(Eq.2.9) 


For the 2 nd isobar step we used the first law of Gay-Lussac, i.e.: 


V~T 


= Vx 
Ti T 0 


v * = V t T * 


(Eq.2.10) 


We then can replace V x in (Eq.2.9), and finally obtain: 


El rp _ T T Po 

y 1 o - e 0 ■ — 

h Pi 


Pi'Vi _ Po'Vq 
Ti T 0 


(Eq.2.11) 


Since our starting conditions “0” as well as our final state “1” have been chosen arbitrarily, we can 
conclude that: 


Pl-Vl 


= const 


(Eq.2.12) 


According to experimental observations, this constant is given by the amount of gas in mole n, and the 
universal gas constant = 8.314 J/(mole K), and one obtains the universal law of ideal gases: 


p■V = n ■R -T 


(Eq.2.13) 
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Example 2.1.: 

Consider a small lab (20 m 2 area, height 3 m), where an automatic C0 2 -fire extinguisher, 
during a sudden fire in one of the hoods, suddenly releases 10 kg of C0 2 . What is the effective 
mass per m 2 acting on the windows of the lab due to the suddenly increased pressure (use the 
ideal gas equation of state for approximation)? 


Solution: The 10 kg of C0 2 (M = 44 g/Mol) correspond to 10000/44 = 227.3 Mol gas. The 
corresponding additional pressure therefore is 

n-R-T 227.3-8.314-293 nooo nooo N 

Pco 2 = — = -^- Pa = 9228 Pa = 9228 - 2 

The effective additional mass per m 2 corresponding to this excess pressure is then given by 


F = m - g => 9228 N = m - 9.81 — 

kg 


=> 


m = 940.7 kg 


Therefore, nearly 1 ton acts on 1 m 2 of our lab windows, so most likely the windows will 
explode. 
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Note that the ideal gas law (Eq.2.13) should be considered as a phenomenological law which is confirmed 
by experience, but it provides no explanation in terms of microscopic molecular motions. In the next 
section “kinetic gas theory” we will therefore discuss the ideal gas based on a microscopic single particle- 
based description developed by Maxwell, Boltzmann and others in the 19 th century 

2.2.2 Kinetic Gas Theory 

For the microscopic description of the ideal gas, we consider particles with mass m but no volume, which 
move without interactions in a strictly ballistic way between the rigid walls of a cubic sample container 
with side length a. Further, the number of gas particles is very large, so we can apply Boltzmann statistics. 
Note that a detailed derivation of the Boltzmann distribution is beyond the scope of this book, so we 
refer the reader to any common text book, section “statistical thermodynamics”. 

For this ideal gas system, the velocity distribution as a function of sample temperature and single gas 
particle mass m is then given, according to Maxwell and Boltzmann, as: 

^ (u)s (s77) 1S ' uZ ' ex p[“-0-] (EqZ14) 

It should be noted here that the exponential in this expression corresponds to the Boltzmann factor. 
Here, the Boltzmann constant is given as k B = 1.38 ■ 10 _23 //s. According to the Boltzmann factor, 
states of lower energy have a higher statistical probability of being occupied by the particles, as shown 
in figure 2.2. Multiplying this factor with the term u 2 > which is introduced in Eq.(2.14) due to a spatial 
integral over all possible directions of particle motion, one obtains a maximum velocity, which in 
dependence of temperature is shifted to higher values (see figure 2.3). 
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T = A E/k, N ± = 0.73 N 


T -> oo, N ± = N 2 = N/2 


Fig. 2.2 .: Boltzmann distribution at different T, 2-state-system 


In general, the ratio of occupation probability of an energetically excited state N 2 compared to the 
occupation probability of the ground state N 1 according to the Boltzmann factor depends on temperature 
and energy difference E 2 — E 1 = A E as: 


N2 

Ni 



(Eq.2.15) 
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Contrary to Fig. 2.2., our ideal gas particles in principle can assume any velocity and not just a quantized 
set of discrete values. However, the Boltzmann factor itself is not necessarily restricted to quantized 
energy states, but describes any probability ratio of populations — separated in energy by A E. 


We will meet this extremely important concept of the Boltzmann factor more frequently throughout this 
book. Here, let me just name one analogous example from kinetics, the well-known Arrhenius-equation: 

k = A ■ exp [- (Eq.2.16) 

with k the reaction velocity constant, and E A the energy of activation. The exponential term here describes 
the probability to reach the excited state via thermal excitation at temperature T. 

Coming back to our ideal gas treated on a microscopic level, let us discuss the shape of the Maxwell- 
Boltzmann velocity distribution at different temperatures in more detail: 



Fig. 2.3 .: Maxwell-Boltzmann velocity distribution at different temperatures (from cold (blue) to hot (red)) 

With increasing temperature, the maximum is shifted towards higher velocities, and the distribution 
becomes much broader. 
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The expression for the velocity distribution (Eq.2.14) allows us to calculate an analytical expression of 
the average squared velocity of a single particle as a function of sample temperature. We also need a 
relation between this average squared velocity and the macroscopic pressure to discuss the equation of 
state of the ideal gas in terms of microscopic molecular motions. Pressure is force per area, and force is 
the time derivative of the momentum, which itself is mass times velocity. Therefore, the pressure enacted 
by a single gas particle colliding with one wall of our cubic box is given as: 



^ d(m-u x ) ^ 


(Eq.2.17) 


The momentum of this particle is only changed at the collision, where it is reverted. Further, the time 
between two collisions is given by the velocity and the distance, and we therefore can express the 
momentum change with time as: 


f d(m-u x ) 

V dt 


) 


A (m-u x ) _ m-u x — m-u x _ m-u x 2 
At 2 d x /u x l x 


(Eq.2.18) 


For all particles within our cubic box, the overall momentum transfer related to the pressure is given 
by summing up the individual effects and replacing the squared velocity of the single particle with the 
average quantity: 

(Eq.2.19) 
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or 


p ■ V = N ■ m ■ ( u x 2 ) (Eq.2.20) 

Since all directions x, y and z have the same probability, the average squared velocity in 3d is given as 
(u 2 ) = (u x 2 ) + (u y 2 ) + {u 2 ) = i ■ <« x 2 ) (Eq.2.21) 


leading to 


p-V = N (Eq.2.22) 

According to the Maxwell-Boltzmann distribution (Eq.(2.14)), the average squared velocity of the 
particles can be calculated, corresponding to an average kinetic energy per particle of 

\-m-(u 2 ) = \-k B -T (Eq.2.23) 

Therefore, for a total number of gas particles N A = 6.02 ■ 10 23 > we obtain: 

p-V = N ■m- { ^Y = N A -k B -T = R-T (Eq.2.24) 


with the gas constant given as R = N A ■ k B = 8.314 J/(mole ■ K)\ 


Example 2.2.: 

Calculate the mean squared velocity of nitrogen molecules at room temperature, and compare 
your result to the velocity of sound. 


Solution: We simply use the mass per molecule m = 28 g/N A , and the following expression 


- • m • (u 2 ) = - ■ k B - T 


=> 


(u z ) 


2 \ _ 


3-feg-r _ 3-1.38-10 -23 -293 ; 


0.028/6.02-10 23 s 2 


= 260799 —t" 


Therefore, the mean velocity is 
Jiu 2 ) = V260799- = 510.7 - 

5 5 


Our result agrees well with the velocity of sound = 300 m/s. Note that in case of a sound wave 
a very large number of mainly nitrogen molecules has to move in the same direction. 
Therefore, this group velocity, i.e. the velocity of sound = 300m/s, has to be smaller than the 
statistical velocity of individual molecules. 
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2.2.3 Real Gases, isotherms and critical conditions 

For ideal gases, we ignored the particle volume as well as interparticle interactions, whereas in reality 
gas particles are molecules or atoms which possess both characteristics. The interactive forces are either 
caused by permanent dipoles or by fluctuating dipoles which depend on polarizability of a given molecule. 
For example, in the case of noble gases the fluctuating dipoles increase from He to Xe resulting in greater 
strength of attractive forces. To take these effects into account, we reformulate the ideal gas equation of 
state into the so-called Van-der-Waals-equation: 

(? + 7r) ‘ (“ “ b ) = (P + *0 ‘ £ “ b ) = R ‘ T (Eq.2.25 a) 

or, with v = — the molar volume: 

n 

(p + Jf) ' (v ~ b) = (p + n) ■ (v — b) = R ■ T (Eq.2.25 b) 

Note that the pressure the real gas enacts on the wall of a container, compared to the ideal gas, is 
reduced due to the attractive forces, i.e. p re = Pid — n < p^, whereas the sample volume itself is 
increased by the molar particle volume, i.e. v re = V i( ^ + b > . If you compare the ideal gas and 

the real gas laws at identical R a T, this means that: (i) the external pressure to keep the volume of the 
ideal gas constant has to be higher than in case of the real gas, and (ii) at given container volume V the 
sample volume accessible by the gas particles is larger for ideal gases than for real gases. The excluded 
volume b corresponds effectively to a short-range repulsive interparticle interaction. For spherical gas 
particles, like noble gas atoms, it can be simply calculated from the single particle volume according to 
the following sketch: 



Fig. 2.3 Excluded volume for one pair of spherical gas particles 
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In Fig. 2.3., the dotted circle corresponds to the excluded spherical volume which contains only a single 
particle at once. This volume is defined by two spherical particles at closest contact, and therefore has a 
radius of twice the single particle radius, and correspondingly eight times the single particle volume. In 
conclusion, since the sketch holds for two particles which exclude each other, the excluded volume per 
single particle corresponds to four times the particle volume, i.e. b — 4 ■ V Atom • 

The overall interaction between two real gas particles is described by the interaction pair potential 
sketched in figure 2.4. or 

*W = 4-*-[(t) 12 -(t)1 (e i- 2 - 26) 
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In Fig. 2.4., the parameter o m corresponds to the interparticle distance o where attractive and repulsive 
interaction energy compensate exactly, whereas the minimum position in the figure corresponds to the 
stable state with maximum net attraction. Therefore, o m is a measure of the particle size of excluded 
volume (see van-der-Waals-Eq.). Note that the depth of this minimum is parameterized by e. This 
parameter therefore provides a measure for the attractive interparticle forces or the internal pressure 
of the van-der-Waals-Eq. In the Eq.(2.26), one should also note that the attractive interaction energy 
based on dipole-dipole-attraction scales with inverse interparticle distance to the power of 6, whereas 
the repulsive interaction caused by direct overlap of the electron clouds of neighboring particles has a 
much shorter range and scales with inverse interparticle distance to the power of 12. Here, one should 
recall that the attractive Coulomb interaction energy between a positive and a negative charge scales 
with 1/a, and therefore has a much longer range than dipole-dipole-attraction. 



Both internal pressure and excluded volume have a strong effect on the shape of the isotherms of a 
real gas. Whereas the isotherm of the ideal gas at all temperatures according to Boyle, Mariotte (see 
above) is a simple hyperbolic function, the van-der-Waals Eq. allows for three different typical shapes of 
isotherms: above the critical temperature, the isotherms are similar to those of an ideal gas, and simply 
hyperbolic. At the critical temperature, the isotherm shows a turning point with slope zero, the critical 
point, whereas at lower temperature the isotherms show, with decreasing pressure and increasing volume, 
first a minimum and then a maximum. 


28 


Download free eBooks at bookboon.com 






Basic Physical Chemistry 


Thermodynamics 



r ftf273l h 273) 



Fig. 2.5 .: Isotherms of a real gas above, at and below the critical temperature (from: Gerd Wedler und Hans-Joachim Freund, Lehrbuch 
der Physikalische Chemie, p. 271,6.Auflage, Weinheim 2012. Copyright Wiley-VCH Verlag GmbH & Co. KGaA. Reproduced with permission.) 


Note that, for isotherms below the critical point, the part from D to C shown in figure 2.5., corresponds 
to an increase in sample pressure with increasing sample volume, which is physically impossible. The 
reason for this non-realistic behavior is that the van-der-Waals-Eq. (Eq.2.25) ignores the phenomenon 
of phase transitions completely. In reality, if we start at point G with a real gas with high volume, low 
pressure and subcritical temperature, and increase the pressure, the volume will decrease strongly until 
we reach point A. From there, in our experiment, in contrast to the theoretical prediction of the van- 
der-Waals-Eq., we will follow the horizontal line. This means the volume will decrease strongly while the 
pressure keeps constant, which is simply the case due to condensation of our gas at constant pressure. 
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Once we reach point B, all material has become liquid, and due to the low incompressibility of the 
condensed phase the pressure now has to increase strongly if we try to further decrease the sample 
volume. Here, it is important that the experimental horizontal connecting the pure liquid at B with 
the pure gas at A can be constructed from the van-der-Waals-curve by keeping the area below the two 
different curves, which corresponds to the overall work, identical. This is obvious since the energy change 
of the system from point A to point B corresponding to the condensation heat has to be independent of 
the process itself. At temperatures above the critical temperature, this condensation is not possible any 
longer, and there is no difference between gas and liquid. 


The macroscopic variables of state of the critical point, (px, V K , T K ), can be expressed as simple functions 
of the microscopic variables of the Van-der-Waals-Eq. (a, b) and vice versa, if you consider that within 
the isotherm at T = T K it is a turning point with slope zero. Therefore, the first and second derivative 
of p(V) have to be zero, that is, 




(Eq.2.27) 
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Solving the Van-der-Waals-Eq. and using, in addition, these derivatives, one obtains the following 
relations between the microscopic van-der-Waals-parameters of a real gas and its critical point: 


b =r v K 

(Eq.2.28) 

a = 9 --R-T K -V K 

(Eq.2.29) 

„ - 3 . R .I]L — — 

P K 8 V K 3 -b 2 

(Eq.2.30) 


In conclusion, it is possible to determine the microscopic characteristics of the real gas, such as effective 
single particle volume b and interparticle attraction a, via a macroscopic measurement of the pressure 
in dependence of volume at constant temperature, and consequent identification of the critical point. 


Example 2,3: 

Compare the volume and pressure of the ideal and real gas at room temperature and normal 
pressure for the two examples He and C0 2 . The respective van-der-Waals parameters are a = 
363.7 1 O' 3 , b = 0.0427 10' 3 for C0 2 , and a = 3.45 10' 3 , b = 0.0237 10' 3 for He (all quantities in 
SI-units). 

Solution: we compare the van-der Waals equation and the ideal gas equation to obtain the 
following expressions for the real and ideal quantities of state, respectively, using n = 1 mole. 

(;Pre + pri) -(y re -b) = R-T and p id ■ V id = R ■ T, 

Therefore 

Pre = Pid- 77-2 and V re = V id + b. 

v re 

Using ideal conditions T = 293 K and p id = 10 5 Pa, we obtain!^ = 8 ' 3 ^ 5 293 m? — 
0.02436 m 3 . Based on this result and the set of van-der Waals parameters, we get for 

He: V re =V id +b = (0.02436 + 0.0237 ■ l(T 3 )m 3 = 0.02438 m 3 

Pr,=m~i Z = (10= - tjUS) Pa = 99994 Pa 
C0 2 : V re = V id +b = (0.02436 + 0.0427 ■ 10 _3 )m 3 = 0.02440 m 3 
Pre = Pid ~ = flO 5 - 363 - 7 ' 10 U Pa = 99389 Pa 

t"re Hid 2 y 0 . 02440 2 ) 

We conclude that the ideal gas equation is quite valid for common gases at room temperature. 
Even for C0 2 which has considerable van-der-Waals attractive forces due to its dipole 
character, at room temperature the pressure is only reduced by much less than 1% compared to 
the ideal gas. 
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2.3 The fundamental laws of thermodynamics 

2.3.1 Definition of temperature 

If two systems A and B coexist in thermal equilibrium, and also B is in thermal equilibrium with a 3 rd 
system C, then all three systems are in thermal equilibrium and have the same temperature. This simple 
principle is the basis for temperature measurements as well as one foundation of phase equilibrium and 
chemical equilibrium. 

2.3.2 The first law of thermodynamics 

The first law of thermodynamics is the law of conservation of energy including, in contrast to conventional 
physics, heat as an additional means of energy exchange between system and environment. It is formulated 
for differential changes or measurable changes in total energy (7, respectively, as: 

dU = 8Q+8W ••• A U = Q + W (Eq.2.31) 

with U the internal energy, Q the heat flow, and W the work, i.e. typically volume work (expansion or 
compression). To define the heat flow, we refer to adiabatic processes, where the system is thermally 
isolated from its environment, and therefore A U( ad ^ = Q + W = W^ ad y therefore Q — —W -\- W( ad y 
Note that heat and work are not quantities of state, in contrast to U, but depend on the process itself, 
whereas A U only depends on the final and the initial state. As a consequence, Q and W also are not 
total differentials, but have to be formulated as or SQ or 8W, respectively. In experimental practice, 
A U becomes important for processes at constant volume (isochors), since then it corresponds directly 
to the heat exchange between system and environment. On the other hand, for isobaric processes the 
heat transfer corresponds to the enthalpy change AH (H = U + fiV, see Eq.(2.33)): 

isochors: 8W = 0 -> dU = 8Q = c v ■ dT (Eq.2.32) 

isobars: 8W = -pdV = —RdT -> dH = d(U + pV) = 8Q = c p ■ dT (Eq.2.33) 

For illustration, let us consider these two processes with a system of 1 mole ideal gas. Note that, at given 
heat exchange Q, the temperature effect AT will be larger in case of the isochoric process, since in case 
of the isobaric part of the energy will be lost due to volume expansion work. As a consequence, the heat 
capacity c v has to be smaller than c p , namely c p — c v = R for 1 Mol ideal gas. 

Importantly, for the ideal gas the energy change for isothermal processes is zero, i.e. in general: 

dU = 8Q + 8W = c v dT = 0 (Eq.2.34) 

This means that the volume expansion work is compensated by heat exchange, or 8Q = — 8W\ 
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For 1 mole ideal gas, this work can be calculated as: 

dW = -pdV = -ydV -> W = -f^ydV = -RT- In (|) (Eq.2.35) 

In contrast, for real gases one also has to take into account the internal work versus the attractive 
interparticle interactions, and the total internal energy change is given as: 

dU = c v dT + ndV = 0 (Eq.2.36) 

Therefore, if a real gas is expanding into vacuum, both heat exchange (Q), volume work (W) and, as a 
consequence of the first law of thermodynamics, the overall change in internal energy (A U) are all zero. 
The temperature will thus drop to compensate for the internal work (ndV). 

(Joules experiment). 


Finally, let us finish our discussion of the first law of thermodynamics in the context of ideal gas processes 
with the adiabatic process. Here, the first law will allow us a very simple derivation of an equation of state 
for adiabatic processes, which will come in handy in the next section about the 2 nd law of thermodynamics. 
For adiabatic processes of 1 Mol ideal gas the first law reads as following: 


dU = SQ + SW = 0 - pdV = c v dT = SW ad 

[jw L y i 
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To derive our equation of state, we insert the ideal gas law for p, and obtain: 
-ydV = c v dT 
— — = — dT 

V R 


(Eq.2.38) 

(Eq.2.39) 


Integration of this differential with finite boundaries (T^V) as initial and (T 2 ,V 2 ) as final state yields 


-In 




(Eq.2.40) 


Or 


v 1 ■ T l Cv/R = V 2 ■ T 2 Cv/r = const (Eq.2.41) 

If you apply the ideal gas law to replace T with p, and consider that for one mole ideal gas c p — c v = R 
then, 


p . y c p/c v — cons t (Eq.2.42) 

This is our equation of state for an adiabatic process, expansion or compression, of the ideal gas, allowing 
us to calculate any missing variable of state. 

Comparing this equation with the equation of state for the isotherm of an ideal gas (Boyles law: 
P ' V = const ), we expect at identical volume expansion from identical initial states that the adiabatic 
p(V) drops more in pressure than the isotherm, i.e. that also the temperature is lowered, since in this 
case the expansion volume work is not compensated by heat exchange. 

2.3.3 The 2 nd law of thermodynamics and the Carnot process 

Whereas the first law deals with the energy, the 2 nd law of thermodynamics is concerned with a new 
quantity of state, the entropy. We should note that in thermodynamics “order” is not just spatial order 
in the usual 3-dim space. It rather refers to a “phase space” which includes momentum (velocity) 
coordinates. For example, in an ideal gas, the spatial ordering of molecules is independent of temperature 
at constant volume. However, the “order” in velocity space decreases with increasing temperature (and 
increasing average molecular velocity in a gas) resulting in an increase of entropy. Quantitatively and 
macroscopically, entropy is defined as the reversibly exchanged amount of heat between system and 
environment at given temperature T 

dS = -> AS = (dT = 0) (Eq.2.43) 
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This means that dU=TdS - dW , and we will show later that the reversible work is the maximum you can 
get, as a theoretical limit. Eq.2.43 is the first part of the 2 nd law of thermodynamics, defining the change 
in entropy as a new fundamental quantity valid for strictly reversible processes. 

To discuss the maximum work one can extract from a reversible cyclic process, we next consider the 
famous Carnot process as a reference. This process describes an ideal reversible and therefore hypothetic 
machine partially transferring heat into work. The principle is illustrated in figure 2.7. 



Fig. 2.7.: thermodynamic machines - principle of the Carnot process 

This machine is based on the reversible heat flow from a reservoir with higher temperature T w to a 
reservoir with colder temperature T k . Note that both reservoirs are considered to be infinitely large, so 
their temperature will not change during the process. Importantly, the amount of heat transferred to the 
colder reservoir Q’ is smaller than the amount of heating taken from the hot one Q, and the difference 
W = Q - Q’ can be used a volume expansion work. Note that the overall entropy change for this reversible 
process is zero, i.e. 

AS to t = - jT + -jr = 0 (Eq.2.44) 

If we are running our machine in reverse direction, which should not change the respective amounts of 
heat and work, since the process is perfectly reversible, we obtain a so-called heat pump, which transfers 
heat from cold to warm upon work input, in the reverse direction compared to the spontaneous heat flow. 

The efficiency of this type of machine transferring heat into work is defined as the ratio of work output 
over heat input, i.e. 

w = (Eq.2.45) 

1 Q T w v H ' 
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Only if T k = 0, a limit which according to the 3 rd law of thermodynamics (see section 2.3.5.) can never be 
reached in practice, the efficiency of our machine may reach one. In general, one cannot transfer heat Q 
into work W without loss, and the highest efficiency possible is given by the reversible machine, a limit 
which can never be reached in practice. Note that, in accordance with energy conservation, it could be 
possible to transfer heat into work directly. However, the 2 nd law tells us that this cannot be achieved: 
meaning that heat somehow is a worse form of energy compared to work from a practical point of 
view. This is directly related to entropy, as we will see. First, let us derive the above formula describing 
our limit in heat-work conversion efficiency by discussion of the Carnot process in more detail. This 
hypothetical ideal circular process is based on isothermal and adiabatic expansions and compressions, 
as shown in figure 2.8. 



The first step is an isothermal expansion from A to B at the higher temperature T^ i.e. the ideal gas here 
remains in thermal equilibrium with the hotter reservoir. During this first step, the gas takes up the heat 
amount Q = Q AB from the hot reservoir to compensate the volume work. Next, the gas expands to C 
adiabatically in thermal isolation, and therefore its temperature has to drop to T k . The third step is an 
isothermal compression to D, this time in contact (or thermal equilibrium) with the colder reservoir. To 
compensate for the compression work, here the heat amount Q’ = Q CD is transferred from the gas to the 
reservoir. Finally, the circle is closed from D to A with an adiabatic compression, where the temperature 
of the gas increases again to T w . The amount of work and heat exchanged between gas and environment 
during each step are shown in the figure. 
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The amount of total work per cycle is given as 


\w\ = m I = | -RT W ■ In g£) + c„ ■ (7* - T w ) - RT k ■ In g) + c y ■ (T w - T k )\ (Eq.2.46) 


or 


I W\ = 


-RT U , ■ In 


w -~'v A 


m- 


Rr ‘- in ©i 


(Eq.2.47) 


Whereas the amount of heat needed to obtain this work, exchanged by the first step, is given as: 


101 = 


RT W 


‘"©I 


(Eq.2.48) 


Because the two adiabatic curves in figure 2.8. connect different states with identical temperature, 
respectively, the volume ratio has to be identical as well, i.e. 


v A _v B ^ v c _ 


Vd V c 


Vd V a 


(Eq.2.49) 


If we insert this equation into Eq.2.46, we obtain for the efficiency of the Carnot machine: 


11 = 


\Zw t \ 


T W~ T k 


101 


RT W - In 


© 


< 1 


(Eq.2.50) 
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Work and heat exchanged during one working cycle between the ideal gas and the environment can 
simply be expressed as definite integrals or areas in a p(V) or T(S) representation, respectively 



T 


A 



Fig. 2.9 .: p-V- and T-S-diagram of the Carnot process, areas indicating the amount of transferred energy 


Another reversible cyclic process with efficiency identical to that of the Carnot process is the Stirling 
process, which corresponds to a Carnot cycle with the two adiabatics replaced by two isochors, 
respectively. 



Fig. 2.10: p-V- circle of the Stirling process 
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Finally, the 2 nd part of the 2 nd law of thermodynamics is related to irreversible spontaneous processes. 
Note that it is not possible to prove the following statements, but, according to common experience, the 
following formulations of the 2 nd law of thermodynamics in practice are always fulfilled. 


i. For any irreversible spontaneous process the overall entropy (system + environment) has to 
increase, i.e. AS > 0. The general definition of the change in entropy at given temperature 
has been given by Clausius as AS > with AS = ^ for reversible processes, and AS > ^ 
for irreversible processes, respectively. For system plus environment, the overall heat transfer 
Q — 0, therefore in case of irreversible processes AS > 0. 

ii. Heat always flows, spontaneously, from warm to cold. 

iii. An irreversible process is always less efficient than a reversible one, meaning that heat is 
converted into work with a lower degree of efficiency. Let us consider a reversible machine 
transferring heat into work (= generator), with an efficiency of: 


ri=j = < 1 (Eq.2.51) 

V l w 

Running the machine in reverse mode, i.e. transferring work into heat flow, we have built a 
so-called heat pump with power conversion defined as: 


8 = 7 ] 


-l 


Q_ 

w 


T w 

T w~ T k 


> 1 


(Eq.2.52) 


If we now compare two given machines or two given heat pumps working with identical 
temperature reservoirs, one reversible and the other irreversible, according to the 2 nd law of 
thermodynamics 


7]rev ^ Vi 


and 


> £n 


Only for the reversible machine it is not important if it runs as a generator or as heat pump, 
in respect to the respective amounts of transferred energy. On the other hand, if we reverse 
the processes of our irreversible machine, not only the directions but also the amount of the 
transferred energies will change, and therefore e. ^ rj. _1 . 

o o irrev • irrev 


Note that these different expressions for the 2 nd law of thermodynamics related to spontaneous irreversible 
processes are all equivalent, as can be shown simply if one considers a coupled machine consisting of a 
generator and a heat pump both working between the same heat reservoirs, where the work generated 
by machine 1 is used completely to run the heat pump (machine 11). First, let us prove via this coupled 
machine, in combination with the 2 nd law of thermodynamics, that all reversible machines must have 
identical efficiency and power conversion. 
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Fig. 2.11 : Coupling of reversible generator and reversible heat pump, assuming different efficiencies 


We first assume that there exist two reversible machines with different efficiencies, and choose the 
generator as the machine with the higher efficiency, while the reversible machine with lower efficiency 
is running in reverse mode as heat pump. In conclusion, we get 


Vn = 


\w\ 

\Qu\ 


<Vi = 


\m 

i<?/i 


(Eq.2.53) 
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or l<?/l < IQ//1 and |Q/ | < |Q// | (Eq.2.54) 

This corresponds to a spontaneous flow of heat from the cold to the hot reservoir, which is physically 
not possible, or at least has never been observed. 

Next, let us consider the coupling of a reversible and an irreversible machine. There are two possibilities: 
either machine I is reversible and the heat pump II is irreversible, or the heat pump II is the reversible 
machine and the generator I the irreversible one. Let us consider the first case: 




Fig. 2.12: Coupling of two processes, one reversible the other irreversible, always leading to spontaneous heat exchange, 
i.e. heat flowing from the hot to the cold reservoir. 

The reversible machine one, if running in reverse mode, should have a higher power conversion than 
machine II. Note here again that only for the reversible machine it does not matter, in terms of amount of 
heat and work exchange, in which mode the machine is running. If you reverse the mode of an irreversible 
machine, as stated by the name, these amounts of energy transfer will change. For our coupled machine 
we now get the following net effect: 

£/ = |W| > £// = ° r ^ > anc * ^ 1 I > I* (Eq.2.55) 

This corresponds to a spontaneous irreversible heat flow from the hot reservoir to the cold one. Also, it 
should be noted that the change in entropy from the point of view of the two infinitely large reservoirs is 
given by the expression in the figure, whereas from point of view of the irreversible machine AS > Q/T. 
As a consequence, the overall change in entropy for our coupled machine AS > Q u /T k — Qu/T k — 0, 
or AS > 0, which is another synonym for irreversibility. 
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Finally, if we couple an irreversible generator I with a reversible heat pump II, 



(Eq.2.56) 


or \Qi\ > | Qu |and \Qi \ > I Qu |. 


(Eq.2.57) 


Again, the introduction of one irreversible process causes the spontaneous irreversible flow of heat from 
hot to cold and an increase in overall entropy 


Example 2.4.: 

Calculate the minimum energy necessary to cool 1 L of water within a fridge (from T = RT 
(room temperature, 21°C) to T = 4 °C). Assume that the heat capacity of water is constant in this 
temperature regime (=4.18 J/(g K)), and use the Camot-process as the best case reference. 

Solution: The amount of heat which will be transferred from the warm water to the surroundings 
(here: the fridge) without additional cooling is given as 

Q = c v ■ AT = 4.18 ■ 1000 ■ 1 7J = 71060/ 

To maintain the temperature in the fridge, this amount of heat has to be compensated following 
the working scheme of the heat pump, i.e. transferring 71060 J plus electrical work into a larger 
amount of excess heat which is transferred from the fridge into the kitchen. In the best case, the 
efficiency of the transfer of electrical energy (= work) to heat is given as 


In addition, energy conservation states that the amount of work plus compensated cooling heat 
(= 71060 J, see above) has to equal the amount of heat transferred from the fridge to the kitchen 
(Q), therefore 



71060 J+W 
W 


(17.29 — 1) -W = 71060/ => 


W = 4362 J 


To cool 1 L of water from room temperature to the 4°C within our fridge, we therefore at least 
need the following amount of electric energy: 


4362 J = 4362 W ■ s = 1.212 ■ 1(T 3 kWh. 
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2.3.4 The free energy and the free enthalpy: 

From the 2 nd law of thermodynamics and the meaning of reversibility, we realize that, to judge the 
efficiency of a given chemical process, we have to introduce two more quantities of state which combine 
energetic and entropic aspects, the free energy A (or Helmholtz energy) and the free enthalpy G (Gibbs 
enthalpy). The total differential of these new quantities is given as: 

dA = d(U — TS ) = TdS - pdV - TdS - SdT = -SdT - pdV 

dG = d(H — TS) = TdS + Vdp - TdS - SdT = -SdT + Vdp 

Note that these formulae have been derived by using the following expression for the 1 st law of 
thermodynamics: dU — TdS — pdV — dQ rev — pdV = dQ rev — dW rev . Therefore, one can show that 
A A = AU — T AS is the maximum work which you can get from a given isotherm process, which is the 
case if the process is reversible. Let us illustrate this with two quantitative examples: 


(Eq.2.58) 

(Eq.2.59) 
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i. We consider an energetically favored, but entropically unfavorable, chemical process, 
where both the energy and entropy of our system (subscript S) are decreasing (for example 
crystallization of the solute out of a saturated solution), and calculate the exchange in heat 
and work with the environment (subscript U) for the reversible and the irreversible case: 

A U s = -1000 ], AS s = -10 J/K, T = 50 K , therefore 
reversible: Q rev — T AS s — —500 J, W rev — AA — AU s — Q = —500 J 

irreversible: Qirrev < ThS s = -600 J, W irrev = AA = AU S - Q = -400 /. 

This means that in case of an irreversible process, the work you can get is always smaller than from the 
reversible process. 

ii. Next, let us consider a process which is both energetically and entropically favored (for 
example burning of fuel). 

AU S = -1000 ], AS s = 10 J/K, T = 50 K, therefore 
reversible: Q rev — TAS s — 500 J, W rev = AA = AU s — Q = —1500 J 

irreversible: Qirrev < TAS s = 400 J, W irrev = AA = AU S — Q = —1400/. 

Again, the work you can get from the process is largest for the reversible case. Note here that in our 
example the work output is much larger than the energy change of the system! 

The change in Gibbs free enthalpy has a similar meaning: it describes the maximum work you can get 
from a process at isothermal and isobaric conditions, any work due to volume change of the system 
excluded. This is derived as following: 

dG T = dH - TdS = dU + d(pV ) -TdS = dW rev + pdV + Vdp (Eq.2.60) 

i i 

dGf'P = dW rev + pdV = (dW rev — pdV ) + pdV = dW rev (Eq.2.61) 

i 

with dW rev the reversible or maximum work one can get with the exclusion of any volume work -pdV. 
For example, this meaning of AG is obvious if you consider electrochemical processes (see section 4): 

AG — —z ' F ' EMK (Eq.2.62) 

where EMK is the maximum voltage you can get if your chemical battery is running reversible, and 
— z ■ F ■ EMK is the maximum output in electric energy or work. 


44 


Download free eBooks at bookboon.com 


Basic Physical Chemistry 


Thermodynamics 



Fig. 2.1 3: difference of work (left) and heat (right), and the 2 nd law of thermodynamics. Note that, if the gas is 
expanding in one direction (= volume work), during this process a regular flow of the particles in this 
direction is needed. In contrast, if the gas is heated at constant volume, only the velocity of the irregular 
motion of the gas particles is increasing. 


Obviously, the 2 nd law of thermodynamics and all its consequences are based on a fundamental difference 
between the two ways of energy transfer, work and heat. As shown in figure 2.13, work usually requires 
some regular motion of all molecules of the system, whereas heat is just based on the average velocity 
of the particles which are moving irregularly As a consequence, heat has a “lower quality” and cannot 
be transferred into the much more regular work by 100%. 



Brain power 


By 2020, wind could provide one-tenth of our planet’s 
electricity needs. Already today, SKF’s innovative know¬ 
how is crucial to running a large proportion of the 
world’s wind turbines. 

Up to 25 7o of the generating costs relate to mainte¬ 
nance. These can be reduced dramatically thanks to our 
(^sterns for on-line condition monitoring and automatic 
luS|kation. We help make it more economical to create 
cleanSkdneaper energy out of thin air. 

By sh?fe|ig our experience, expertise, and creativity, 
industries can boost performance beyond expectations. 

Therefore we need the best employees who can 
kneet this challenge! 


Power of Knowledge Engineering 


Plug into The Power of Knowledge Engineering. 
Visit us at www.skf.com/knowledge 


45 

Download free eBooks at bookboon.com 



Basic Physical Chemistry 


Thermodynamics 


Summary of thermodynamic quantities of state and their practical meaning: 

(i) internal energy: 

The change in internal energy of a process is corresponding to the heat exchange with the 
environment at isochor conditions: A U v = c v ■ A7, 

with c v the respective heat capacity. 

(ii) enthalpy: 

The change in enthalpy of a process is corresponding to the heat exchange with the 
environment at isobar conditions: AH p = c p ■ AT , 

with c p the respective heat capacity. Note that for gases (1 mol particles), c p — c v = R. 

(iii) entropy: 

The change in entropy corresponds to the reversibly exchanges heat divided by the 
temperature, if the process is run isotherm. More general, 

AS = j^ dT 

Note that, in practice, all processes are irreversible. Therefore AS is only a limit which 
can never be measured directly, but be determined by interpolation under different 
experimental conditions approaching more and more the ideal reversible case (for 
example, chemical batteries running at very low electrolyte concentrations). 

(iv) free energy: 

The change in free energy is the maximum work output one can get from a process at 
isotherm conditions: AA T = W rev 

Again, this is a limit which can never be measured directly but only be determined by 
interpolation of a series of experimental data. 

(v) Free enthalpy: 

This is the maximum work output one can get from a process at isotherm and isobar 
conditions, any volume change effect excluded. For example, this could be 
electrochemical work: AG r p = W rev = —z ■ F ■ EMK 


2.3.5 The 3 rd law of thermodynamics 

Whereas the 1 st and 2 nd law deal with the change of quantities of state, i.e. internal energy and entropy 
the 3 rd law provides the basis to determine one absolute value: it states that the entropy of a perfect single 
crystal at T = 0 K equals zero, based on the fact that there is only one set of microscopic variables to 
describe this perfectly ordered system. 
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2.4 Heat capacities 

2.4.1 Heat capacities of gases 

We have seen in the last section that, if we try to determine changes in energy or enthalpy from 
experimentally detected changes in temperature, the respective heat capacities are needed. In experimental 
practice, these are typically determined by calibration. For pure systems, however, it is possible to predict 
the heat capacities based on physical models of the microscopic effect of increasing heat on the single 
particle motion. 

For ideal gases, we already mentioned that the difference between isochor and isobar heat capacities in 
case of 1 Mol particles is simply the gas constant R, which is easily derived as following: 

Q p = A H p = AU p + pAV = c v -AT p +R- AT p = (c v + R) ■ AT p = c p ■ AT p (Eq.2.63) 


or 


( c v + R) = c p . (Eq.2.64) 

To get an estimate for C v , we consider that a gas molecule consisting of N atoms is described by 3N 
Cartesian coordinates x ,y and z, which gives, for the microscopic motion of a molecule, 3N different 
degrees of freedom. 3 always are translation of the whole molecule in x, y and z direction, 3 more are 
rotation for non-linear molecules, whereas linear molecules only have 2 rotational degrees of freedom 
because there is no momentum of inertia if this type of molecule is rotating around its molecular main axis. 


The remaining 3N - 6 or 3N - 5 degrees of freedom, respectively, are attributed to molecular vibrations. 
Importantly, each degree of freedom of translation or rotation contributes with R/2 to the molar heat 
capacity , while a vibrational degree of freedom contributes with R. This difference can be explained, 
in a simplified way, by the different types of energy: translation and rotation are purely kinetic energy, 
whereas vibration contains both kinetic and potential energy, since it describes the motion of atoms in 
a force field. As an illustrative example, let us consider the two 3-atomic molecules C0 2 and H 2 0. 


Degrees of freedom: 


Translation: 

Rotation: 

Vibration: 


C0 2 : 

3 

2 

9-5 = 4 


Therefore, the total molar heat capacity of each species is given as: 
Cv,co 2 = 3-^ + 2-^ + 4-/? = 6.5-f? 
c V ,h 2 o = 3-| + 3-| + 3- J R = 6i? 


H 2 0: 

3 

3 

9-6 = 3 


(Eq.2.65) 

(Eq.2.66) 
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The heat capacities differ by nearly 10% just because one rotational degree of freedom is transferred to 
a vibrational one, if you go from H 2 0 to C0 2 . 

Importantly, these heat capacities are limiting values reached at high sample temperature. To address 
the vibrational degrees of freedom, you need a certain sample temperature already, typically several 
100-1000 K. Rotation is more easily addressed already at room temperature, whereas translational 
motion is for free, i.e. the 3 translational degrees of freedom are excited already at the evaporation 
temperature irrespective of its value. As a consequence, the heat capacity of the gas increases stepwise 
with temperature, and the position of the steps or characteristic excitation temperatures depend on 
specific molecular parameters such as mass of atoms, bond lengths and bond strengths. 



2.4.2 Heat capacities of solids 

The atoms within a solid cannot move nor rotate but only undergo vibrational motion. We therefore 
expect a limiting heat capacity of 3R (law of Dulong-Petit), since for a solid consisting of N atoms all 3N 
degrees of freedom correspond to vibrations, and therefore contribute with R per Mole to the heat capacity. 

Cv, molar O' -» oo) = 3R (Eq.2.67) 

However, we have learned in the last section that vibrational modes are not accessible at very low 
temperature. Therefore, experimentally the heat capacity of solids is given by the following expression 
(Debye): 


Cy,molar 



(Eq.2.68) 


with 0 D the so-called Debye-temperature that is characteristic for a given solid material. 
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The thermal accessibility of the vibrations of a solid has theoretically first been addressed by Einstein, 
who assumed that all atoms have the same oscillatory frequency According to the Boltzmann factor 
(see above, section 2.2.2) one then would expect a molar internal energy given as 


U = 3N l ■ h ■ v E ■ 


[ -p (zS) 1 


(Eq.2.69) 


And correspondingly the heat capacity is the temperature-derivative, i.e. 


C V,molar 



(Eq.2.70) 


At very high temperature, this complicated expression yields (Taylor-series-expansion of the exponential!) 


Cv, molar (T -» °°) = 3R ■ 


fh-v E \ 2 . 

h-v E 

k-T 

V k-T ) 



~ 3 R 


(Eq.2.71) 


According to Eq.2.71, the Einstein model fulfills the Dulong-Petit-law valid for every solid consisting 
of 1 Mole atoms at very high sample temperature. However, at low temperature the molar heat capacity 
predicted by Einstein is smaller than the one experimentally found. 
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The model by Debye describes the experimental result. In contrast to Einstein, he assumes that the 
vibrations within a solid cover a whole spectrum of different oscillatory frequencies, starting from zero 
to v E , with the probability (or better: degrees of freedom) of a given frequency as: 

P(v) = 9N (Eq.2.72) 

Note that the total number of degrees of freedom, given as the integral of P(v), equals 3N. This distribution 
of frequencies is plotted in figure 2.15: 



Fig. 2.15: frequency spectrum of the oscillations of a solid 

The consequence of this improved model developed by Debye is that oscillations of lower frequencies 
are already accessible at lower temperatures. Therefore, the heat capacity at lower temperatures is larger 
than predicted by Einstein, and in agreement with the experimental results. Physically, the existence of 
lower frequencies is feasible if we consider the coupling of multiple atoms to one oscillator. The frequency 
of a harmonic oscillator depends on the mass which is moving, and the force which is constraining this 
movement. If you consider the coupling of atoms to an oscillating aggregate, the mass scales with the 
total number of atoms whereas the force only scales with the number of surface atoms. Consequently, 
the oscillatory frequency will decrease with increasing number of coupled atoms, vibration of the whole 
solid body showing the smallest frequency possible. 

2.5 Phase equilibrium 

2.5.1 Pure components, p-T-diagrams 

The Gibbs phase rule states that, in any given system in equilibrium, the maximum number of intensive 
variables (= variables which do not depend on the size of the system, i.e. temperature, pressure) which can 
be freely chosen without altering the state of the system is given as F = K - P + 2, with K the number of 
chemical components and P the number of phases in coexistence. F is also called “degrees of freedom”. For 
example, for an ideal gas we get P =1 and K =1, therefore F = 2, meaning you can chose temperature T and 
pressure p independently without changing the state of the system, i.e. ideal gas. On the other hand, for 
real matter at its triple point we get K =1 and P =3, therefore F = 0, i.e. no intensive variable can be chosen. 
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The p-T-diagram of a real substance is shown in figure 2.16. It consists (at least!) of three different regions 
for the different phases solid, liquid and vapor, and three different equilibrium curves which pair-wise 
separate these regions. All three curves meet in one point, the triple point, where all three phases are 
in equilibrium. 



Let us start with the liquid-gas-equilibrium, where F = K - P = 1, meaning we can either chose the 
pressure or the temperature independently. To obtain the relation between these two variables p(T), 
the vapor pressure curve, we start with the equilibrium condition that the free Gibbs enthalpies of the 
coexisting phases (1 = liquid, 2 = vapor) have to be equal, choosing (arbitrarily) 1 mole substance for 
each phase. Here, one should note that the equilibrium conditions p, T will not depend on the amount 
of matter considered in each phase. E.g., at normal pressure water always boils at 100°C irrespective of 
the amount: 


G 1 = G 2 (Eq.2.73) 

Note that Eq.(2.73) is directly related to the 2 nd law of thermodynamics: for spontaneous processes or 
non-equilibrium AS > y. This yields, in case of isobaric and isothermal conditions or at given p, T that 
AS > or AG = AH - TAS < 0, i.e. the condition of stationary equilibrium. If we follow the vapor 
pressure curve shown in figure 2.16, the respective free enthalpies of the two coexisting phases will 
change with temperature and pressure. However, since at any given point of the curve the stationary 
phase equilibrium is maintained, these respective changes in free enthalpy have to be equal, leading to: 

dG 1 = —SidT + Vidp = dG 2 = —S 2 dT + V 2 dp (Eq.2.74) 
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Taking into account that the molar volume of the vapor phase exceeds that of the liquid phase by 3 
orders of magnitude, and inserting the ideal gas equation for V 2 , we obtain, after separating the variables: 


dp= S -^dT = 


T^-T-S, A V H A V H dT 

T-V 2 T-V 2 R T 2 P 


(Eq.2.75) 


Integration from the triple point coordinates (7 0 , p 0 ) to an arbitrary point on the vapor pressure curve 
leads to the Clausius-Clapeyron-equation: 



A yH ( 1_ _ r\ 
R \T 0 t) 


(Eq.2.76) 


Alternatively, any point on the vapor pressure curve (T 0 , p 0 ) can be chosen as reference, for instance the 
boiling temperature at standard pressure 1 atm. The Clausius-Clapeyron-equation allows one to predict 
the boiling behavior of a given pure substance. Note that boiling, corresponding to crossing the vapor- 
pressure curve from liquid to vapor regime, can be achieved in two different ways: either increasing the 
temperature at given pressure, or lowering the pressure at a given temperature (vacuum distillation). 
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Sublimation, the phase transition from solid to vapor, can be described exactly the same way as the 
liquid-vapor transition, since again the molar volume of the vapor phase exceeds that of the condensed 
phase by several orders of magnitude. We simply consider 1 to be the solid phase, replace the evaporation 
enthalpy by the sublimation enthalpy, and get in analogy to the Clausius-Clapeyron-equation: 

Ht)=*¥■■(?-%) 2 - 77 > 

Note that, if our reference point is the triple point, this time we chose the variable point as the lower 
integration boundary, since the sublimation curve runs from the triple point to the left in the p-T-diagram. 

Finally, let us consider the solid-liquid equilibrium. Here, we have to consider both the volume of the 
liquid and the solid phase. At not too high pressure, these volumes will be constant due to the very 
low compressibility of a condensed phase, and the equilibrium condition dGi = dG 2 (here: 1 = solid, 
2 = liquid) after separation of the variables leads to: 

dp = w'T (Eq.2.78) 


or: 


p = p» + tf |n © <E< i' 2 ' 79) 

Again, we chose the triple point as the reference ( T 0 , p 0 ) since it is common to all three phase equilibria 
curves. Note that the variable point ( T, p) on the melting curve, as shown in the figure, has to be above 
the triple point or at higher pressure. However, if the molar volume of the liquid phase is smaller than 
that of the solid phase, as for example in case of water, the melting temperature T is smaller than 7 0 
(anomalous behavior). As a consequence, one can melt water by applying pressure (see ice skating), or 
fortunately for us (and for the fish) a lake or river freezes from the top, thereby partially isolating itself 
from the cold air of the environment and rarely freezing completely. 


We can simplify the melting curve expression if we consider that the melting temperature changes only 
slightly with increasing pressure, and using a Taylor series expansion for the logarithm: 


. A S H . f T\ , A S H , ( A , T-T 0 \ , A S H 

P = P‘ + U-' n k) = P‘ + IJ-' n ( 1 + —)”P« + U 


T-Tp 

To 


(Eq.2.80) 


This is a linear expression for p(T). Note that the slope of the melting line depends on A S V which is 
typically positive (volume of the liquid phase is larger than that of the solid one), but can be negative 
(anomalous behavior, for example water!). 
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2.5.2 Partial molar quantities, the chemical potential 

Before we discuss the phase behavior of multi-component systems, for example solutions or mixtures, 
we have to introduce a new type of physical-chemical quantities, partial molar quantities. We will mainly 
use the partial molar free enthalpy or chemical potential, but the difference between pure components 
characterized by molar quantities, and mixture characterized by their partial molar correspondents, is 
best illustrated via the molar volume. 

(i) Partial molar quantities, the partial molar volume 

In the ideal case, the total volume of a binary mixture is given by the respective sum of the molar volumes 
of the pure components times the mole number of each component, respectively: 

V = n 1 -V 1 ° + n 2 - V 2 ° (Eq.2.81) 

However, in most cases specific interactions between the two types of molecules lead to a volume change 
after mixing, which can be expressed by an excess term as: 

V = n 1 -V 1 ° + n 2 - V 2 ° + AV mix = n 1 -V 1 +n 2 -V 2 (Eq.2.82) 

Here, the excess term corresponds to the total difference of the molar volumes of the pure components and 
the so-called partial molar volumes, which have to be used to describe the total volume of a real mixture: 

w mix = m • (V, - V ! 0 ) + n 2 • (V 2 - V 2 °) (Eq.2.83) 

The practical meaning of these partial molar volumes is illustrated by the following fictional experiment: 
take a given binary mixture to which we add an infinitesimally small amount of the two components 1 
and 2, respectively. The corresponding change in volume then is given as 

dV = (£ 7 ) dni + (£ 7 ) dn 2 = ^i dn i + v 2 dn 2 (Eq.2.84) 

The partial molar volumes therefore correspond to the partial derivatives of the total volume over the 
molar amount of the respective component. Next, we consider a macroscopic volume consisting of 
many of these infinitesimal volume changes dV , where the composition cZn 1? dn 2 is identical for all dV: 

V = £ dV = V 1 £ dn x + V 2 Z dn 2 =V 1 n 1 + V 2 n 2 (Eq.2.85) 

Here, it should be stressed that the partial molar volumes only depend on the relative composition of 
the mixture, and not on the absolute amounts of the components! 
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Taking the differential of this expression we get: 

dV = d(y r n x + V 2 n 2 ) = dV x n x + V 1 dn x + dV 2 n 2 + V 2 dn 2 
On the other hand, we have shown before that 

dV = dnx + d n 2 = V 1 dn 1 + V 2 dn 2 

Combining these two expressions for dV , we obtain the important Gibbs-Duhem-equation 
dV x n x + dV 2 n 2 = 0, or n 1 dV 1 = —n 2 dV 2 


(Eq.2.86) 


(Eq.2.87) 


(Eq.2.88) 


This equation enables us to calculate the dependence of the partial molar volume of one component 
on the composition of a binary mixture, if we measure this dependence for the other component. To 
measure the partial molar volume of component 2, for example, one has to measure the total volume of 
the mixture while only varying the amount of component 2, keeping the other variables n 1? pressure 
and temperature constant. The partial molar volume V 2 then is obtained as the tangential slope of a plot 
of total volume V versus n 2 , i.e. 


V 2 = 

Van 2 J p> T >ni 


(Eq.2.89) 
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Note that this slope could be negative, as for example in case of pure water to which we add NaOH: here, 
the formation of a more densely packed hydration shell leads to an overall contraction of the volume if 
the amount of added NaOH is still small. 

(ii) The chemical potential 

Fundamentally much more important than the partial molar volume is the partial molar free enthalpy, 
which is also called the chemical potential, and provides us with the criterion for phase equilibrium of 
multi-component systems: 



(Eq.2.90) 


Note that for a pure component ( n 2 = 0), this chemical potential is simply the molar free enthalpy 
The chemical potential is the key quantity in the discussion of both phase equilibrium and chemical 
equilibrium of mixtures, like the free enthalpy for the pure component phase equilibrium as we have 
shown before. The equilibrium conditions for mixtures are^ = i u 1 (stationary) and d^ , 

respectively. The expression for the dependence of the chemical potential on the sample composition of 
a binary mixture can easily be derived if we consider the mixing process of two different ideal gases A, 


B. Note that therefore we will ignore any enthalpic contributions or specific interactions, and therefore 


the difference between molar free enthalpy G and partial molar free enthalpy [jl will only depend on 


the mixing entropy! 

(iii) Mixing free enthalpy of two ideal gases and the chemical potential 

Before we will treat the phase equilibria of multicomponent systems or mixtures, we will derive an 
expression for the chemical potential by considering the following mixing process of two different ideal 
gases, as sketched in figure 2.17. 


n A , V A , p 

nB, Vb, p 


n A +n B , V a +V b , p=p A +p B 


Fig. 2.17: Spontaneous mixing of two different ideal gases to deduce the free enthalpy of mixing 

Before mixing, both gases with molar amount n A and n B , respectively, shall have the same pressure p. 
The chemical potential, or molar free enthalpy, of gas A before mixing then is given as: 


A Ia*(p) = A*/’ 0 + Jp P 0 vdp = p/' 0 + fp 0 ydp = G Amolar (p) 


(Eq.2.91) 
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Here, the symbol * signifies pure component, and signifies standard conditions, i.e. here it refers to the 
standard pressure ofp = 1 bar. The total free enthalpy of our system before mixing (figure 2.17, left) is 
therefore given as 

G = n A fi A * + n B n B * = n A (u /' 0 + RT ■ ln^) + n B (ji B *’ 0 + RT ■ ln^) (Eq.2.92) 

After mixing, the two gases assume their respective partial pressures given by the molar fractions, i.e. 

v ‘ = -^' v = x ‘' p • p ‘ = ^r B ' p = x ‘- p < E i- 2 - 93 > 

The total free enthalpy after mixing therefore is given as 

G = n A (m /’ 0 + RT ■ \nfy + n B + RT • In^f) (Eq.2.94) 

and the free enthalpy of mixing is 

A m G - G — G - n A RT ln^- + n B RT\n^~ = nRT(x A \nx A + x B lnx B ) (Eq.2.95) 

Note that the expression in brackets is negative. Therefore the free enthalpy of mixing is always < 0, 
which is obvious since the mixing of two ideal gases is an irreversible spontaneous process. 

Since we have considered here only entropic effects, the enthalpy of mixing has to be zero, and the 
mixing entropy is simply given as: 

= —nR{x A In X 4 + x B lnx B ) (Eq.2.96) 

This is always > 0 , which is obvious for an irreversible spontaneous process according to the 2 nd 
fundamental principle of thermodynamics. 

In the following, we will use this general expression for the chemical potential of an ideal purely entropic 
mixture irrespective of its state (gas, liquid or solid): 

p A = Pa* + RT\nx A (Eq.2.97) 

Note again that /^refers to the pure component at identical pressure and temperature as the mixture, 
and therefore is the molar free enthalpy G A mo i ar . 
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2.5.3 Phase behavior of binary systems - liquid-vapor-transition 

We start with a liquid mixture containing only one evaporative component, the solvent, and will derive 
a formula for the increase in boiling temperature in respect to that of the pure solvent. This system is 
formally sketched in figure 2.18. 



Fig. 2.18: Phase equilibrium between a binary liquid mixture 
and the pure vapor of the volatile component A 


At equilibrium, the chemical potentials of component A, which coexists in both the vapor and the liquid 
phase, have to be identical, respectively, or 


Ha = Ha 


(Eq.2.98) 
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Since we are interested in the dependence of boiling temperature on concentration of component B in 
the liquid phase, we have to use the derivative, or: 

dn' A = dii' A (Eq.2.99) 

To avoid differentiation of the expression RT In x A which is part of the chemical potential of the binary 
liquid phase, we consider alternatively 

d^- = d^- (Eq.2.100) 

d^r + R • d\nx A = d^r (Eq. 2 . 101 ) 

Note that the pressure in this case is kept constant, the only intensive variable being the boiling 
temperature T and the concentration or molar fraction of the solvent in the liquid phase x A . 


To solve this differential equation and thereby obtain a quantitative relation between phase transition 

Q 

temperature and solute concentration, we need the total differential of - , which is given as: 

d(j) = -^ dT + ^ d P (Eq.2102) 

Noting that n' A *= G Amolar = G a *, we obtain 

-^dT + R-d\nx' A = -^j^dT (Eq.2.103) 

- 7r) dT = ^dT = -R ■ d lnx' A (Eq.2.104) 

with A v H a the molar evaporation enthalpy of component A. Integration with the boundaries pure solvent 
(T = T a * , x A = 1, p = 1 bar ) and an arbitrary mixture (7, x A , p = 1 bar), we finally get: 

— {~t + T7* j = - ln ^ (Eq.2.105) 

Note that T > T A *, i.e. the mixture has a higher boiling temperature than the pure solvent. Also, this 
effect is purely of entropic origin, since our derivation is based on the chemical potential using the formula 
we previously have deduced from the mixing free enthalpy of ideal gases. Therefore, the effect also does 
not depend on the chemical nature of the solute component B, but only on the relative concentration, 
or better, the molar fraction of the solvent in the liquid phase x A . 
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This type of entropy-based effect found in binary mixtures is also called a colligative phenomenon, 
further examples are lowering of the freezing temperature (see next chapter), lowering of the vapor 
pressure which is the direct correspondent of the increase in boiling temperature (law of Raoult), and 
finally osmotic pressure. All respective mathematical expressions are derived from the same principle, 

iii 

i.e. the equality of the chemical potential of the solvent component in both coexisting phases [i A = Ta- 
This leads us to an alternative derivation of the increase of boiling temperature formula merely based 

in 

on this stationary phase equilibrium \i A — Ta and avoiding the explicit mathematical solution of any 
differential equation: 

We start with the stationary phase equilibrium given as: 

liA * + RT ■ 1 nx' A = h'a* (Eq.2.106) 

and therefore, using the Gibbs-Helmholtz-expression A V G A = A V H A — TA v S a : 

Ta * — ft a * = A v G a — AyH A — TAyS A — RT ■ lnx^ (Eq.2.107) 

with A v H a the molar evaporation enthalpy of the pure solvent, and A v S a the molar evaporation 
entropy The evaporation entropy, defined as the reversible heat uptake normalized by the evaporation 
temperature, can simply be replaced if we consider this equation at the boundary case of the pure solvent 
(T = T a * , x A = 1, p = 1 bar ): 


A v H a - T a *■ A V S A = RT ■ In 1 = 0, or A V S A = 


(Eq.2.108) 


and we immediately obtain our above result 


by H a 
R 


(-? + ?^)=- in ^ 


(Eq.2.109) 


Note that this expression can further be simplified if we consider very dilute solutions, i.e. comparatively 
small values of x B = 1 — x A , and use the Taylor series expansion of the logarithm: 

+ = * X B (Eq.2.110) 


Also the left hand side of this equation can be simplified considering that the change in boiling temperature 
AT — T — T a *will be comparatively small for such dilute systems: 


*v Ha (_ I , M = *yH A ( T-T a * \ _ A v H a ( AT \ = / 
R V T^~T a *J R \T a *-T ) ~ R \T a * 2 J B 


(Eq.2.111) 
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Finally, we will get the following simple formula for the colligative phenomenon of boiling temperature 
increase: 


, m B 

A v H a ( AT \ _ > _ Mb_ _rnB_ m MA_ 

R \T a * 2 ) B n A ^4 M b m A 

M A 


Or 


(Eq.2.112) 


. = RJ\a^Ma 

A V H A 


m B . 

M b m' A 


(Eq.2.113) 


Here, n B is the molar amount of solute in the mixture, m B the solute mass and M B the molar mass of 
the solute, index A accordingly for the solvent component in the mixture. The effect AT is therefore, in 
case of dilute solutions, directly proportional to the molal concentration of the solute B , and it does not 
depend on the chemical nature of the solute but only on that of the solvent, with the proportionality 
factor, also called ebulioscopic constant, given as R Ma . 
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Next, we consider the liquid-vapor-equilibrium of a binary mixture where both components A and B 
are volatile, as sketched below: 



Fig. 2.19: Phase equilibrium between a binary liquid mixture 
and the vapor mixture in case both components are volatile 


The vapor pressure of the liquid phase then is given by Raoults law, i.e. summing up the reduced (in 
respect to the pure liquid) vapor pressures of both components A and B 

p = x'a'Pa* +x b ■ p B *= x A -p A * +(1 - x A ) ■ p B *= p B * +(p A * -p B *) ■ x A (Eq.2.114) 

This linear equation describes the boiling line, with the two intercepts Pa * and Pb * the vapor pressures 
of the pure liquid, respectively. 

For the composition of the vapor phase we use Daltons law, which describes the partial pressure of 
component A as: 


Pa ~ x a ' P (Eq.2.115) 

At equilibrium for two coexisting phases, the overall pressure as well as the pressure per component 
must be identical for liquid and vapor state. Therefore, we simply may insert the linear equation for p 
into Daltons law, and obtain either a relation for the composition of the vapor phase in respect to that 
of the liquid phase (= coexistence curve), or for the gas pressure as a function of vapor composition 
(condensation curve): 


x A = 


x'a'Pa* 

Pb*Hpa*~Pb*>x a 


(Eq.2.116) 


p = 


Pa*'Pb* 

Pa*Hpb*-Pa*>x'a 


(Eq.2.117) 
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Figure 2.20: Coexistence curves - compositions of liquid and corresponding vapor phase 

The figure shows some examples: if both components have the same volatility, i.e. Pa *— Pb *> the 
composition of the liquid and the vapor phase is identical, and distilling leads to no fractionation of 
the liquid mixture. On the other hand, if one component has a higher volatility, e.g. Pa *> Pb the 

n r 

more volatile component A becomes enriched in the vapor phase upon vacuum distillation ((x A > x A ). 

If you combine the boiling line and the condensation curve (Eqs.2.114 and 2.117), you get the isotherm 
phase diagram of boiling for an ideal, i.e. purely entropic, binary mixture. 



Fig. 2.21: Isothermal boiling curve of an ideal binary mixture, and the principle of distillation 
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The average composition of the whole system comprising both liquid and vapor is given as 


_ N''X A +N"'x" A 

^4 AT 


(Eq.2.118) 


with Nthe total amount of particles (both A and B!) in the liquid phase, and N" correspondingly the 
total amount of particles in the vapor phase. Boiling line and condensation curve enclose the two-phase- 
region, where liquid and vapor phase coexist, with respective compositions given by the horizontal 
intersection with the boiling line or condensation curve, respectively. 


Finally, the arrows inserted in figure 2.21 illustrate how a vacuum distillation proceeds: if the pressure 
becomes smaller than the limit defined by the boiling line, we get the formation of vapor enriched in the 
component A of higher volatility. Further lowering the pressure leads to a decrease in component A in the 
liquid phase, but, if all vapor is collected in the same flask, also a decrease in A in the vapor phase compared 
to the first vapor formed, until all liquid has been evaporated and the overall vapor composition corresponds 
to the composition of the original liquid mixture. This also illustrates the obvious fact that, if distillation is 
used for purification of a liquid, one has to collect the vapor in different flasks successively (fractionation). 
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For fractionated distillation, you can either lower the pressure, in which case you need the isothermal 
boiling phase diagram to predict the process, or you can increase the temperature. This leads us to the 
isobar boiling phase diagram sketched in figure 2.22. The corresponding equations are derived from 
the stationary equilibrium condition fi A = ji A , analogous to the colligative phenomenon of increase in 
boiling temperature derived before (see Eq.2.98). 



The only difference is that this time both the vapor and the liquid phase are a mixed phase, therefore: 


[i A * +RT -\nx A = ii A * +RT ■ In x A 


(Eq.2.119) 


or 

li A * —\i A *= A v G a = A v H a — T ■ A v S a = RT -\nx A — RT - \nx A =RT ■ ln^r (Eq.2.120) 

x A 


Using again the pure liquid A as boundary condition, we can replace the entropy: 


AyH A T a *■ AyS A — RT A *■ Iny — 0 , AyS A 


A vHa 
t a* 


(Eq.2.121) 


And finally we get: 

X A _ Ay Ha _ T—T A * 

x' A R T-T a * 

And analogous for the 2 nd component B: 
j n 2A — . t ~ t b* 

X B R T'Tb* 


(Eq.2.122) 


(Eq.2.123) 
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Inserting x B = 1 — x A , x B = 1 — x A , and replacing , we finally obtain the following relative 
complicated expression for the isobar boiling curve: 



(Eq.2.124) 


Importantly, T A * and T B * in this formula are the boiling temperatures of the respective pure 
components at the pressure of the experiment, which is not necessarily the standard pressure 1 bar! 


Note that we can derive this expression in a much simpler way if we insert the Clausius-Clapeyron- 
equation, connecting vapor pressure and boiling temperature of the pure components, i.e.: 


f px*m \ 

i _ h v H x _ i 

( i 1\ 

\Px*(T x *)) 

R ' 

^ 1 
* 

* 1 

1 

HI 


into the isothermal boiling curve, with p = Vx * (Tx *) the actual lab oratory pressure or boiling pressure 
(i.e.: at T = T x * we have P — Px * and therefore the pure liquid X would be boiling) 


x A = 


P~PB* 

Pa*~Pb* 


(Eq.2.126) 


Finally, we consider a non-ideal mixture, where the interparticle attraction between molecules of 
component A and B is weaker than the intermolecular attraction between molecules of the same species 
(A or B). In this case we find a minimum in the boiling curve (see figure 2.23). 



Fig. 2 .23: Isobaric boiling curve of an azeotropic binary mixture (azeotropic minimum, e.g. water/ethanol) 


66 


Download free eBooks at bookboon.com 


















Basic Physical Chemistry 


Thermodynamics 


This type of interactions destabilizes the liquid phase, and therefore shifts the boiling temperature 
towards lower values in comparison to the boiling point of pure components A or B, respectively. As a 
consequence, the boiling curve intersects with the condensation curve at the minimum, and a knot is 
formed in the two-phase-region, also called azeotropic point. At this composition, the vapor upon boiling 
shows the identical composition as the liquid phase, and no purification by distillation is possible. An 
example of this type of maximum-azeotrop is the mixture H 2 0/EtOH. 

If the interactions between molecules of type A and B are stronger than the attractive interactions A-A 
or B-B, the azeotrop is a maximum in the condensation curve, and the boiling curve forms the knot. 
An example of this type of azeotrop is a mixture of H 2 0 and HN0 3 . 
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2.5.4 Phase behavior of binary systems - solid-liquid-transition 

We start again with the simplest case: a binary mixture of components A and B in equilibrium with a 
pure solid of component A. 



Fig. 2.24: Phase equilibrium between a binary liquid mixture 
A, B, and the pure solid of the component (A) 


As in case of the liquid-vapor equilibrium discussed before, the chemical potentials of component A, 
which coexists in both the liquid and the solid phase, have to be identical, respectively, or 


V-A = V-'a (Eq.2.127) 

Since we are interested in the dependence of melting temperature on concentration of component B in 
the liquid phase, we have to use the derivative, or: 

d\i A = d\i A (Eq.2.128) 


As before in case of the increase in boiling temperature of a binary liquid mixture, we get at constant 
pressure the relation: 



RT 2 _ RT 2 
H A*~ H A * h s H A 


(Eq.2.129) 


with A s H a the molar melting enthalpy of component A. Integration of this differential equation leads to 


*sH A 

R 


(-f + fh) = ln ^ 


(Eq.2.130) 


Note that this expression again can further be simplified if we consider very dilute solutions, i.e. 

r r 

comparatively small values of x B = 1 — x A , and use the Taylor series expansion of the logarithm: 

(Eq.2.131) 


68 


Download free eBooks at bookboon.com 









Basic Physical Chemistry 


Thermodynamics 


Also the left hand side of this equation can be simplified considering that the change in melting 
temperature A T = T — T A * will be comparatively small for such dilute systems: 


Wa ( 1 ■ 1 \ = a s h a ( t-t a * \ _ a v h a /-a t\ _ 

R V T^Ta*) R \T a *-T ) ~ R \T a * 2 ) 


(Eq.2.132) 


Finally, we will get the following simple formula for the colligative phenomenon of melting temperature 
decrease: 


, m B 

A s H a ( AT \ _ n B _ m b _m B M a 

R \T a * 2 ) B n A ^a_ M b m a 

M A 

or 


(Eq.2.133) 


_ R-T A * 2 -M A _ rn B 1 
A s H a Mb m A 


(Eq.2.134) 


Here, n B is the molar amount of solute in the mixture, m B the solute mass and M B the molar mass of 
the solute, index A accordingly for the solvent component in the mixture. The effect AT is therefore, in 
case of dilute solutions, directly proportional to the molal concentration of the solute B , and it does not 
depend on the chemical nature of the solute but only on that of the solvent, with the proportionality 
factor, also called cryoscopic constant, given as R Ma . 


Consider the following experiment to measure the melting temperature: a beaker containing the liquid 
sample, either pure solvent or solution, is embedded with a cooling bath, say a water/ice-mixture 
including salt. In case the cooling process is very slow, the sample freezes at equilibrium conditions, and 
the following cooling curves are obtained for pure solvent and solution, respectively: 




Fig. 2.24: Experimental cooling curves of pure liquid (left) and a binary mixture containing only one crystallizable 
component (right) 
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In case of the pure solvent, first the liquid is cooled until its freezing temperature is reached (1). At T = T A *, 
the solid phase is formed, and the temperature is kept constant until all liquid is frozen (2). Formally, 
at this stage the cooling is compensated by the heat of freezing, leading to a constant temperature. If all 
liquid is frozen, the pure solid is further cooled (3). 

In case of the solution, first the freezing temperature T A is lower than in case of the pure solvent due 
to the colligative phenomenon of lowering of the freezing temperature. In addition, the temperature is 
further decreasing once the freezing has started, since freezing of the solvent causes a further increase 
of the dissolved solute in the liquid phase, enhancing the colligative effect (2). Finally, once all solvent 
has frozen, a solid phase is further cooling (3). 

Next, let us consider the case where both components A, B are freezing, but not forming a mixed 
crystal in the solid phase. In this case, the phase behavior is a combination of the system we have just 
discussed, i.e. freezing of only the solvent, where the role of the solvent is played either by component 
A or B, respectively. 
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Fig. 2.25: Phase equilibrium between a binary liquid mixture and the solid phase in case both components may crystallize but do 
not form a mixed crystal 



0 x A ‘,x A “ 

Figure 2.26: Isobar melting diagram of a binary mixture with Eutectic point 
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At small volume fractions of component B (left side of the diagram), A plays the role of the solvent: at 
x B = 0 the transition from liquid to solid takes place at the freezing temperature of the pure component 
A (T = T a *), and it is decreasing with increasing x B due to the increase in the colligative phenomenon 
of freezing temperature suppression. In this region, only component A is found in the solid phase! On 
the other hand, at large volume fraction x B component B plays the role of the solvent: at x B = 1 the 
freezing temperature is that of pure component B (T = T B *), decreasing with decreasing x B (or increasing 
x A = 1 - x B ) again following the increase of the colligative freezing point suppression. 

Both freezing point suppression curves intersect at an intermediate concentration (in the figure at 
x B = 0.6), the Eutectic point. Formally, this point correspond to the azeotropic knot in the phase diagram 
of boiling of a binary liquid mixture: at this temperature, both components A and B are freezing, but not 
mixing within the solid phase on a molecular level! At this point, therefore, the liquid mixture cannot 
be purified by crystallization. 

Another aspect is the number of thermodynamic degrees of freedom, or independent intensive variables, 
also given in the figure: in the region of the pure liquid, F = 3, namely pressure, temperature and 
composition x B . In the two-phase regions, F is only 2: at given temperature, the compositions of the liquid 
and the solid phase are defined, or, at given composition of the liquid phase, the melting temperature is 
defined. Therefore, x B and T cannot be chosen independently without leaving the state of two coexisting 
phases. Note that in this region the composition of the solid phase is either simply x B = 1 or x B = 0. 
Finally, in the pure solid region, F also is only 2, since the composition of the two coexisting solid phases 
is x B = 1 and x B = 0, respectively. 

This melting diagram leads to the following experimental cooling curve: 



Figure 2.27: Experimental cooling curve for a binary mixture with Eutectic point 
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In region (1), the pure liquid is cooling until at T = T A the pure component A (x B = 0) starts to freeze. 
In region (2), component A continues to freeze, leading to an increase of the solute component B in the 
liquid phase and consequently a further decrease of the freezing temperature, until, at the eutectic point, 
also component B starts to freeze. In region (3), both components A and B freeze and the composition 
of the liquid phase is kept constant at that of the eutectic point, wherefore in this region also the freezing 
temperature remains constant, until all liquid has frozen. Finally, in region (4), the binary solid mixture 
of pure components A and B is further cooling. 

We close this section on the liquid-solid-phase behavior of binary systems by showing some typical 
phase diagrams: if both components form a mixed crystal in the solid phase, the isobaric melting phase 
diagram looks identical to the isobaric boiling phase diagram. The formal treatment via the chemical 
equilibrium then is exactly identical, replacing boiling temperatures and boiling enthalpies with freezing 
temperatures and freezing enthalpies, respectively. 
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On the other hand, if the two components are only partially miscible in the solid phase, various different 
types of phase diagrams are found, as shown in figure 2.29: 



Figure 2.29: Experimental isobar melt diagram of Ag/Pt (from: Gerd Wedler und Hans-Joachim Freund, 
Lehrbuch der Physikalische Chemie, p.396, 6.Auflage, Weinheim 2012. Copyright Wiley-VCH Verlag GmbH 
& Co. KGaA. Reproduced with permission.) 


74 


Download free eBooks at bookboon.com 














Basic Physical Chemistry 


Thermodynamics 


2.5.5 Colligative phenomena 

In this section on colligative phenomena, we will first combine the effects of increase in boiling 
temperature and decrease in freezing temperature by discussing the chemical potential in dependence 
of temperature for pure liquids and dilute solutions, respectively. 



Fig. 2.30: The temperature-dependence of the chemical potential of a pure solvent and of a dilute 
solution, explaining the colligative phenomena of "decrease in melting temperature" and "increase in 
boiling temperature" 

If we regard colligative phenomena, both the solid and the vapor phase are no mixtures in the 
thermodynamic sense, and therefore the chemical potential of the pure solvent and the dilute solution 
only differ at the liquid state by the term RT ■ In x A . Therefore, the liquid line is shifted to lower chemical 
potential for the liquid solution. 

Note that the slope dp/dT is given by the entropy of the respective phase, since dG = —SdT + Vdp. 
The disorder increases from solid to vapor, and consequently the slope becomes steeper from left to right. 
This explains the colligative phenomena also in a simple graphical way: for the solution, the intersections 
with the solid line and the vapor line are, in respect to the pure solvent, shifted to the left and right, 
respectively. Further, this shift has to be more pronounced for the melting transition, explaining why 
you often find an experiment about freezing point suppression in the basic practical course physical 
chemistry, but rarely an experiment about the increase of the boiling temperature: the latter effect is 
simply too small for dilute solutions to be measured in a practical course. 
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Example 2.5: 

Consider two aqueous solutions containing 1 wt% of salt (NaCl) or sugar (saccharose), 
respectively, and calculate the boiling and melting temperatures using the approximation for 
very dilute solutions. The kryoscopic constant of water is K c = 1.86 K ■ kg /mole, the 
ebuliscopic constant is K e = 0.521 K ■ kg/mole. 

Solution: the respective colligative phenomena are given as 

71 

AT — K x -, with m s the mass of the solvent water (= 1 kg), and n the molar amount of 

m s 

dissolved particles. With the molar masses of salt (M(NaCl) = 58.5 g/Mol, note that 2 ions are 
formed per NaCl upon dissolution in water) and sugar (M(saccharose) = 342.3 g/mole), we 
calculate the following: 

decrease in melting temperatures (compared to T = 0 °C): 

2 10 10 

AT = 1.86 ■ = - 0.636 K for salt, AT = 1.86 ■ -p 1 = - 0.054 K for sugar, 

increase in boiling temperatures (compared to T = 100 °C): 

2 10 10 

AT = 0.521 ■ -p = + 0.178 K for salt, AT = 0.521 ■ = + 0.015 K for sugar. 

Note that the effects are much stronger for the salt solution, and also that the change in melting 
temperature is nearly 4 times stronger than the change in boiling temperature. Typically, the 
freezing temperature decrease is much larger than the boiling temperature increase, which is 
one explanation why typically only freezing point suppression is quantitatively investigated in 
a Physical-Chemical lab course. 


Finally, we are still missing one last colligative phenomenon, which also is the most pronounced and 
therefore the most suitable for experimental practice, i.e. accurate measurement of the solute molar 
fraction: the osmotic pressure. 
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Fig. 2.31: Osmosis, pure solvent in the left chamber (I), dilute solution in the right chamber (II) 


As sketched in figure 2.31, at the beginning of our experiment solvent and dilute solution are separated 
by a semipermeable membrane, which only allows permeation of the small solvent molecules. Also, 
both liquid phases shall have the same hydrostatic pressure. To reach equilibrium, the system tries to 
equilibrate to solute concentration in both cylinders, which can only be achieved by migration of solvent 
from left to right through the membrane. This, however, causes an increase of the hydrostatic pressure 
in the cylinder containing the solute, and a corresponding decrease in pressure in the cylinder with the 
pure solvent. Finally, the system reaches chemical equilibrium at a defined pressure difference n which 
depends on the actual solute concentration. 
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Note that, in experimental practice, the amount of solvent flowing through the membrane should be 
negligibly small in comparison to the overall amount of solvent in the solution in the right chamber. 
Otherwise, the starting concentration of the solute would change due to dilution with the inflowing 
solvent, and one would have to recalculate the solute concentration in the right chamber at stationary 
conditions accordingly. In experimental practice, one is therefore often using bulky chambers with 
very thin necks, so the pressure difference is only created by a small amount of inflowing solvent. More 
sophisticated, one can also measure an external force per area needed to completely suppress any solvent 
flow from left to right if solution and solvent are otherwise kept at identical pressure. In this case, the 
osmotic pressure is compensated by this external force. 

Formally, the osmotic equilibrium is treated as following: 

Pa(.x a = l,p)=fi A * (p) = I*a(.Xa < l,p + JT) (Eq.2.135) 


This means that the chemical potential of the pure solvent at lower pressure is identical to the chemical 
potential of the solvent in a solution at higher pressure p + n. Explicitly inserting the dependence of 
chemical potential on molar fraction of the solvent ( p A ( x A < 1, p + n) — p A * (p + n) + RT ■ In x A ) 
and pressure ( = v A,mol dp =>Pa * (P + n) = p A * (p) + J p+7r V Amo1 dp), 

respectively, we finally get: 

p A (x A < l,p + n) = p A * (p) + f* +n V A , m oi dp + RT ■ In x^ (Eq.2.136) 

Or, for the equilibrium between pure solvent and solution: 

RT ■ \nx A = — J p P+ir V Aimol dp (Eq.2.137) 

For dilute solutions we again use the Taylor series expansion to replace the logarithm simply with the 
molar fraction of the solute x B . Also, the molar volume of the liquid phase V A mo i is nearly independent 
of pressure, and this leads to a fairly simple expression: 

nV A ,mol = RT-x b (Eq.2.138) 


Considering that 

v _ n B _ n B 
n A +n B n A 

we end up with an expression similar to the ideal gas equation, namely: 
nV = n B ■RT 

called the vant Hoff equation of osmotic pressure. 


(Eq.2.139) 


(Eq.2.140) 


For illustration of the osmotic effect and its practical relevance, consider the following problem: 
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Example 2.6: 

Compare the minimum amount of energy needed to prepare 1 m 3 (= 1 ton) of drinking water 
from sea water via (i) simple boiling, and (ii) reversed osmosis (= squeezing sea water through 

a membrane to remove the salt). The salt (mainly: NaCl) content of sea water is 3.5 wt%, the 

k j 

molar evaporation heat of water is A H V mol = 40.66 • 

Solution: 

(i) For simply boiling the water and collecting the pure water vapor, you need at least the 
energy amount: 

le 6 g kj 

Q — n - A H Vmol = — — 7 ——- • 40.66 • —^ = 2.26e9 kj 
v v,moi lg g / mo i mol J 

(ii) For osmosis, the minimum amount of energy corresponds to the reversible work AG, which 
is given by the difference in chemical potentials of pure water and sea water, respectively: 

HC=n-0i i .-l. i ) = 

The molar fraction of water in sea water is given as: 

7i// 2 o 1000/18 

%A = n Hz o + 2 ■ n WaC , = 1000/18 + 2-35/(58.5) 

Inserting this in the equation for AC, we finally get for the minimum energy amount needed 
for removing the salt from 1 ton of sea water by so-called reversed osmosis: 

AC = W rev = 2.99e6 kj = 0.83 kWh 

In the internet, you find that in technical practice about 1.5-2 kWh of energy is needed to 
prepare 1 ton of drinking water by reversed osmosis, which is about twice the amount needed 
for the reversible process (= 50% efficiency!). Note that still, in terms of energy consumption, 
reversed osmosis is about 500 times more efficient (and economical) than simple boiling! 

Even if the energy to achieve the pressure needed to squeeze the sea water through the 
membrane is obtained via photovoltaics, a solar cell of efficiency of only 10 % still leads to an 
overall energy consumption efficiency osmosis : boiling = 50 : 1 !!! 


79 


Download free eBooks at bookboon.com 










Basic Physical Chemistry 


Thermodynamics 


Summary of colligative phenomena for thermodynamically ideal dilute solutions 
(only the entropy of mixing is considered (athermal mixtures!)): 

(i) vapor pressure suppression: 

For a dilute solution, the vapor pressure is lower than that of the pure solvent, as 
described by the law of Raoult: P = x a ' Pa * 

(ii) boiling temperature increase: 

A dilute solution shows a higher boiling temperature than the pure solvent. This effect 
depends only on molar solute concentration and the solvent properties expressed by the 

11- • , A m R-Ta* 2 -Ma mu 1 

ebulioscopic constant: AT -- ■ — 

A v h a m b ™a 

(iii) melting temperature decrease: 

A dilute solution shows a lower melting temperature than the pure solvent. This effect 

depends only on molar solute concentration and the solvent properties expressed by the 

. m R-T a * 2 -M a rnu 1 
cryoscopic constant: AT -- ■ — 

A s h a m b ™-a 

(iv) osmotic pressure: 

The osmotic pressure of a dilute solution is given by the van’t Hoff equation: 

nV=n B ■RT 

All these phenomena are, in principle, useful to “count the number of dissolved particles”. Note 
that osmotic pressure is the most sensitive, but restricted to comparatively large solute particles 
which are unable to pass the membrane. 


We conclude this section on phase equilibrium of two-component systems with a simple thermodynamic 
model of saturated solutions, neglecting contributions by solvation enthalpy, i.e. A H mix = 0: 


We can treat this equilibrium exactly the same way as the lowering of the freezing temperature in 
solutions, but this time the solute B is found in the pure solid phase, therefore: 


A s H b 
R 


K + ^) =ln *» 


(Eq.2.141) 


Here, the solubility parameter A S H B /(RT B *), with A S H B the melting enthalpy and T B * the melting 
temperature of the pure solute B, defines the solubility of a given component B. Note that the solubility 
of B for our ideal mixture reaches 100% at the melting temperature, i.e. we have a pure liquid of B in 
this case. At lower temperature, the solubility continuously decreases, reaching 0 at T = 0 K. Also, the 
solubility at given reduced temperature T/T B * is larger the smaller our solubility parameter. 
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2.6 The chemical equilibrium 

If we consider a spontaneous chemical reaction at constant pressure and temperature, it takes place, like 
any spontaneous process, only if the change in Gibbs free enthalpy AG < 0! Considering the Gibbs- 
Helmholtz-equation AG = AH —T AS, there are different possibilities for such so-called exergonic 
spontaneous reactions: the best case is a reaction releasing heat to the environment and creating 
disordered products, like for example burning of fuel: here, both the entropy of the reactants and of the 
environment are increasing, which is very favorable according to the 2 nd principle of thermodynamics. 

Another possibility is an exothermic process, i.e. heat release from the system to the environment, while 
forming ordered products. If the overall change in entropy is still positive, i.e. AS tot = AH/T — AS > 0, 
this process may also be spontaneous, e.g. crystallization of a solute from a saturated solution. Similarly, 
an endothermic process is possible if the increase in disorder compensated the loss in entropy of the 
environment due to the heat transfer, i.e. A S tot = AS — AH/T > 0. A simple example of this type of 
process is the dissolution of salt in water, often leading to a decrease in temperature (endothermic process). 

Finally, there is no possibility of a spontaneous endergonic process with 100% conversion (conversion = 
molar percentage of products formed from the original chemicals (educts), 100% conversion therefore 
means the reaction is complete!), i.e. endothermic and creating ordered products at the same time. 
The four different combinations of reaction enthalpy and reaction entropy are summarized in the table 
(tab.2.1.): 


change in enthalpy 

change in entropy 

spontaneous reaction ? 

exotherm (AH < 0) 

AS > 0 

yes, AG < 0 

exotherm (AH < 0) 

AS < 0 

yes, if |T • AS| < \AH\ 

endotherm (AH > 0) 

AS > 0 

yes, if T ■ AS > AH 

endotherm (AH > 0) 

AS < 0 

no, AG > 0 


Table.2.1: Conditions for spontaneous (exergonic) chemical reactions 


If the conversion is smaller than 100%, however, and a mixed phase of starting materials and products is 
formed at the chemical equilibrium, also an endergonic process may take place, at least up to a conversion 
rate of 50%. Thermodynamically, this is based on the fact that, in addition to the free enthalpies of the 
pure products and educts, we also have to consider now contributions from the free enthalpy of mixing. 
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Figure 2.32. illustrates the total free enthalpy of our reaction system versus conversion £ in case of an 
endergonic process (free enthalpy of the products is larger than that of the educts), and for a process 
with AG = 0 . Note that the stable state, i.e. the minimum in free enthalpy, is reached at a conversion 
of 0.5 (or 50%) for the latter case, since here the minimum in free enthalpy of mixing (or maximum 
in entropy of mixing) is the driving force for the process. On the other hand, an endergonic process 
still may proceed until the mixing free enthalpy exactly compensates the corresponding reaction free 
enthalpy, that is, to conversions between 0 and 0.5. 




Fig. 2 .32: Reaction free enthalpy versus conversion and the chemical equilibrium 

Mathematically, we can treat the chemical equilibrium exactly in the same way as the phase equilibria. 
Let us consider the following reversible chemical reaction: 

A + B ^ C + D (Eq.2.142) 

This reaction reaches its equilibrium conversion if the chemical potentials of educts and products are 
identical, or: 

Ha + Hb = He + Hd (Eq.2.143) 

If the equilibrium state is a homogeneous reaction mixture of all components, we get: 

Ha * +RT ■ In X 4 + ii B * +RT -lnx B = He * + RT ■ 1 nx c + Hd * + RT ' l n x d (Eq.2.144) 


Or 


—RT - In x c — RT • \nx D + RT • \nx A + RT ■ \nx B = He * * ~Ha * ~Hb * 

-RT -ln x f^=-RT- ln K x = A R G * (Eq.2.145) 

X A X B ± 
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With K x the constant of chemical equilibrium, which, according to the last equation, depends on reaction 
temperature and the difference in free enthalpy of pure products and pure educts A R G *. For neither 
exergonic nor endergonic reactions, i.e. A R G *= 0, we obtain K x = 1, which, as we have already 
discussed, corresponds to a conversion <f = 0.5. For exergonic reactions A R G *< 0, K x > 1, and the 
conversion at equilibrium is larger than 50%. 


Let us conclude this chapter about the chemical equilibrium with a brief visit of the famous principle of 
Le-Chatelier: the equilibrium of a chemical process reacts on the change of an external parameter like 
temperature or pressure such as to minimize the change. Simply spoken, if we consider an exothermic 
chemical reaction and increase the reaction temperature, the endothermic back reaction is favored to 
compensate, leading to a shift of the chemical equilibrium towards smaller conversion. Formally, this is 
fairly simple to understand if we discuss the above equation in a little more detail: 


—RT ■ In K x — * —fA^S * 


In K x 


&R H * , &rS* 
RT R 


d\nK x _ A r H* 
dT ~ RT 2 


(Eq.2.146) 

(Eq.2.147) 

(Eq.2.148) 


If the process is exothermic, i.e. A R H *< 0, K x therefore has to decrease with increasing temperature. 
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Another example of the principle of Le-Chatelier is the pressure dependence of the Haber-Bosch-reaction, 
the formation of NH 3 from the elements H 2 and N 2 (see Eq.(2.149)). Since formally this reaction leads 
from five gas molecules to two, an increase in pressure favors the process. 

3 H 2 + 2 N 2 ^ 2 NH 3 (Eq.2.149) 


Example 2.7: 

Consider a simple chemical equilibrium of type A ^ B. Calculate the relative change in 
conversion of A to B, if the reaction temperature is increased from T = 25°C to T = 100°C, and 
the reaction enthalpy A R H = —80 kJ/Mol. 


Solution: the variation of the chemical equilibrium constant, which provides a measure for the 
conversion, with temperature is given as 


d\nK x _ ArH_ 
dT ~ RT 2 


or 


\nK x (T 2 )-\nK x (T 1 ) = \n^- ) 


*rh r i n 

R V7i T 2 J 


k x (t 2 ) 

K x (Ti) 


n 

1\ 

r-8oooo 

/ 1 

_ 1 \1 

V7i 

t ) 

= ^ [ 8.314 ' 

V298 

373/ J 


This means, for example, that if at room temperature we have a conversion A to B of 50 % 
( K x = 1), at T = 100 °C the conversion is only about 0.15 % ( K x = x A /x B ~ 0.00151/1). Our 
result is also a confirmation of the principle of Le-Chatelier: if you increase the reaction 
temperature for our strongly exothermic chemical reaction, the equilibrium is shifted in favor of 
the endothermic backreaction. 

As an additional exercise (not typical for Physical-Chemical problems, but related to general 
chemical calculus), let us determine the exact conversion corresponding to K x = 0.00151: 

K x = — = — = 0.00151 => x = °' 00151 = 0.001508 

* x B 1-x 1 + 0.00151 


The result proves our approximation x B « 1 to be acceptable in this case. 
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2.7 Reaction energy 

As we have learned in the last chapter, the temperature dependence of the chemical equilibrium depends 
on the reaction enthalpy. To measure this reaction enthalpy, one can simply determine the change in 
temperature if the reaction is carried out at isobar conditions and 100% conversion, i.e. A R H *= c p - AT. 
However, if it is not possible to monitor the reaction itself directly, one can apply the principle of Hess, 
which is based on the fact that the reaction enthalpy or the reaction energy both are quantities of state, 
that are independent of the process itself. Consider, for illustration, the reaction enthalpy of the hydration 
of ethylene to ethane. Instead of observing this process, one could measure the heat of burning of ethylene 
and hydrogen as the educts, or pure ethane as the product, respectively, with pure oxygen (see fig. 2.33): 




Fig. 2 .33: The principle of Hess - conservation of energy/enthalpy for chemical reactions irrespective of the reaction pathway 


The hydration enthalpy is then simply the difference between the two heats of burning AH = Q 1 — 
Q 2 . Note that in this way any reaction energy or enthalpy can be determined, taking into account 
also that the enthalpy of a pure element in its stable modification at standard conditions (T = 298 K, 
p = 1 bar) is defined as zero. For illustration, consider the burning of hydrogen with oxygen to water at 
standard conditions H 2 + V 2 ^2 H 2 0. The heat of burning in this case is directly the reaction enthalpy 

A r H*=H H2 o-H H2 - 1 / 2 H 02 . 

According to the convention, we can also write A R H *= H Hz0 . In this way, one can determine, based on 
some simple model reactions, a whole set of reaction enthalpies mainly on calculations. More importantly, 
one can also determine the enthalpy of individual chemical groups like CH 2 or COOH, leading to the 
total enthalpy of more complicated organic molecules without the necessity of more experiments. 
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To complete this puzzle, we need to know how to transfer our reaction enthalpy from standard temperature 
to any reaction temperature. Again, we utilize the fact that the enthalpy is a quantity of state irrespective 
of the process itself. Consider two possibilities for the reaction from A(T) to B(T+AT): either you can 
first heat the educt A and then carry out the reaction at temperature T+AT, or you can carry out the 
reaction first at T and then change the temperature of the product to T+AT. Either way, the overall 
change in enthalpy should be the same. This leads us to the important Kirchhoff law: 


£c p (A) ■ 

AT + A r H t+at = A r H t + £ c p (E) ■ AT 

(Eq.2.150) 

ArH t+at 

— A r H t — 

(£ c p (E) — £ c p (A)) • AT 

(Eq.2.151) 

d(A R H) 

= (£ c p (E) 

— £ c p (A)) ■ dT 

(Eq.2.152) 

^(A R H)^ 

= £ c p (E) - 

- £ c p (A) 

(Eq.2.153) 

^(ArH)) 

£ ^iCpi 


(Eq.2.154) 


That is, the change in reaction enthalpy with temperature is simply given by the difference in isobar heat 
capacities of products and educts. 


7 + i 


T 


1 



II 


/ 





J 


I 



2 



Fig. 2.34: Kirchhoff law, two different routes from educt A(T) to product B(T+AT) with identical overall 
reaction enthalpy 
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3 Kinetics 


The velocity of a chemical reaction is defined as the change in concentration (for gas reactions: partial 
pressure) with time, taking into account the stoichiometric coefficients of the respective components. 
In experimental practice, this concentration can be quantified in various ways, for example: photometry 
(measurement of absorbed light at specific wave length following the law of Lambert-Beer), conductometry 
(measurement of electric conductance, in case either ions are formed or consumed during the chemical 
reaction), or polarimetry (measurement of the rotational angle of polarized light in case optically active 
components are involved). 

In the simple case of elementary reactions, that is, reactions with one single reaction step from educt 
to product (involving an excited state, see also figure 3.3 at the end of this chapter!), the decrease in 
concentration of component A with time for the reaction 

n^A + n B B —> tiqC (Eq.3.1) 

is given as: 


- d J± = k .C A nA. CB 


kb 


(Eq.3.2) 


More general, the reaction velocity for this reaction is given as: 

_ 1 dc a 1 dcR 1 dec 

ua dt n# dt nc dt 


(Eq.3.3) 


The velocity constant k is dependent on temperature, and may be increased by adding a catalyst, as we 
will see when we introduce the Arrhenius equation below. The total order of the reaction is given by the 
sum of all stoichiometric coefficients of the educts, i.e. n = n A + n B . Here, it should be stressed again 
that this reaction order is only identical with the molecularity of the reaction, i.e. the total number of 
molecules reacting, if the whole process consists of a single reaction step! 
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Importantly, for more complex reactions which usually consist of multiple reaction steps, one has to 
distinguish clearly between the order of the reaction and the molecularity of the individual reaction 
steps: the order corresponds to a macroscopic experimental quantity, namely the exponent with which 
the reactions velocity scales with the respective concentration (see Eq.3.3). On the other hand, the 
molecularity of a single reaction step corresponds to the microscopic mechanism of this reaction step. As 
a consequence, a chemical reaction which is experimentally found to be of 2 nd order must not necessarily 
be a bimolecular reaction, since a combination of several reaction steps of different molecularities can lead 
to an overall reaction order of 2. Actually, a variety of possible combinations of elementary reaction steps 
(= proposed mechanisms of a chemical reaction) may lead to the same experimental result concerning 
the order of a chemical reaction! 

In the first section of this chapter about chemical kinetics, we will mathematically treat the kinetics 
of the simplest case possible, elementary reactions with only a single educt species (section 3.1.1), or 
bimolecular reactions with two different educt species A, B (section 3.1.2.). 
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3.1 Elementary reactions (one single reaction step, molecularity = order) 

3.1.1 Elementary reactions of type A -> B 

In case of a single educt A, the differential equations can easily be resolved. Here, we distinguish reactions 
of orders 0 to 3: 


Order 0: Note that this is no elementary reaction in the strict sense, since the actual reaction mechanism 
is a two-step process, 1 st step the combination of educt and a catalyst, and 2 nd step the formation of the 
product which is released from the catalyst: 


A + Cat. -> A ••• Cat. -> B (Eq.3.4) 

In case of high concentrations of the educt, the catalyst is completely covered, and the reaction velocity 
therefore is constant: 


~^r = k 

dt 

(Eq.3.4) 

d C 4 — kdt 

(Eq.3.5) 

f CA(t) dc A = ff -kdt 

J c, 4,0 A J 0 

(Eq.3.6) 

£ 

r~t- 

II 

£ 

O 

1 

(Eq.3.7) 

i.e., the concentration decreases linearly with reaction time. 


Order 1 : A typical example of unimolecular elementary reactions are radioactive decay processes. The 
reaction velocity here is proportional to the actual educt concentration, and therefore decreases with 

time as: 

1 

gr S' 

II 

(Eq.3.8) 

^ = -kdt 

ca 

(Eq.3.9) 

fC A (t)dc A=j t 

J C A ,0 C A J 0 

(Eq.3.10) 

In c A (t ) — In c A 0 = —k ■ t 

(Eq.3.11) 

C; l(t) = C A q ■ exp(-fc ■ t) 

(Eq.3.12) 


The educt concentration decreases exponentially with time. 
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Order 2: Examples are any bimolecular homologous reactions, for example NO + NO -> N 2 0 2 . The 
reaction velocity in this case scales with the educt concentration squared: 


_d£A_ _ k 2 
dt ~ K ° A 


(Eq.3.13) 


dc A 


CA 


2 


= —kdt 


(Eq.3.14) 


rC A (t) dc A 
C A, 0 C A 2 



(Eq.3.15) 


1 + — =-k-t 


C A (t) C A n 


The educt concentration varies hyperbolically with time. 


(Eq.3.16) 


Order 3: Because they afford the simultaneous collision and reaction of three molecules, elementary 
reactions of this type are highly improbable. Nevertheless, we will treat this elementary reaction as an 
additional exercise how to solve simple differential equations. The reaction velocity scales with the educt 
concentration cubed: 


d c A _ Is • r 3 

~ at ~ K Ca 


= -kdt 

CA 3 


jCAVdc± = ft kdt 

J C A) 0 c A 3 J 0 


■ + -- - — —k ■ t 


2-c A (t) 2 2-c a> o 2 


(Eq.3.17) 

(Eq.3.18) 

(Eq.3.19) 

(Eq.3.20) 


To compare the conversion curves for these different elementary reactions of reaction orders 0, 1, 2 and 
3, we consider identical starting concentrations c A 0 = 1 mol/L and identical magnitudes of the velocity 
constant k = 1 (in respective units). For 0 order, this means k — 1 mol /(L ■ 5 ), for 1 st order k = 1/s, 
for 2 nd order k — 1 L/(mol ■ 5 ), and for 3 rd order k — 1 L 2 / (mol 2 ■ 5 ). In this case, the initial reaction 
velocities should all be equal to — ^ = 1 mol/(L ■ s), and the velocity (or decaying slope) decreases 
with time for reactions of order >0, with the curvature more pronounced the higher the reaction order. 


Note that, for reactions of 1 st order, we can define a half-conversion-time, i.e. the reaction time it takes 
until 50% of the educt amount is consumed: 


t 


_ In 2 
V2 - — 


(Eq.3.21) 
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All 4 conversion curves are compared in figure 3.1, including an indication of t m for the reaction of 
1 st order: 



t 


Figure 3.1 : Time-evolution of the normalized educt concentration for elementary reactions of orders 0,1,2 and 
3 (from left (= linear decay) to right), including t m for the reaction of 1 st order as arrows. 
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Note here that t 1/2 does not make sense for elementary reactions of orders other than 1, which can easily 
be seen from a simple numerical example: for instance, if a reaction of order 0 takes 30 minutes for the 
conversion of 50% of the original amount of the educt, it will take only 15 minutes to reach 25%. On 
the other hand, if a reaction of 2 nd order takes 30 minutes for the conversion of the first 50% of educt, 
25% will be reached after another 60 minutes. 


Example 3.1: 

Consider a radioactive material with halftime 30 years. After how much time has the radiation 
level decayed to below 1% of its initial value? 

Solution: We use the following equations for an elementary reaction of 1 st order: 
c A (t) = c A (t = 0) ■ exp [—k ■ t], and k = 

1 1/2 


Therefore, k = — = — = 0.0231 y~ x . With 

*i/2 30 y y 


CA ^. = — = exp[—0.0231 y 1 ■ t] we finally get 

c A (t=0) 100 J J J b 


In 0.01 
-0.0231 y 1 


199.4 y. 


It will take about 200 years or several generations until the radiation level is below 1%! 


3.1.2 Elementary bimolecular (and 2 nd order!) reactions of type A + B -> C 

Many elementary reactions of 2 nd order comprise two different types of educts A, B. In this case, the 
general solution of the differential equation is not any longer trivial, as in case of single component 
elementary reactions. For an elementary 2 nd order reaction of type A + B —> C the velocity is given as: 


d c A 
dt 


dcr j 

— = k ■ c A ■ c B 


(Eq.3.22) 


To solve this equation, which contains three variables at this stage, we express the actual educt 
concentrations via their respective starting concentrations minus the concentration of product already 
formed, i.e.: 


- k ■ { C A,0 - c c) ■ ( C B ,0 “ C C ) 


(Eq.3.23) 


or 


_ dec _ 

0 c A,0-Cc) m ( c B,0~CC ) 


k ■ dt 


(Eq.3.24) 
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This seemingly complicated differential equation can easily be solved via so-called partial fraction 
integration. First, we split the left side of the equation as: 


dc c 


K r dc c + _K 1 dc £ _ =k . dt 


0 c A, 0 -Cc)'(cb, 0 -Cc ) ( CA,0~ec ) ( C B,0-Cc) 


(Eq.3.25) 


If we compare the two split fractions with the original fraction, we obtain ■ (c 5 0 — c c ) + K 2 ' 
(c A)0 ~ c c ) = 1. The constants and K 2 can be determined via the boundary conditions, i.e. atc c = 0 
we obtain ■ c B0 = 1 — K 2 ■ c A>0 . At c c = c A 0 , we get 


K, = 


1 ( c B,Q- c A,0 ) 


(Eq.3.26) 


and therefore 


K ? = 


(ca,o~ c b,o) 


(Eq.3.27) 


The general quite complicated solution of this equation is then given as: 


, ln ^.o-(c B ,o- CC ) = k . t 


( c B,0~CA,o) c B,0 m ( c A,0-Cc) 


(Eq.3.28) 


There exist two procedures to obtain a mathematically more simple solution of the differential equation, 
using special starting conditions: 

i. If the starting concentrations of both educts are identical, i.e. c A 0 = c B>0 , the differential 
equation (Eq.3.24) is simply given as: 


dec 


— k ' dt 


0 C A,0-Cc ) 

In this case, we can easily integrate the equation with appropriate boundary conditions 


(Eq.3.29) 


fo C 77 l£ ^ = fo k - dt 


(ca.o-cc) 0 


(Eq.3.30) 


And obtain: 


( c A,0~ c c) C A,0 


= k ■ t 


(Eq.3.31) 
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ii. If we choose one component in large access, e.g. c A 0 » c B 0 , the concentration of this 
component will nearly remain constant during the whole reaction, and again the differential 
equation (Eq.3.24) becomes fairly simple: 


dcr i i. 

r \ — fC ' Ca n " 

Ob,o-cc) A 

(Eq.3.32) 


(Eq.3.33) 

c c = c B ,o - c B ,o ■ exp [-k ■ c Afi ■ t] 

(Eq.3.34) 


This 2 nd approach is valid for any solvolysis reaction, i.e. a chemical process where component B reacts 
with solvent A in dilute solution, for instance methanolysis of carbonic acid chlorides. 
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3.2 The kinetics of more complex multistep chemical reactions 

3.2.1 The Lindemann-mechanism as an example of a seemingly unimolecular decay process 

In practice, for a seemingly unimolecular decay process of type A —> B in the gas phase, sometimes 
you encounter reaction orders between 1 and 2, depending on conversion. The simplest mechanism 
to explain this strange phenomenon was formulated by Lindemann: a combination of a bimolecular 
equilibrium reaction and a monomolecular final reaction step of an activated species. The Lindemann 
multistep mechanism reads as following: 

A + A A + A *-> B (Eq.3.35) 


A* is the activated species formed by a bimolecular collision of two molecules of the educt A. This species 
can either deactivate again by bimolecular collision, or form the product species B by a monomolecular 
decay. To solve this kinetic problem, we can formulate the following differential equations: 

dp B i 

-^T = k i ' Va* (Eq.3.36) 


for the monomolecular product formation from the activated species, and 
= k 2 ■ va ~ k - 2 ■ Va ' Va* ~ ' Va * 


(Eq.3.37) 


for the change in partial pressure of the activated species with time. Note that, in contrast to the previous 
kinetic equations, for the gas phase reaction considered here we have replaced concentrations c with 
partial pressures p. To solve this equation system even without integration, we may treat the partial 
pressure Va * as constant, since the activated species is formed and consumed at equal rates of a certain 
time regime of our complex reaction (= stationary state). In this case, 

^f 1 = k 2 ■ va ~ k -2 ■ Va ' Va* -h-p A * = 0 (Eq.3.38) 


or 


Va* 


^2-va 2 

k-2'VA+ki 


and consequently 

dp_B_ _ kykrVA 2 
dt k-2'PA+ki 


(Eq.3.39) 


(Eq.3.40) 


Consider the denominator of the right side of this rate equation for the formation of product B in more 
detail: at high pressure of the educt Va or at an early stage of the reaction, k 2 ' Va » k 1 • and therefore 


dps_ _ kykrPA 
dt k-2 


(Eq.3.41) 


95 


Download free eBooks at bookboon.com 







Basic Physical Chemistry 


Kinetics 


In other words, we find a reaction of 1 st order at low conversion. This is plausible since at high Pa all 
bimolecular reaction steps due to the high collision probability are rather fast, and therefore the slowest 
reaction step, the monomolecular decay of activated species A*, determines the overall reaction order. 

At low pressure of the educt Paox at a later stage of the reaction, k_ 2 ■ p A « /q, and therefore 

dPB i 2 

^T = h-PA Z (Eq.3.42) 

In this case, we find a 2 nd order reaction, since now the collisions become rarer, and therefore the 
bimolecular collisions as the slowest reactions steps define the overall reaction order. 

The Lindemann mechanism is the simplest form of a combined multistep reaction to explain the change 
in reaction order from 1 to 2 upon proceeding conversion. Note here that the elementary reaction of 
order 0 described above also is a two-step reaction, explaining in this case the change from reaction order 
0 at high educt concentration or low conversion, to reaction order 1 at lower concentration or higher 
conversion, in which case not any longer all reaction sites of the catalyst will be occupied! 

3.2.2 The Michaelis-Menten mechanism as an example of a multistep enzymatic reaction 

A very famous multistep reaction similar to the Lindemann mechanism plays an important role in 
biological enzyme reactions, the Michaelis-Menten process. In this case, an enzyme is serving as a catalyst, 
and the starting reaction is a bimolecular collision leading to a combined enzyme-substrate-complex. This 
activated complex can either form the product or be cleaved without changes to the substrate, releasing 
in both cases the free enzyme E which then may react again. In analogy to the Lindemann mechanism, 
we may formulate this mechanism as 

E + S^ES^E + P (Eq.3.43) 

Let us name the reaction velocity constants related to the equilibrium reaction E + S ^ ES/q and /c_ 1? 
and the constant related to the formation of the product P k 2 . The maximum reaction velocity or highest 
turnover rate is reached if all enzyme exists in form of the enzyme-substrate complex ES. Next, we define 
the so-called Michaelis-constant, which corresponds to the substrate concentration at half-saturation 
of the enzyme, i.e. c E = c ES = c Eq ■ 0.5, with c Eq the total enzyme concentration. This constant, in 
case k 2 « k_ 1 (= Michaelis-Menten approximation), can be derived as following. Analogous to the 
Lindemann mechanism, we consider the rate equation for the stationary state: 

= 0 = k 1 -c E -c s ~ (fc_! + k 2 ) ■ c ES ~ k x ■ c E ■ c s - k_ 1 ■ c ES = 0 (Eq.3.44) 


or 


k-i . c es 


ki c E 


(Eq.3.45) 
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Therefore, the Michaelis-constant is simply given as 


Km C S • (v max /2) - —- 

*1 

The product formation reaction velocity in general is given as: 
^- = k 2 -c ES = k 2 -c s -c E - K m ~ x 


(Eq.3.46) 


(Eq.3.47) 


and 


kz ' c es — k 2 ■ ( c Eq — c E ) — v max — k 2 ■ c E 


(Eq.3.48) 


Therefore, 


k 2 ' c s ■ c E ■ K m 1 — v max — k 2 ■ c E 


(Eq.3.49) 


or 


_ _ v max _ 

E {k2'Cs'K m ~ 1 +k 2 ) 


(Eq.3.50) 
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Inserting this expression in the reaction velocity finally yields the following Michaelis-Menten relation 
describing the product formation velocity in dependence of substrate concentration, Michaelis-constant 
and maximum reaction velocity only: 


~ ^2 ' C S ' C E ' Km 


v max '^2 ' c S'Km _ v max ' C S 

( k 2 -Cs-K m 1 +k 2 ) c s +K m 


(Eq.3.51) 


This reaction velocity is plotted in the following figure 3.2, illustrating the practical meaning of important 
parameters such as maximum velocity or Michaelis-constant for biological enzyme-substrate-reactions. 



3.3 Activation energy 

So far, we have not considered the dependence of reaction velocity on temperature. Obviously, an increase 
in temperature leads to an increase in velocity, a concept which was first addressed by Arrhenius (on 
a purely phenomenological basis!): he assumed that all reactions afford a transitional state located, in 
respect to the educts, higher in energy by the amount of the activation energy. The ratio of activated to 
non-activated species at given temperature now is given by the Boltzmann factor (see chapter 2, Eq.2.15), 
and according to Arrhenius the temperature dependence of the reaction velocity constant is given as: 

fe = i4-exp[-^] (Eq.3.52) 
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RT is the thermal energy per mole, and E A the activation energy per mole of reacting educt particles. 
Note that in this simple purely energetic approach A is merely a phenomenological fitting parameter 
corresponding to the maximum possible reaction velocity, i.e. k- A at T -> oo. A more sophisticated 
theory following Eyring, and beyond the scope of this textbook, shows that A is of entropic origin, a 
concept however not yet known to Arrhenius. 


energy 



reaction coordinate 

Figure 3.3: Energy profile of a chemical reaction, and activation energy (Arrhenius concept) 


The origin of the activation energy barrier E A is easy to understand if we consider that, for the formation 
of a new reaction product, first chemical bonds of the original educt molecules have to be broken. This 
bond energy is in the order of 40-60 kj/mole, explaining also the famous rule of thumbs that if one 
increases the temperature of a chemical reaction from room temperature by 10°C, the reaction velocity 
typically doubles. 
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Example 3.2: 

Consider a chemical reaction with activation energy E A = 50 kj/mole. Which temperature do 
you need to enhance the reaction velocity, compared to that at room temperature, by a factor of 
(i) 10 or (ii) 100, respectively? 


Solution: we use the Arrhenius equation k = A ■ exp[— E A /(RT)]. Note that we will not need A 
to solve this equation, since we consider only a relative change in velocities in this problem, 
and not absolute values of the velocity constant: 


k(T=293 K) _ J_ _ A-exp[-E A /(R-293 K)] _ exp[-E A /(R-293 K)] 
k(T ) ~~ 10 ” A-exy[-E A /{RT )] “ exp[ -E A /(RT)] 


exp [- —■ (— — -)] 

^ I R \293K TJ\ 


In 0.1 


50000 _ / 1 _ 1\ 

8.314 \293 K t) 


- = In 0.1 • ——— H—-— = 0.003030 => T = 330 K = 56.9 °C 

T 50000 293 K 


Using the same approach we get for an increase in velocity by a factor of 100: 

- = In 0.01 ■ —— H—-— = 0.002647 => T = 377.8 K = 104.7 °C 

T 50000 293 K 
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4 Electrochemistry 

Electrochemistry plays a very important role both in technology (for example chemical sources of 
electric energy) and chemical analytics (for example potentiometric or conductivity measurements, both 
providing a quantitative measure for the concentration of charged solute particles). This chapter is divided 
in two parts: (1) the transport of charged particles in an electrolyte solution (= ion conductivity), and (2) 
the electrochemical equilibrium. In the 2 nd part, we will address the different types of electrochemical 
electrodes, and finally combine a pair of electrodes to obtain an electrochemical Galvanic cell. 

4.1 Electric Conductivity 

An electrochemical cell used, for example, for electrolysis, consists of the source of an electric voltage, 
as well as of two chemical electrodes immersed typically in an aqueous solution of an electrolyte. The 
electric circuit consists of negative charges, the electrons, which migrate through all metal parts of this 
cell, and the ions both positive (cations) and negative (anions) migrating through the aqueous solution 
towards the respective electrodes. At the electrodes, the ions either release an electron or accept an 
electron, leading to a respective chemical reaction (oxidation if an electron is released from the ion to 
the electrode, or reduction if the ion accepts an electron from the electrode). 
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For example, let us consider two copper electrodes immersed in a dilute aqueous hydrochloric acid 
solution (this cell is formally expressed as Cu/HCl/Cu, / signifying phase boundaries between the solid 
copper electrodes and the liquid electrolyte solution, respectively), connected to an external electrical 
power source to enforce an electrolytic process. The chemical reactions taking place at the two electrodes 
are then given as: 

cathode: 2 H + + 2 e" —► H 2 1 

anode: 2 Cl" —> Cl 2 T + 2 e" 

cell reaction 2 H + + 2 Cl" —» H 2 f + Cl 2 1 

Note that the number of electrons exchanged per formula equals z = 2 for this example. 

Michael Faraday discovered the following two laws quantitatively describing the electrolysis: 


1. m H2 ,m C i 2 ~I ■ t 


(Eq.4.1) 


The mass of the material formed at each electrode is proportional to electric current multiplied 
by time, or total electric charge. 


2 mH 2 _ M H2 

™ci 2 m ci 2 


(Eq.4.2) 


The mass ratio of materials formed at the respective electrodes is identical to the ratio of the 
molar masses. 


Introducing the Faraday constant, we can combine these two equations to express the total charge 
exchanged during the electrolysis: 

Q = I -t = z-N A -e-^ = z-F (Eq.4.3) 

meaning that the molar amount of material formed at the respective electrodes ( m /M) corresponds to 
the transport of z mole electrons (or the charge zF Coulomb, with the Faraday constant F = 96484.6 C 
mol x ) from electrode to electrode via the metal wiring of the electrolysis cell.. 

To discuss the conductivity and ion migration in a quantitative way, let us consider a more schematic 
representation of the electrochemical cell: 
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Two electrodes at distance 1 (within the electrolyte solution) both of area A are connected to an electric 
power source with voltage U, and immersed in an electrolyte solution. In case the area of the electrodes 
is larger than the electrode distrance squared ( A»l 2 ), a homogeneous electric field E of magnitude 
U/l is created. Due to electrostatic interactions, the ions of charge z t ■ e are accelerated towards their 
respective electrodes by an attractive force F E = z t ■ e ■ E . This force is balanced by the frictional force 
F r = 6n ■ r] ■ Ri ■ v t ,, with rj the viscosity of the solvent (for water at room temperature, r\ - 1.00 mPa s), 
R the radius of the ion, v t and the velocity of the migrating ion (index i referring to ions of species i, 
for example H + oder Cl ) . At force balance, the acceleration is zero and the ions are migrating with the 
constant velocity: 


Vi = 


Zj-e-E 

6n-7]-Ri 


(Eq.4.4) 


The ion velocity therefore also depends on the experimental setup, i.e. the electric field strength. One 
therefore defines, as a more general quantity just depending on material properties, the so-called 
ion mobility: 


u t 


Vi _ z E e 
E 6n-rj-Ri 


(Eq.4.5) 
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This microscopic property is related to the electric conductance (or resistance) of the electrochemical 
cell, which can more easily be measured than the migration of individual ions. To derive a relation 
between microscopic ion mobility and macroscopic ion conductance, we consider the electric current 
within the electrolyte solution, consisting of positive and negative ions migrating in opposite directions: 


F dQ + , dQ~ dN + -e , dN~-e 

I = =-I- 

dt dt dt dt 


(Eq.4.6) 


The amount of negative or positive charges (dN + ■ e, dN~ ■ e) migrating within a time step over a 
certain distance dx towards an electrode of surface area A is given by the average ion concentration and 
the migration velocity as: 


dN + -e 

dt 


d J^. A -(^) = F ■ c + ■ A • (—) 

A-dx \dt J- 1_ \dt J .|_ 


= F ■ c + ■ A ■ v + 


(Eq.4.7) 


with c + = n + • z + - c. c is the molar concentration of the electrolyte, n + the number of cations per 
formula, and z + the elementary charge number of these cations. For example, for the salt MgCl 2 , n + = 1 
and z + = 2. 


In total, we obtain for the electric current: 


/ = F ■ A ■ (c + ■ v + + c • v ) = F • A • (c + ■ u + + c * u ) ■ y 


(Eq.4.8) 
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Comparing this equation with Ohms law 


i = Z =u .±.i =u .±. K 

R Ip l 


(Eq.4.9) 


we see directly that the resistance R depends on the geometry of our electrochemical cell (= cell 
constant A/l) and, via the specific resistance p or its reciprocal, the specific conductivity k, on the ion 
concentration and material properties of the ions (charge, size). Comparing Eqs. (4.8) and (4.9), we 
find that the specific conductivity, as a macroscopic material property, depends on the microscopic ion 
mobilities and the electrolyte concentration as: 


k = F ■ (c + ■ u + + c ■ u ) 


(Eq.4.10) 


To eliminate the concentration dependence, we define the molar conductivity of the electrolyte as: 

A = — F ■ (n + • z + • u + + n~ • z~ • u ~) (Eq.4.11) 

Note that, by measuring the electric resistance of an electrolyte solution, we determine the sum of the 
ion mobilities, but not the mobility of one ion species. Importantly, we have to use alternating voltage 
for this experiment. Otherwise, electrolysis will change the ion concentration with time. In this case, we 
also would need a certain voltage before the chemical reaction can take place. One experimental setup 
to measure the electric resistance of an electrolyte solution very accurately is the compensation setup 
or Wheatstone-bridge (found in common textbooks on experimental physics). 

So far, we have shown that the measurement of the electric resistance and calculation of the molar 
conductivity only yields the sum of cation and anion conductivities 

A= F ■ (n + ■ z + ■ u + + n~ ■ z“ ■ u ~) = A + + A_ (Eq.4.12) 

To directly determine the ion mobility of one ion species, two different experimental approaches may 
be used: 

i. If the ions absorb visible light, one can directly observe the migration of a colored front 
moving in an electric field (a special electrochemical cell (formed like a “U”) is carefully 
filled with a solution containing the colored ions, and a colorless electrolyte solution on top 
in contact with the two electrodes, see textbooks for more details). A famous example is the 
measurement of the migration of the colorful Mn0 4 -ion. 


105 


Download free eBooks at bookboon.com 


Basic Physical Chemistry 


Electrochemistry 


ii. If the ions are colorless, the method developed by Hittorf may be used (see fig. 4.3.). A 
special electrolysis chamber consisting of three different compartments which are separable 
is used, and the balance of the respective change in electrolyte amount within the two 
electrode chambers after a defined amount of electric charge has been applied to the cell is 
determined, respectively, for example by volumetric titration. 



Fig. 4.3: Hittorf electrolysis chamber system (schematic, dark grey arrows = migration of cations, 
light grey arrows = migration of anions), K = cathode chamber, A = anode chamber 


The relative part of the electric current due to cation and anion migration, respectively, is given as 

i+ Q+ xi-\- A-l 

7 - Q++Q- ~ u + +u_ - "a” - + (Eq.4.13) 

/_ Q- U- A_ 

/ — 0++0_ — u++u- ~ A ~ (Eq.4.14) 


with £ + , t_ the Hittorf numbers of the given electrolyte. Note that £+ + £_ = 1. Let us consider the 
change in electrolyte amount in each electrode chamber if we apply one mole charge to an aqueous 
solution of HC1 for illustration: 


cathode chamber 

anode chamber 


- 1 mole H + 

- 1 mole CL 

change due to electrolysis 

- £_ mole CL 

+ t_ mole CL 

change due to ion migration of Cl' 

+ £+ mole H + 

- £ + mole H + 

change due to ion migration of H + 

-1_ mole HC1 

+1_ mole HC1 

overall change in HCl-amount 
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Determination of the overall change in HC1 therefore directly yields the Hittorf numbers. If you 
independently measure the overall conductivity of the electrolyte, you then can calculate the molar 
conductivities (and ion mobilities) of the individual ions. 

One should note that once you know the total molar conductivity of the electrolyte and the molar 
conductivity of one ion species, you always can calculate the molar conductivity of the 2 nd species. This 
concept can be used to determine the molar conductivity of the colorless Nations, which are also not 
suitable for direct investigation via the Hittorf method by electrolysis in aqueous solution, as following: 

1 st step: Determine the molar conductivity A of the salt NaMn0 4 by measurement of the electric 

resistance of an aqueous solution. 

2 nd step: Determine the ion mobility of the colored Mn0 4 - ions by direct observation of the ion 

migration within the U-cell, and calculate the molar ion conductivity as A_ = F ■ u_. 
3 rd step: You now may calculate the molar ion conductivity of Na + as A + = A — A_. 

Finally, let us discuss the effect of the hydration shell of associated water molecules formed in aqueous 
solution on the cation mobility. In dilute aqueous solution, the ion mobility decreases in the order 
u(H + ) > u(K + ) > u(iVa + ) > u(Li + ). This is surprising, since one would expect the Li + - ions to 
be smaller than the K + - ions, for example. However, in aqueous solution the ions are surrounded by a 
hydration shell of water molecules, and the number of these water molecules is the larger the higher 
the surface charge density: approximately, the hydration shell ofL/ + - consists of 12, that of Na + of 8, 
and that of K + of 4 water molecules, rendering the K + -ion effectively smaller and therefore more mobile 
in aqueous solution. 



Figure 4.4: The GrotthuG mechanism, explaining the high ion mobility of H + and OH' 
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An exception is H + where a special conductivity mechanism takes place (see fig. 4.4): the ion is not itself 
migrating, but electron pairs are shifted between the ion and neighboring water molecules, leading to 
an effective charge mobility about 4 times larger than that of migrating ions. A similar mechanism can 
be formulated to explain the high mobility of OH -ions in aqueous solution. 

Finally, the ion mobility also depends on temperature, since the viscosity of the solvent and therefore 
the frictional resistance is decreasing with temperature. Accordingly, ion mobility times solvent viscosity 
should be independent of temperature (“Waldens Rule”). 

4.1.1 Electric conductivity of weak electrolytes 

For so-called weak electrolytes the molar conductivity is strongly decreasing with increasing concentration 
of the substrate, for instance acetic acid dissolved in water. The reason is that the weak electrolyte is not 
fully dissociated except at infinite dilution, and the degree of dissociation is strongly decreasing with 
increasing concentration. Let us consider the example acetic acid in aqueous solution in more detail. 
The dissociation of the weak acetic acid is determined by the following chemical equilibrium: 

HAc ^ H + + Ac~ (Eq.4.15) 
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The according dissociation constant and degree of dissociation are defined as: 

[H+HAc-] 

Hdiss,c [HAc] 

_ [Ac~\ _ [Ac~] _ [H + ] _ c 0 -[HAc] 

[Ac~]+[HAc] c 0 c 0 c 0 

Inserting the degree of dissociation a in to K dissc , one obtains: 

_ a-cp-a-co _ a 2 -cp 
diss ,c (i -a)-cp (1 -a) 


(Eq.4.16) 

(Eq.4.17) 


(Eq.4.18) 


This quadratic equation may be used to calculate the degree of dissociation as a function of electrolyte 
concentration c 0 . According to Ostwald, the concentration dependence of the molar conductivity of the 
weak electrolyte, which is determined by the dissociated charges only, then is given as: 


A(c 0 ) = a ■ A ( 


(Eq.4.19) 


with Aoo the limiting molar conductivity measured at infinite dilution or 100% dissociation a — 1. 


Experimentally, conductivity measurements can be used to determine the dissociation constant. 
According to Ostwald s law, K dissc is given as: 


Kdiss,c 


( A (cq) V 

V Aoo / 


■co 


MS 2 ) 


(Eq.4.20) 


or 


1 _ Cq-A(cq) _j_ J_ 

^t(co) Kdiss.c'^ OO Aoo 


(Eq.4.21) 


Measuring the molar conductivity at various concentrations, and plotting 1/A(c 0 ) versus c 0 ■ A(c 0 )> 
therefore leads to a straight line yielding A^ from the intercept, and finally K diss ;C from the slope. Note 
that this method has the major disadvantage that a large experimental error in the intercept causes an 
even larger error in K diss c -K diss c mzy therefore be determined more accurately ifAoo is measured 
independently, for instance by a clever combination of the molar conductivities of fully dissociated so- 
called strong electrolytes: 

Aoo (//Ac) = A(//CZ) + A (NaAc) - A(NaCl) = A (//+) + A(Ac“) (Eq.4.22) 
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So far, we have ignored the effect of interionic interactions, which cause a decrease in ion mobility with 
increasing ion concentration. For the weak electrolyte, this contribution is taken into account if in our 
formulation of K diss we switch from concentrations to so-called activities, i.e. 


K, 


diss 


O-Ac + 
a HAc 


f Ac ' c Ac '/// + ' c // + 
f HAc ' c HAc 


r 2 

/+ ^Ac "£// + 
f HAc ' c HAc 


f± 

f H Ac 


K, 


diss ,c 


(Eq.4.23) 


r 2 

J± is the mean-squared ionic activity coefficient, which decreases (from 1 at infinite dilution) with 
increasing ion concentration. f HAc is the activity coefficient of the uncharged acetic acid molecules, which 
is close to 1 and nearly independent of concentration. Note that, since / + 2 depends on concentration, 
Kdiss,c is not a constant! 


Finally, also the molar conductivity itself has to be corrected for interionic interactions, and Ostwald s 
law therefore has to be modified: 

A(c 0 ) = / A ■ a ■ Aoo (Eq.4.24) 

As we will show in the next section, according to Poisson-Boltzmann-theory f A is given by a square- 
root dependence on concentration: 

/a = 1 B • Vc 0 (Eq.4.25) 
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Example 4.1: 

Consider two aqueous solutions of a weak organic acid with concentrations 0.1 mol/L and 0.01 
mol/L, respectively. The ratio of the electric resistance for these two solutions is 1:3. Calculate 
the dissociation constant of the acid. 


Solution: This problem is related to the conductivity of a weak electrolyte. We therefore have to 
solve the following set of equations: 


A oc ■ Aqq , 



R~ 1 -A/l 


and 


= 


1 -a 


The ratio of Ohm resistances R yields, therefore, a ratio of dissociation constants: 

a(c=0.1mol/L) _ =0.1 mol/L) _ /fl -1 (c=0.1 mol/L)-A/l\ //fl -1 (c=0.01 mol/L)-A/l\ 

a(c=0.01 mol/L) A(c=0.01 mol/L) V 0.1 mol/L )'\ 0.01 mol/L J 


a{c=0.1mol/L) , 

f 3 ) 

i/I 

f 1 1 

a(c=0.01 mol/L) 

V0.1 mol/L) 

\0.01 mol/LJ 


We insert this result in K s and obtain: 

g _ a 2 -c _ a(c=0.01mol/L) 2 -0.01 Mol/L _ a(c=0.1 mol/L) 2 -0.1 mol/L _ (o.3- a(c=0 .01 mol/L)) 2 -0.1 mol/L 
5 1 -a l-a(c=0.01 mol/L) l-a(c=0.1mol/L) l-0.3-a(c=0.01mol/L) 


a(c=0.01 mol/L) 2 -0.01 mol/L _ ( 0.3-a(c=0.01mol/L)) 2 -0.1mol/L 


l-a(c=0.01 mol/L) 


l-0.3-a(c=0.01 mol/L) 


0.01 mol/L ■ (l - 0.3 ■ a(c = 0.01 mol/L)) = 
0.3 2 ■ 0.1 mol/L ■ (l — a(c = 0.01 mol/L)) 


0.01 - 0.01 ■ 0.3 ■ a(c = 0.01 mol/L) = 0.3 2 ■ 0.1 - 0.3 2 ■ 0.1 ■ a(c = 0.01 mol/L) 


or 


a(c = 0.01 mol/L) = °f ' 0ml Q ' Q1 = 0.04762 , and 

v ' 0.3 2 -0.1—0.01 -0.3 

a(c = 0.1 mol/L) = Q '°^ 762 = 0.01587 , and therefore finally: 


Ks 


a 2 -c _ 0.04762 2 -0.01 mol/L _ 0.01587 2 -0.1 mol/L 
1-a ~ 1-0.04762 _ 1-0.01587 


2.4 ■ 10" 5 mol/L 
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4.1.2 Interionic interactions and the Poisson-Boltzmann-formalism of the Debye-Huckel-theory: 

The Debye-Hiickel theory is a quantitative model allowing to calculate the mean-squared ionic activity 
coefficient /+ (see Eq. (4.23) ) based on pure electrostatic interactions and the Boltzmann probability. 
Qualitatively, we consider the interactions of a central positive ion with a surrounding cloud of negative 
ions (see fig. 4.5). Note that the overall net charge of the cloud has to be -1, since otherwise the principle 
of electroneutrality would be violated. 



Figure 4.5: The Debye-Huckel model of the counterion cloud 

The chemical potential of a 1,1-electrolyte, for example NaCl in aqueous solution, including interionic 
interactions, is given as 


M = M(id) + 2RT ■ In f+ 


(Eq.4.26) 


with the chemical potential in case interionic interactions may be ignored, that is, at infinite 
dilution of our electrolyte solution. Physically, the interionic interaction per single ion can be expressed 
as electrical work w to create a central positive charge in presence of the elecrostatic potential created 
by its surrounding ion cloud, or w = e ■ cp w (r = 0), where <p w (j = 0) is the electrostatic potential of 
the ion cloud at the position of the central ion. Note that this work is negative in value, i.e. this charging 
process leads to a gain in energy! Using this physical concept In f± , can simply be expressed as 


In f± 


V'-V'jid) _ e-(p w (r=0) 
2 RT ~ 2 RT 


(Eq.4.27) 
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Our problem how the ionic activity coefficient (and also the molar ion conductivity) depends on 
electrolyte concentration is quantitatively solved if we are able to calculate cp w (r = 0) as a function of ion 
concentration. To obtain a general expression for the potential of the ion cloud (p w (r), we consider the 
potential of the single ion cp 1 (r) (positive charge 1+) and that of the single ion screend by the surrounding 
counterion (negative charge 1-) cloud <p 2 ( r )- Note that we still limit our formalism, for simplicity of 
the resulting mathematical expressions, to a 1,1-electrolyte. Having derived the final expression for the 
ionic activity coefficient, we will finally present the general expression valid for i,j-electrolyte solutions, 
with i,j > 1. The electrostatic potentials for our 1,1-electrolyte are given as: 

= ST7 (E( 1' 4 ' 28) 

^W = i'r ex p(-^) 

£ is the permittivity of the solvent, and r D is the screening length of our sceened Coulomb repulsion, 
at this stage still a non-specified parameter. The electrostatic potentials given by Eqs. (4.28) and (4.29), 
respectively, are compared in fig. 4.6. Next, we will calculate the Debye screening length r D , which is 
directly related to the ionic activity coefficient /+ . 
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r/a.u. 


Figure 4.6: electrostatic potential of a single ion (blue), and a single ion including the counterion cloud (red) 


We start our calculation of r D with the Poisson-equation, one of the Maxwell equations connecting a 
radial-symmetric charge density profile Pi( r ) with the corresponding electrostatic potential <Pi(r): 


— . A. ( r 2 . d( Pi ( r ) 1 _ _ Pi ( r ) 

r 2 dr l dr ) s 


(Eq.4.30) 


Inserting our approach for the screened Coulomb potential ((Eq.4.29)) into the Poisson-equation (Eq. 
(4.30) ) leads to: 


<Pi(r) _ pj(r ) 
r D 2 £ 


(Eq.4.31) 


This equation can be solved if we consider also the Boltzmann equation (see Chapter 2, Eq.(2.15)) to 
express the local charge density of ions of species j at distance r from our central ion i i(Nj (r)). The 
local ion concentration Nj (r) depends on the ratio of electrostatic interaction energy to thermal energy 
k B T (with the Boltzmann constant k B = R/N A , R = universal gas constant, N A = Avogadro number, 
see also chapter 2) as: 


Nj(r) 

(Nj) 


= exp 


( 


e-(pj(r) \ 

k B T J 


(Eq.4.32) 


(Nj ) is the average concentration of ions of species). Note that at very high temperature, Nj (r) =(Nj), 
i.e. the concentration of ions of species j is constant all over the electrolyte solution, whereas at lower 
temperatures ions and counterions form a regular alternating spatial arrangement as sketched in fig. 4.5 
(= the ion cloud model). 
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Considering both positive and negative ions, we get a local charge density distribution as a function of 
distance r in respect to our central cation given as: 


e-cpjir) 
k B T . 


') - (N) ■ e ■ exp (Eq.4.33) 


Pi (r) = N + (r) ■ e — N_(r ) ■ e = ( N ) ■ e ■ exp ( 


Note again that, for simplification, here we have considered ions of charges +1 or -1, only, and 
(iV + ) = (N_) = (N ) (that is, the average particle number concentrations are identical for cations and 
anions). 

Usually, the electrostatic interaction energy is small compared to k B T (weak perturbation), and using 
a Taylor series expansion of the exponential function the local charge density can therefore be written 
simply as: 


Pi O’) = {N)-e-(l 


e-(Pi(r ) 

k B T 


■) —(N) -e-(l + 


e-(Pi(r ) 

k B T 


■) 



(Eq.4.34) 


To treat the general case of multivalent ions, it is better to express Pi (r ) in dependence of the ionic 
strength instead of in dependence of the average ion number concentration (N): 



(Eq.4.35) 


p is the density of the solvent, N A the Avogadro number, and m 0 a dimensionality factor (rrft = 
1 mole/kg).. The ionic strength itself is defined as following: 



(Eq.4.36) 


with rrtj the molality of ions of species j. As an illustrative example, consider 0.5 m NaCl: 


1 = y {0.5 • l 2 + 0.5 • (-1) 2 } = 0.5 


(Eq.4.37) 


With 


(N) = p-N A -m 0 ■ I 


(Eq.4.38) 


we get an average ion number concentration of 0.5 ■ N A per liter for both ionic species. 
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2 e 2 / \ 

If we relate the expression Pi(r) =- 6 ^ r ■ p ■ N A ■ m 0 ■ / (Eq.4.35) with the result of the Poisson 

equation < £i^l - _p^LL (Eq.4.31), we obtain the following expression for the ion concentration 
dependence of our screening length r D : 


r D 


2 


e-RT _ e-RT 

2-p-(e-N A ) 2 -I-m® ~ 2-p-F 2 -I-m 0 


(Eq.4.39) 


with the Faraday constant F = eN A — 96484.6 C/mole. 


Our final goal was to express not the screening length but the mean-squared ionic activity coefficient 
/ + 2 as a function of ion concentration (or ionic strength). For this purpose, we still have to show how 
<Pw ( r = 0) depends on the Debye screening length r D . The electrostatic potential of the ion cloud (p w (r) 
is simply given by the difference of the potentials <Pi(r) and <P2( r ) (see Eqs.4.28, 4.29), i.e.: 

(Pw (r) = q>2 (r) - <P 1 (r) = ^ f (exp (- - l) (Eq.4.40) 


For small distances r from the central ion, we can use a Taylor series expansion for the exponential in 
(Eq.4.40), and obtain: 


«’»' (r = 0) = sr'r ( 1 



_e _ 1_ 

4n£ r D 


(Eq.4.41) 
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Formally, the electrostatic potential of the ion cloud therefore corresponds simply to the Coulomb 
attraction potential of a single anion at interionic distance r D ! The electrostatic work W to charge one 
mole of cations from zero charge to charge e against this electrostatic potential of the counterion cloud 
is then given as: 


W = N A -tf e <P w (r — 0)dQ — — 


Na rQ=e 

4n£'rp ^0 


QdQ = 


Na 


4ns-r D 



F 2 

8ns-Ny[-rj) (Eq.4.42) 


In (Eq.4.42), the infinitesimal work to charge a single ion in respect to a constant electrostatic potential 
created by the surrounding ion cloud of counterions is given as dw = cp w (r = 0 )dQ, and therefore the 
total work per ion to achieve a final charge of e is w = e (p w (r = 0 )dQ. 


2 

Finally, the mean-squared ionic activity coefficient f+ of the electrolyte is given by considering both 
cationic and anionic species, leading to: 


In f± 


N + -W+N_-W 

(N + +N_)-RT 


F 2 

8ns -Na -r^-RT 


= —£■ — = — k ■ V7 = —k" ■ Vc 

rD 


(Eq.4.43) 


For aqueous solutions at T = 25°C, Eq.4.43 yields: 


In/+ = -0.509 ■ \z + ■ z_| ■ V7 


(Eq.4.44) 


This square-root dependence of the activity coefficient on ionic strength also leads to a square-root 
dependence of the molar conductivity of strong electrolytes in solution (an exact derivation is beyond 
the scope of this book), as described by the famous Vc-law by Kohlrausch: 

A (c 0 ) = Aoo - B ■ Vc 0 (Eq.4.45) 
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In summary, interionic interaction increasing with ion concentration therefore leads to a systematic 
decrease in the molar conductivity. Qualitatively, and as a conclusion to this mathematically quite 
elaborate section, we consider the interionic interaction of a central ion with its surrounding counter 
ion cloud via two simplified pictures (see figs. 4.7, 4.8): 

(i) relaxation of the ion cloud: 



Consider the migration of the positively charged central ion to the left, and the migration of the negatively 
charged ion cloud to the right according to the electrostatic attractions caused by the respective electrodes. 
As a consequence, a dipole is created, or charge neutrality is destroyed on a local level. This leads to a 
retardation force, i.e. the migrating anions are pulling back the cation towards the center of the counter 
ion cloud, and vice-versa. 


(ii) electrophoretic effect: 



Figure 4.8: electrophoretic effect: 


The frictional force is enhanced due to the fact that hydration shells of ions moving into opposite 
directions are in close contact. Therefore, the frictional force determined by the relative velocity difference 
is stronger than that of a single ion moving through solvent molecules which on average do not move 
in a specific direction. 
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To conclude this section about ion conductivity in electrolyte solutions, fig. 4.9 shows the concentration 
dependence of the molar ion conductivity of two strong acids and one weak acid in aqueous solution, 
illustrating the difference between weak and strong electrolytes. 
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4.2 The electrochemical potential and electrochemical cells 

4.2.1 Chemical electrodes 

A system consisting of a metallic rod immersed in an aqueous solution containing the corresponding 
metal cations shows the following chemical equilibrium: 

M z+ (aq) + ze~ Af (s) (Eq.4.46) 

Cations are migrating to the metal rod (or metal electrode), take up electrons from the metal and adsorb 
as metal atoms (reduction), or - the opposite reaction - metal atoms release z electrons to the metal rod 
and migrate into solution as metal ions (oxidation). Which process is favored depends in practice on the 
counter electrode, a 2 nd chemical electrode which closes the electric circuit (see fig. 4.1)! 

The reduction process causes the electrode to assume a positive charge, whereas the oxidation leads to 
a negatively charged electrode, and correspondingly a negative electric potential (electrode respective to 
solvent) A cp. In Eq. (4.46), (aq) means dissolved in water, and (s) means solid phase. 

For a given chemical electrode, its electric potential A<p, related to the difference in electric charge 
between the aqueous and the solid phase, is determined by the electrochemical potentials of the products 
and educts of the above redox reaction. In general, an electrochemical potential is defined as chemical 
potential plus electrostatic interaction energy, i.e. 

JTi= Hi + N A z t ecp = /it + z t Fcp (Eq.4.47) 

Our electrochemical equilibrium therefore is given as: 

Fm(s) ~ F-M z+ ( < aq) z Fe~{s) (Eq.4.48) 

with Fm(s) the chemical potential of the uncharged metal electrode, and the two electrochemical 
potentials given as: 


FM z+ (aq) FM z+ (aq) zF<p(uq) 


(Eq.4.49) 


Fe ( 5 ) Fe ( 5 ) F<p(s) 


(Eq.4.50) 
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In Eqs.4.48 and 4.49, z is the charge number of the metal ions in solution, for example for Cu 2+ z-2. 
The electrostatic potential difference A cp = cp(s) — <p(aq)between the solid and the aqueous phase in 
electrochemical equilibrium is then given by inserting Eqs.4.49 and 4.50 into Eq.4.48, and resolving 
for A cp: 


A (p = (p(s) - (p(aq) = — 






or 


A (p = <p(s) - (p(.aq) = ^ ■ {iU M z + (a(/) 0 + RT ■ In a M z+ + zfi e - (s) - q M(s) ] (Eq.4.51b) 

M M z+ (aq ) 0 is the chemical potential of the metal ions in aqueous solution at standard conditions, i.e. 
concentration 1 mole/kg solvent excluding interionic interactions (pseudo-ideal solution, not a real 
system since at such high concentration, as shown in the previous section, Poisson-Boltzmann-theory 
predicts strong interionic interactions). a M z+ is the activity of the metal ions in the solution. Per definition, 
all chemical potentials of solid phases are constant and therefore defined as zero (or, more accurately, 
included in the standard electric potential difference of the respective electrode). In conclusion, the 
electrode potential therefore depends on the materials and the ion concentration as: 

A cp = Acp® + — ■ In a M z+ (Eq.4.52) 


with the electric standard potential of the chemical electrode A<p 0 given as 

A(p0 = 7 f' k« z+ (aq) 0 + z Re-(s ) “ Rm(s)} (Eq.4.53) 

i.e. corresponding to an idealized aqueous solution at standard conditions, i.e. concentration 1 mole/kg 
solvent excluding interionic interactions, therefore activity a M z+ = 1. 

More general, for any redox reaction of type OX + ze~ ^ REDwe can define an electrode potential 
difference as: 


A(p = A(p 0 + --ln^ 

zF a RED 


(Eq.4.54) 


The electrode itself is formulated as RED, OX/M. Note that only components in the solute phase, typical 
aqueous solution, or gases will explicitly contribute to the concentration-dependent term in the logarithm. 
An alternative formulation is based on the chemical equilibrium, i.e. 


RT 

RED,OX/M ~zp 


A cp — A <p 0 nz7n / n/r In 


a RED 


&OX ' a e~(sY 


RT 

RED,OX / M ~zp 


= A(p 0 nI , n ^ v/ „ - ^-ln K 


(Eq.4.55) 
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Note that, per definition, a e -( s ) = 1. 

To illustrate the meaning of this general equation (4.55), let us consider in detail examples for the most 
important types of chemical electrodes: 

i. gas electrodes 

A famous very important example is the hydrogen electrode H + (aq)/H 2 (g)/Pt, which is also used as the 
reference system, i.e. A<p 0 p t /// 2 ( < g)/// + (aq) electrode is set to zero. Note that electrode potentials 

can only be measured within a closed electric circuit, that is, in reference to each other. Therefore, there 
exists no stand-alone absolute electrode potential, but all data are based on the H + (aq)/H 2 (g)/Pt reference 
system, leading to the electrochemical potential series found in all textbooks of Physical Chemistry. In 
this system, the electrochemical standard potential of any given chemical electrode corresponds to the 
electric potential difference one would measure in respect to the hydrogen electrode, using standard 
conditions, i.e. p = 1 bar, T = 298.15 K and activities of ions and gases = 1. 
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Nobler metals, like copper, silver etc., in respect to this hydrogen standard show a positive electrode 
potential, whereas less noble metals, like iron, zinc etc., show a negative potential. This means that, at 
standard conditions (ion concentration in the aqueous phase 1 mole/kg, pressure of gases contributing 
to the electrode potential 1 bar), a zinc electrode combined with the hydrogen electrode to one electric 
circuit will show spontaneous oxidation of zinc and reduction of hydrogen ions, whereas a copper 
electrode combined with the hydrogen electrode will show reduction of copper ions and oxidation of 
hydrogen gas. 

The H + (aq)/H 2 (g)/Pt electrode is based on the following redox equilibrium, using a Pt wire to transfer 
the electrons between hydrogen gas and protons: 

H + (aq ) + e~(s) 1 / 2 H 2 (g) (Eq.4.56) 

The electrochemical equilibrium is defined via the respective chemical or electrochemical potentials as: 

V 2 ^2 (s) = ^H+Caq) + Me-(s) (Eq.4.57) 

According to the general scheme, the electrostatic electrode potential difference then is given as: 

A<p = (pis') - (piaq) = i ■ + RT ■ In a H+(a(?) + zp e - {s) - Ph 2 ^) 0 ~ RT ■ ln/tf 2 ( 5 )}(Eq.4.58) 

Here, fn 2 (g) * s fugacity of the hydrogen gas (or effective pressure), defined as 

y-p H 

fn 2 {g) = (Eq.4.59) 

withy < 1 an activity coefficient (taking into account the difference between ideal and real gas, see 
chapter 2), and P 0 the standard pressure 1 bar. Combining all concentration-independent terms then 
leads to: 


A<P A( P 0 pt/H 2 (g)/H+(aq ) 



a H + (aq) 
f 0.5 
J # 2 ( 0 ) 


ii. electrodes of 2 nd type 


(Eq.4.60) 


In case the cation concentration is determined via the solubility of a nearly insoluble salt, using 1 mole/kg of 
the counterions as the standard condition, an electrode of 2 nd type is obtained. Examples are the Calomel- 
electrode Cl~(aq)/Hg 2 Cl 2 (s > )/Hg(V) or the silver/silver chloride electrode Cl~(aq)/AgCl(s)/Ag( < s ). 
For the later, the following redox equilibrium is the basis of the electrodes potential: 


AgCl(s) -he (s) Ag(s ) + Cl (aq) 


(Eq.4.61) 
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We use the general expression for the electrode potential of any chemical redox electrode 

A cp = A cp 0 + ^ ■ In (Eq.4.62) 

zt a RED 

and consider that on the “OX-side” there are only solid components with activity 1, whereas on the 
“RED-side” the only non-solid component is Cl~(aq). Consequently the electrodes potential is given as: 

- _ a 0 , RT_ i 1 _ A 0 AT , 

l\(p - !\(p C r{aqVAgCl{s)/Ag{s) + f acl - ~ ^ Cr{aq)/AgCl( < s)/Ag{s)~ F ' [U a Cl~ (Eq.4.63) 

This electrode of 2 nd type provides a nice example how the solubility constant of can be calculated from 
standard electrode potentials. The CZ“(aq) /AgCMjs) /Ag(s)-electrode can formally be considered also 
as a simple Ag + (aq)/Ag(s ) -electrode, where the effective concentration of Ag + -ions is determined 
by the solubility product, 

K L,AgCl - a Ag + ' a cr (Eq.4.64) 

The electrode potential therefore can be expressed in two different ways: 

Acp = Acp\r {aqV A g ci{s)/AgM ~ T ' ln acl ~ = A( P 0 Ag + (a q )/Ag(s ) + T' ln a Ag + (Eq.4.65) 

Eq. (4.65) yields for the solubility constant K L>AgC i : 

A( P Cr(a q )/A g Cl( S )/Ag(s) ~ A( P Ag+(a q )/A g (s) = T' ( ln + ln a Cl~) = T ' ln K W (Eq.4.66) 
With A (P*ci- {aq)/Ag ci { s)/Ag{s) = °- 22 V and A( f* 0 Ag + (a q )/Ag(s) = 0 ' 80 V ( fr0m the SerieS ° f Standard 

electrode potentials) we obtain K LAgC i « 1.5 ■ 10 -10 mole 2 /kg 2 . 
iii. redox electrodes 


As we have seen, every chemical electrode is based on a redox reaction. However, in practice you call 
an electrode a redox electrode if both the OX-species and the RED-species are found in the solute 
phase, typically an aqueous solution. Examples are the Fe 2+ (aq), Fe 3+ (aq)/Pt(s)- electrode or the 
Hydrochinon^aq), Chinon(aq)/Pt(s) -electrode. Note that there is no phase boundary between the 
OX and RED species, wherefore the / is replaced by a. The respective redox reactions and electrode 
potentials are given as: 


Fe 3+ (aq) -he (s) ^ Fe 2+ (aq ) 


(Eq.4.67) 


A(p A (p f e 2 + ^ ag ^p e 3 + ^ ag yp t ^ 


(Eq.4.68) 
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and 


Chinon(aq) + 2 H + (aq) + 2e ( 5 ) Hydrochinon^aq ) 


(Eq.4.69) 


A cp — A cp 


0 

Hydroc hinon ( aq),C hinon ( aq)/Pt(s ) 


_|_ RT_ a a Chinon (aq)' a // + (gq) 
2 F CL Hydroc hinon ( aq ) 


(Eq.4.70) 


Note that the Chinon/Hydrochinon (or Chinhydron) - electrode can also be used to measure pH-values 
of aqueous solutions, since its electrode potential depends on ci H +^ aq > } 2 . This electrode is much better 
suited for this purpose than the Pt/// 2 ( < g)/// + (aq)-electrode, where the potential depends on the 
hydrogen gas pressure which cannot accurately be adjusted without experimental difficulties . 


iv. Diffusion potentials (or membrane potentials) 


Another means to create an electric potential, without a chemical redox reaction, is based on osmotic 
pressure and a semipermeable membrane. Let us, for example, assume that two chambers are separated 
by a semipermeable membrane which only allows the smaller cations of a given salt to permeate, and 
the left chamber contains a lower concentration of salt in aqueous solution (see fig. 4.10). The difference 
in osmotic pressure then causes the smaller cations to migrate through the membrane from the right 
chamber of higher concentration to the left chamber, until the osmotic pressure difference is balanced 
by the electrostatic interaction, or, in other words, the electrochemical potentials of the cationic species 
in both chambers are equal. 



Deloitte 




Discover the truth at www.deloitte.ca/careers 


© Deloitte & Touche LLP and affiliated entities. 


125 

Download free eBooks at bookboon.com 





Basic Physical Chemistry 


Electrochemistry 


a,<p a > 0 

P.Vp < o 

O 

© ® 
o 

("■ \© © 
o 

© 

G© © 

°© ®o 

O °o @ ° 

.y...o.. 


Figure 4.10: Membrane potential 


Ma = Ma + F ■ <p a = Up = Up + F ■ <p p 

(Eq.4.71) 

h(p = (p a - (pp = Up - Ha = Y ■ ln^ > 0 

(Eq.4.72) 


Note that the concentration is not changing much by the ion migration, and therefore the electric potential 
difference is given by the ratio of the original salt concentrations. This type of membrane potentials is 
the basis of our neurons and signal conductivity in our bodies. Semi-permeability of the neurons here 
is achieved via special proteins embedded within the cellular membrane which are called ion channels. 

In this context, one should note that two electrodes are often connected to an electric circuit via a so- 
called salt bridge, which contains an electrolyte not participating at the redox reaction but just closing 
the circuit by ion migration. To avoid charge separation and diffusion potentials, which in this case 
would partially compensate the electrode potential difference (and therefore would lead to a loss in the 
performance of our chemical battery), the cations and anions of the salt bridge should have identical 
ion mobilities. This is the case for the salt KC1, which therefore is used in such salt bridges. In the next 
section, we consider such an electrochemical circuit consisting of two electrochemical electrodes, the 
Galvanic cell. Fig. 4.11 (see below) shows how the electric circuit is closed in this case by the salt bridge. 
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4.2.2 The electrochemical Galvanic cell 

If you combine two of the chemical electrodes we discussed in the previous section to a closed electric 
circuit, you obtain a Galvanic chain or electrochemical cell. As an electrochemical battery, this setup 
is used to convert the chemical energy of a spontaneous redox reaction into electric energy. The 
spontaneous process can also be reversed by applying an external electric voltage (electrolysis chamber). 
For the spontaneous process corresponding to a chemical battery, per convention the electrode showing 
oxidation is placed left, and the electrode showing reduction is placed right. To close the electric circuit, 
typically a salt bridge containing an aqueous solution of KC1 connects the two electrodes or half cells, 
as already discussed. The Gibbs free enthalpy of the redox reaction of the electrochemical cell, and the 
electrochemical potential difference (or electromotive force E) of the Galvanic cell are related as (see 
also definition of AG as electric work in chapter 2 !): 

A R G = -z-F -E (Eq.4.73) 

In general, the redox reaction of the cell Pt(s)/REDi, OXJ/RED r , OA r /Pt(s) is given as: 

REDi + 0X r —> 0X t + RED r (Eq.4.74) 


Here REDi,RED r , are the reduced species in the left and right part of the Galvanic cell, OXi,OX r 
correspondingly the oxidized species, and / / the salt bridge connecting the two chemical electrodes or 
half cells. The electromotive force of this cell is then given via the difference of the electrode potentials, i.e. 


E = A (p r — A cpi = A 


. RT . CLOXy 

+ — ■ In- 

zF a RED r 


-Vz+f ■ In 


a ox l 

a RED : 


= g 0 —-• (Eq.4.75) 

zF a OX r ' a REDi 


In chemical equilibrium, our Galvanic cell has no electromotive force, and the respective activities of 
the components of the chemical redox reaction are determined by the chemical equilibrium constant 
K. Therefore, we obtain 


E = O = E 0 - — -InK (Eq.4.76) 

zF 

With the general relation for chemical equilibrium (see Chapter 2, A R G* = —RT ■ In K (Eq.2.146)), we 
finally get: 

A r G* = —z • F • E® (Eq.4.77) 
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To illustrate the concept of a Galvanic cell, we consider the Danielle-element as a first example of a 
chemical battery. This Galvanic chain connects a zinc-electrode and a copper electrode via a salt bridge, 
and therefore is formulated as Zn(s) / Zn 2+ (aq)// Cu 2+ (aq)/ Cu(s ). The cell is sketched in fig. 4.11: 
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In this electrochemical cell, we find the following chemical redox reactions: 

left chamber: Zn -> Zn 2+ + 2e~ 

right chamber: Cu 2+ + 2e~ -> Cu 


cell reaction + Cu 2+ —> Zn 2+ + Cu 


The less noble zinc therefore is oxidized, while the nobler copper is reduced. The electromotive force E 
of this cell then is then given as: 

E = h<Pcu 2+ /Cu - ^ ( Pzn 2+ /Zn ~ ^Vcu^/Cu® ~ ^Vzn^/Zn* ~ a" 2 + (Eq.4.77) 

As a 2nd example, we consider an electrochemical cell where a hydrogen electrode and a silver/ 
silver chloride electrode are immersed in a single chamber containing dilute hydrochloric acid, 
Pt(s)/H 2 (g)/HCl(aq)/AgCl(s)/Ag( < s ). In this case, we expect the following chemical reactions at the 
respective electrodes: 

left chamber: 1/2 H 2 (g) -> H + (aq) + e~ 

right chamber: AgCl(s) + e~ -> Cl~(aq ) 


cell reaction 1/2 H 2 (g ) + AgCl(s) -> // + (aq) + Cl (aq) 


The noble metal silver therefore is reduced, whereas the less noble hydrogen is oxidized. The corresponding 
electromotive force E of this cell is given as: 


E = A(pcr/AgCl/Ag ~ A <Ph+/h 2 = Cl~ / AgCl/Ag* ~ b<pH+/H 2 * ~ T ' ln (Eq.4.78) 

Since the standard electrode potential of the hydrogen electrode is the has is of the electrochemical potential 
series and defined as A cp H +/ H ® = 0 V, this electromotive force E is simply given as: 


F — Am 0 RT .] n a //+' acr 

t - b<Pcr/AgCl/Ag ~Y ln / 0.5 


(Eq.4.79) 


with A (Pcr/Agci/Ag 0 = +0.2223 7. 
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In practice, electrochemical cells are used in physical chemistry mainly for three different purposes: 


i. determination of the mean ionic activity coefficients f ±: 


Consider for illustration our 2 nd example of an electrochemical cell Pt(s)/H 2 ( < g)/HCl( < aq)/ 
AgCl(s)/Ag(s ). The electromotive force E of this cell can be expressed as: 


E - A(Pcr/AgCl/Ag 0 ~ ~ ‘ 

ln/± 


Q'll + 'Q'Cl 

—— 

fu 2 


A 0 RT 

^Pcr/AgCl/Ag ~ V ' 


In 


™ H + , ™-Cl—f± 


Vh 2 


— EMK th -— 


2 RT 
F 


(Eq.4.80) 


In /+ is obtained if you compare the experimentally determined electromotive force E with the 
theoretically expected value E th calculated from the standard electrode potential A (pcr/AgCl / Ag 0 and 
the concentration of the hydrochloric acid (in mole/kg !). 


ii. a second type of application is the electrochemical measurement of pH-values, as discussed 
before, using a reference electrode and a sensor electrode whose electrode potential depends 
on the concentration of H + . 

iii. more general, one can determine the concentration of many ionic species in aqueous 
solution if one uses the appropriate chemical electrode. This quantitative analytical method 
is used, for example, in potentiometric titration, where you monitor the electromotive 
force E in dependence of added amount of analyte. One example could be the titration of a 
Pe 2+ -solution with Ce 4+ , using a aCZ - /AgCl/Ag -electrode as reference. If the reaction 
flask is an electrochemical cell containing a Pt-wire, we then formally measure via the 
electromotive force E the electrode potentials of the following redox electrodes depending 
on the present species: either Fe 2+ (aq),Fe 3+ (aq)/Pt(s ) or Ce 3+ (aq ), Ce 4+ (aq)/Pt(s). 

Note that, to accurately determine ion concentrations via measurements of the electromotive force, the 
electrochemical cell has to be used without any current to avoid electrolytic processes at the electrodes 
during the measurement, which would change the concentration of interest. In practice, this can be 
achieved either by compensating the electromotive force E with an external electric potential (so-called 
Poggendorf compensation), or by using a high resistance within the electric circuit to suppress the 
charge transport. 
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Example 4.2: 

Consider an electrochemical battery based on two copper electrodes with concentrations 
c Cu 2 + = 1 Mol/L and c Cu 2 + = 1 ■ 10 -5 Mol/L. The volume of the electrolyte solution of each 
electrode is 1 L, and the electric current of this battery is 50 mA. Determine the maximum time 
this battery will last, and the variation in electric voltage with time. 

Solution: This battery is “empty” if the copper ion concentrations in both chambers are equal. 
In this case, formally about 0.5 Mol of copper ions have been transferred, corresponding to the 
exchange of 1 Mol electrons. The charge of this Mol electrons is directly given by the Faraday 
constant, Q = F = 96484 C/Mol. Therefore, this battery will last for: 

O 96484 

Q = I ■ t => t = J = 5 = 1929680 5 = 22.33 d 

Next, let us calculate the electromotive force E in dependence of time. In general, the voltage 
of our battery is given as: 

F = FL . /„ C i( c “ 2+ ) 

2 F c 2 (Cu 2+ ) 

The concentration changes by the following rate per second: 

, Q 0.05 C 0.05/F . _ 0.025/F ,, , „ 2 + 

I = - = -= —— Mol e = - — Mol Cu l+ 

t IS 1 5 1 5 

, since we need two electrons per copper ion. 

Therefore, we obtain the following general expression for the time dependence of the voltage 
of our battery: 

9 _i_\ 0.025 „ 

RT 7 Ci(Cu 2 + ) RT , 1 - 

— 7 -vi ^ «= — In- L 


E — — In 


2 F ‘'"'c 2 (Cu 2 +) 2 F ‘“io-5+M2L t 

F 


The following table gives us an impression how rapidly the voltage decays: 


Time/ s 

0 

1800 

3600 

7200 

18000 

36000 

Voltage/ V 

0.145 

0.0966 

0.0879 

0.0792 

0.0677 

0.0589 


Note that, although this battery will last for more than 20 days, the voltage has reached 50% 
after less than 5 hours (18000 s, see table)! The performance can be strongly enhanced, i.e. the 
voltage kept nearly constant and the maximum duration doubled, if one “catches” the copper 
ions released in the chamber of lower copper ion concentration either by complexation or by 
precipitation. 
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We conclude our chapter on electrochemistry by presenting the chemical redox reactions for some of 
the technically most important electrochemical batteries: 

i. lead accumulator 

oxidation: Pb(s) -> Pb 2+ (aq ) + 2e~ 

reduction: Pb0 2 (s) + 2e~ + 4 H + (aq) Pb 2+ (aq ) + 2H 2 0(l) 

cell reaction Pb(s ) + Pb0 2 (s ) + 4 H + (aq) -> 2Pb 2+ (aq ) + 2H 2 0(l ) 

This cell reaction underlines the importance of the pH of the sulfuric acid for the performance of the 
lead acumulator! 

ii. Ni-Cd-accumulator 

oxidation: Cd(s ) + 2 OH~{aq) -» Cd(OH) 2 (s) + 2e~ 

reduction: 2 NiO(OH)(s) + 2H 2 0 + 2e~ -+ 2 Ni(OH) 2 (s) + 20H~(aq ) 

cell reaction Cd(s) + 2Ni0(0H)(s ) + 2tf 2 0 ^ Cd(0//) 2 (s) + 2Ni(OH) 2 (s ) 
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5 Introduction to Quantum 

Chemistry and Spectroscopy 

In this chapter, I will try to briefly introduce the concept of quantum chemistry, which is essential to 
understand spectroscopic methods, an analytical tool extremely important in chemistry. Some common 
spectroscopic techniques will therefore also be addressed in this chapter. Note that this chapter can serve 
as a rather superficial introduction to quantum chemistry only: especially, I refer from any detailed 
mathematical representation of the topic for the reason of readability (and lack of space), and the quantum 
chemical concept of the chemical bond is also not treated here. 

5.1 Models of the atom 

Historically, our picture for the nature of matter developed from the rather philosophical description of 
the ancient Greeks (atoms = simple small geometrically different objects, like cubes, spheres etc.), via 
Bohr’s model of negative charges migrating on concentric circles around a positively charged nucleus 
like planets orbiting the sun, finally to a new physical interpretation called quantum mechanics. Let us 
briefly repeat the most important historic steps of this development: 

Demokrit, 5 th century BC, stated that atoms (from Greek: ocGopocr = inseparable) are the smallest parts 
of any matter, and differ in geometrical shape and color. It took a very long time until Dalton, in the 
19 th century, interpreted the law of constant proportions to the conclusion that some matter consists of 
small groups of different atoms, called molecules. For example, the smallest parts of the substance water 
are molecules which consist of 2 atoms of hydrogen and 1 atom of oxygen. 

From the late 19 th century, our atomic models changed comparatively rapidly, i.e. within several decades, 
based on new experimental results. This underlines the importance of the actual method of scientific 
research: develop any model only based on experimental results, as well as check existing models based 
on their predictions of experimental results! First, the discovery of the negatively charged electron within 
the cathode-radiation, and the experimental prove of its particular character within the Wilson cloud 
chamber (see fig. 5.1), lead to the Thomson model visualizing atoms analogous to a cake, where the 
negative charges are embedded like resins within a positively-charged lump of comparatively low mass 
density (see fig. 5.2) (Thomson 1897). 
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Figure 5.1 : Cathode radiation and deviation of the electron beam in an electric field proofing the negative charge of electrons 
(top), and proof of its particular character (a (3-source emits electrons!) via the Wilson cloud chamber (bottom) 



Figure 5.2: The Thomson model of the atom (electrons like "resins in a cake"). Note that 
the "cake" itself (grey area) is positively charged to compensate for the electrons charge, 
since the atom itself is charge-neutral. 

Based on the low mass density of the atoms as well as the small size of the electrons, according to the 
Thomson-model one would expect that, if you direct a ray of small positive charges called a-radiation 
(= He 2+ -ions) onto a thin foil of gold atoms, the projectiles should pass through the foil more or less 
unhindered (Geiger and Marsden experiment, 1909). Although this was true for most of the a-particles, 
the experiment also showed that a small fraction of the projectiles was reflected by the atoms of the gold 
foil (see fig. 5.3). This was in contradiction with the Thomson model, and lead to a new model suggested 
by Rutherford, who postulated that the positive charge of the atom should be allocated within the center 
of the atoms within a tiny highly massive core (see fig. 5.3). 
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screen 


Figure 5.3: Geiger-Marsden experiment (top), and interpretation of the experimental results via a new atomic model 
suggested by Rutherford (bottom). 


The problem of the new Rutherford model was to explain why the negatively charged electrons obviously 
are not collapsing into the positively charged atomic core, i.e. what keeps the volume of the atom stable, 
since the atom should mainly consist of “nothing” as proven by the fact that most of the a-projectiles 
in the Geiger and Marsden experiment passed through the gold foil unhindered. Another experimental 
finding not explained by the Rutherford model are spectroscopic measurements showing that atoms 
absorb light of specific wave lengths only Niels Bohr therefore suggested in 1913 a new atomic model 
to explain the spectroscopic observations. 

If we assume that the absorption of light changes the state of the electrons, we may conclude that an 
electron bound within an atom only assumes certain levels of energy. Bohr developed a new model 
still based on classical physics, which postulated that the electrons are moving with constant velocity v 
around the positively charged core of the atom on stable circular orbits of radius r, like planets moving 
around the sun. In this case, the Coulomb attraction between the electrons and the core is compensated 
by the centrifugal force: 

6 ^ m 

F Coulomb ~ ~ ( 4 7rf0 ). r 2 — —Fcentr. = ~ (Eq.5.1) 
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The kinetic energy of these moving electrons is then given as E kin = 1/2 ■ m ■ v 2 . Note that, according 
to the Maxwell equations, such an electron moving within the electric field originating from the 
positively charged atomic core should lead to the continuous radiation of an electromagnetic field. As a 
consequence, the electron would also continuously loose some of its kinetic energy, which finally should 
lead to the collapse of the atom. Bohr “solved” this problem, and thereby also provided an explanation 
for the experimental results from spectroscopy that an electron bound to an atom only assumes specific 
energies (or specific orbits), by postulating that the electron is confined to specific orbits r n fulfilling the 
following relation concerning their rotational momentum: 

|f| = |f X p| — m ■ v ■ r n = n ■ y- (Eq.5.2) 


with n an integer (1, 2, 3, ...), and h = 6.626 ■ 10 -34 / ■ s Plancks constant. With this boundary 
condition, you obtain the following allowed values for the energy of the electron encircling the core of 
the hydrogen atom, i.e. charge = +le: 


j-, _ m e -e 4 1 _i o o 

bn ~ ~ (47 T £ o ) 2 '2 W(2tO)2 ' n -2 ’ n - L ’ Z ’ 6 ’' 


(Eq.5.3) 


m e is the mass of the electron, 9.109 ■ 10 -31 kg. The energy difference between two of these orbits, and 
therefore the frequency of light necessary to lift an electron from one inner orbit to a larger one of higher 
energy, then is given as: 


A E = E n2 - E ni = 


m e -e 4 /I 1 \ _ , 

(47T£o) 2 ‘2-(/i/(27z:)) 2 \n 2 2 ni 2 ) ^ 


(Eq.5.4) 


Note that, if the electron is lifted to defined energetically excited states by, for example, external electrical 
energy (see for example gas lamps), it will go back to its original energetic state while emitting light of 
specific wavelengths, corresponding in energy to the differences given by Eq.5.4 (emission spectrum of 
atoms, see fig. 5.4). 


136 


Download free eBooks at bookboon.com 




Basic Physical Chemistry 


Introduction to Quantum Chemistry and Spectroscopy 


A energy of the electron 




Figure 5.4: Energetic transitions responsible for the hydrogen emission spectrum, and the Bohr atomic model 
(red, blue arrows indicating the Balmer-series found in the visible regime of the electromagnetic spectrum) 
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Emission (or absorption, see further below) spectroscopy therefore can be used to quantitatively measure 
the differences of the discrete energetic states of atoms (or molecules, see below). Bohr could show that 
quantum physics, a new concept introduced first by Plack( 1900) and Einstein (1905), and classical physics 
are in agreement at large quantum numbers n (principle of correspondence). However, at small quantum 
numbers Bohr encountered certain difficulties to match his model onto experimental observations on 
a sound mathematical and physical basis. This problem could only be resolved by a completely new 
theoretical approach, modern quantum mechanics. 

Before we approach the concept of quantum mechanics (which was formulated in the mid-1920s by 
Werner Heisenberg, Max Born and Pascual Jordan (matrix mechanics), Louis de Broglie and Erwin 
Schrodinger (wave mechanics)), let us briefly sketch the nature of light as it is described within the 
classical physical model. Classically, light was considered as an electro-magnetic wave quantitatively 
described by the Maxwell equations and migrating with lights speed (c = 3 ■ 10 8 m/s ). The oscillating 
electric field amplitude vector of light of wavelength A migrating, for example, in x-direction, is given as: 

E(x,t) = E 0 ■ cos (^p - ■ t) (Eq.5.5) 

The electric field amplitude is showing periodic oscillations both as a function of space and time, and 
the oscillation frequency is given by the light velocity and the wave length as: 

v = \ (Eq.5.6) 

To conclude this section, note that at this stage of our description of the atom still based on classical 
physics, we encounter at least two fundamental problems: first, there is no physical reason why an electron 
should be limited to specific orbits only (early quantum theory postulated by Niels Bohr). Second, within 
the classical physical picture the moving electron should loose energy and therefore collapse into the 
atomic core, as already mentioned above. As we will show in the next section, a new concept of physics, 
modern quantum mechanics, is necessary to solve our problem. Interestingly, this new concept is also 
relevant in respect to the nature of light: classically, as we have just described, light is an electromagnetic 
wave. However, Einstein has discovered in 1905 in his famous publication explaining the photoelectric 
effect that light also must have a particle character. To resolve these problems, a new concept of physics 
was introduced in the 1920s (Heisenberg, Schrodinger, 1925), quantum mechanics, which attributes a 
wave character to everything, even solid matter. 

5.2 The wave character of matter, or the wave-particle-dualism 

Analogous to water waves, in some experiments small moving particles also show diffraction phenomena 
typical for a wave character. The general wave equation describes a harmonic oscillation of a physical 
quantity ^(x, t) with space and time, i.e. 

0 = % ■ sin (If- - (Eq.5.7) 
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In eq.5.7., T'q is the amplitude or maximum value of our physical quantity, A the wavelength and r the 
duration of one oscillation. If a circular wave, for instance light, is passing through two neighboring slits 
of adequate separation and size, these slits act as origins of two new circular waves with well-defined 
phase difference, and a screen behind the slits shows a regular interference pattern of dark and bright 
stripes (see fig. 5.5). On the other hand, if a beam of light is shining through a single small slit the screen 
does not simply show a sharp image of the slit, but also light of lower intensity on the edges, which looks 
as if the light is partially twisted while passing through the slit (= diffraction). 


\ 


light source 


J 


A 

h 




interference pattern visible 
on a white screen 


Figure 5.5: Diffraction of light wave on a two-slit setup 


Diffraction and interference are experimental phenomena proving the wave character of light. 

However, there exist experimental results which can only be interpreted if light is consisting of individual 
particles called photons, the photo-electric effect discussed 1905 by Albert Einstein. The photoelectric 
effect is the experimental result that electrons released from metal by incident light are assuming a kinetic 
energy (or velocity) independent of the light amplitude, but only dependent on the wavelength of light. 
On the other hand, the number of released electrons increases proportional to an increase of the power 
(or amplitude) of the light. These results can only be understood if one single light particle (= photon) 
releases one single electron, and the energy of the electron is depending on the energy of this photon, 
which itself depends on the wavelength or frequency of light. The number of released electrons then 
depends on the number of incident photons, which is given by both the power of light and the frequency. 
The energy of a single photon is given as: 

E — h ■ v — h ■ j (Eq.5.8) 
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Figure 5.6: The photoelectric effect, experimental setup (left) and resulting photo electric current (right) 


The experiment itself, as sketched in figure 5.6, works as follows: the electrons, released by incident light, 
are kept in the cesium plate if their energy is too low to migrate to the counter electrode which is kept 
at an adjustable electric potential (or voltage). Thereby, this potential provides a direct experimental 
measure for the energy of the photoelectrons, whereas a measurement of the electric current at lower 
potentials provides a measure for the number of released electrons. In this experimental setup, the 
wavelength and power of the incident light are varied, and the electric current in dependence of the 
voltage is measured. The electric current divided by the elementary charge then directly corresponds to 
the number of incident photons per second, whereas the electrostatic energy e • U plus the ionization 
energy of the metal per atom, typically cesium, corresponds to the energy of a single photon. 
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In conclusion, light behaves, depending on the experiment, like a wave or a particle, a phenomenon 
which is called wave-particle-dualism. This behavior is also found for objects which in classical physics 
typically were regarded as particles, like electrons or even particles of much larger masses. Davisson and 
Germer (1927), for example, showed that an electron beam hitting a small gold foil or a nickel crystal in 
reflection, that is, a periodic lattice on an atomic length scale, shows interference phenomena. 



Figure 5.7: Diffraction of an electron beam - the Davisson-Germer experiment 
and Bragg diffraction 

The Davisson-Germer experiment can be explained by the constructive interference of two electron waves 
reflected at the outermost atom layer and an inner layer of the nickel crystal, respectively. The inner 
reflected beam migrates a longer distance, and constructive interference occurs if this longer distance is 
a multiple of the light wavelength, leading to the famous Bragg law of diffraction (see fig. 5.7 and Eq.5.9): 

n ■ X — 2d ■ sin 9 (Eq.5.9) 

So far, we have learned that, depending on the actual experimental observation, light or even matter 
either exhibits particle or wave character. This fundamental problem of the wave-particle dualism can 
be resolved if we introduce the concept of standing waves. For example, an electron migrating around 
an atomic nucleus on a stable orbit can be regarded as a resonant standing wave, like the strands of a 
guitar. Note that for one-dimensional standing waves the amplitude at the edges have to be zero, limiting 
the allowed wave lengths (see fig. 5.8), and therefore leading to discrete energy levels or quantization of 
energy, as we have seen already in the Bohr model, where classical physics could provide no satisfactory 
explanation for the quantization. 
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Figure 5.8: 1-dimensional standing waves 


The description of very small particles as a standing wave will lead us directly to the mentioned new 
concept of physics, quantum mechanics. For a particle which is more or less localized, we usually have 
to use not a single standing wave but packages of interfering waves. These wave packages provide an 
explanation of Heisenbergs relation A p x ■ Ax > h, i.e. it is not possible to simultaneously determine 
the momentum Vx and the position x of any moving object with a (combined) accuracy A p x - A x < hi 
Another important aspect of the wave description of moving particles is the de-Broglie equation: 

V = h/A (Eq.5.10) 


This important equation relates the momentum p (= particle quality) to the wavelength A (= wave quality), 
and therefore providing an excellent conclusion to the wave-particle dualism. Here, it should be noted 
that especially for particles moving with velocities close to the speed of light c, the momentum has to 
be calculated in a relativistic manner, i.e. 


m-v 



(Eq.5.11) 


Note her that for photons of mass zero, we then get formally =0/0! However, from the general de 
Broglie equation (Eq.5.10.) we know that also for the photon p = h/A\ 


Let us go back to Heisenberg’s uncertainty principle: as the two extreme cases, one can either determine 
exactly the momentum p x or the position x, and correspondingly the position or the momentum have to 
be totally undefined (see fig. 5.9): in the first case, the moving particle is described as a single standing 
wave of given wavelength (or wave vector of magnitude k — 1/A), and therefore its momentum is exactly 
defined by the de-Broglie equation, whereas the position is totally undefined. In the second case, the 
position is exactly defined, and correspondingly the wavelength or momentum of the particle is totally 
undefined. In most cases, both momentum and position are defined within a finite uncertainty. In this 
case, as mentioned above, a wave package with wavelengths within a certain interval is necessary to 
describe our moving particle. 
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—► 

k = p/h 



Figure 5.9: The Heisenberg uncertainty principle - the two extreme cases 
(either momentum (left) or position (right) are exactly defined) 
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Example 5.1: 

Calculate the wave lengths for the following cases: (i) neutrons thermally excited at T = 500 K, 
(ii) electrons accelerated by a voltage of 50 kV 

Solution: 

We use the de-Broglie relation, = ~ , and determine the momentum via the respective kinetic 

1 9 p 2 

energies, E kin = -mv = —, which have to be identical to the thermal energy (case (i)) or the 
electrostatic energy (case (ii)), respectively. 

(i) For the thermal excitation of the neutrons, we find 



p= /- kT ■ 2m = I- ■ 1.38 ■ l(r 23 ■ 500 ■ 2 ■ = 5.874 ■ 10“ 24 — 

r 2 6 - 10 23 m m 

=> A -- - m - 1.128 ■ 10 -10 m - 0.113 nm 

p 5.874-10 -24 

(ii) For the electrons, we find 


p = y/V^lT^2rn = Vl.6 ■ lO" 19 ■ 50000 ■ 2 ■ 9.109 ■ KT 31 — = 1.207 ■ 10“ 22 — 

m m 

=> A = - = 6 ' 626 ' 10 ~^ m = 5.490 ■ 10~ 12 m = 0.005 nm 

p 1.207-10 -22 

The very small wave length of the electron corresponds to the atomic resolution of the electron 
microscope. On the other hand, neutrons with wave length comparable to the length of chemical 
bonds are used for example to explore the molecular structure of nanoparticles and polymer 
chains by so-called neutron scattering experiments. 


We continue with the mathematical solutions of some fundamental problems of quantum mechanics 
in the next section. 
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5.3 Mathematical solutions of some simple problems in quantum mechanics: 
particle in a box, harmonic oscillator, rotator and the hydrogen atom 

The basis to mathematically solve quantum mechanical problems is the Schrodinger-equation 

fiv (x) =E- X ¥ (x) (Eq.5.12) 

with the Hamilton operator, ^(x) the wave function describing the probability (given as Iff^x)! 2 see 
Eq.5.12) to find a particle or system at position dx, and E the discrete energy values the system may 
assume. Here, it should be noted that at very high energies, the regime of classical physics is reached 
where the spectrum of energies a system may assume is a continuum. The operator H is derived from 
classical physical expressions for kinetic and potential energy following a comparatively simple procedure, 
which will not be discussed in detail in this book. Note that in mathematics any “operator” operates on 
the function written behind it, e.g., the differential operator d/dx operates on a function f(x) (d/dx f(x)). 
We will limit our treatment of quantum chemistry instead to the presentation and brief discussion of the 
quantum chemical results for ^(x) and E for the most fundamental simple models. A very important 
result of the Schrodinger equation is the quantization of energy, which already can be understood 
qualitative from the concept of standing waves (see section 5.2). Note that the wave functions ^(x) 
can assume positive, negative or even complex numbers. Therefore, the probability to find a system at 
position x, which has to be a real number, is generally defined as 

V(x)-V(x)* = |V(x)| 2 (Eq.5.13) 

(i) Our first example is the particle in a box. This model describes any moving particle of mass m limited 
to a very small defined volume by infinitely high potential energy barriers. Within these barriers, the 
particle has only kinetic energy. The potential energy therefore has only the effect to restrict the particle 
position. This fundamental simple model can be used, for example, to explain the change in light 
absorption wavelengths of aromatic molecules from benzene to anthracene (see fig. 5.10): in this case, 
the moving particles are the tc- electrons, and the box dimensions are defined by the size of the respective 
molecule. The larger the molecule or the larger the box, the longer the wavelengths of the standing waves 
to describe the position of the electrons within the box. Consequently, to lift an electron from its ground 
state to an excited state by the absorption of light, a photon of lower energy or larger wavelength is 
needed the larger the molecule. Therefore, from benzene to anthracene the absorption maximum shifts 
to larger wavelengths (or, considering colors within the spectral regime, from UV to blue). 


Figure 5.10: Particle in a box model and the n-electron system in condensed aromatic rings. The "box dimension" 
increases from benzene (left) to anthracene (right), and accordingly the absorption spectrum is red-shifted. 
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In case the potential barrier is infinitely high, the probability to find the particle at each wall is zero, 
and V(x) is given by a perfect standing wave (or sine-, cosine-function). With decreasing wavelength, 
the energy increases, since, according to classical physics and the de-Broglie equation (Eq.5.10.), the 
kinetic energy is given as: 

E = -mv 2 = — = —— (Eq.5.14) 

2 2m 2m 

Since for our particle in a box the wavelengths from the ground state to the excited states decrease as 1, 
1 /2,1/3,..., we expect the energy levels to scale as 1,4,9,..., even without exactly solving the Schrodinger 
equation. The wave functions V(x) and the corresponding E- eigenvalues are given in figure 5.11: 
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Figure 5.11 : wave functions and energy Eigenvalues of 1 -dimensional particle in a box-problem 


Here, + and - define the symmetry of the wave function in respect to the center of the box: + is 
symmetric, - antimetric. 

According to the respective changes in wavelength of the standing waves within the box with increasing 
quantum number (see above), the energies depend on an integrate quantum number E n ~n 2 , n = 1,2 ,.... 
Note that the lowest energy level is not zero, which is obvious from fig. 5.11 since the corresponding 
wave function has a finite wavelength. This finite energy of the ground state is also in agreement with 
the Heisenberg relation, since zero energy would correspond to an exactly defined momentum. In 
combination with a particle position limited within the box, this would, in violation of the Heisenberg 
uncertainty relation, lead to A p x ■ Ax = 0. 


For a 3-dimensional box, the energy levels depend on a set of three independent quantum numbers, the 
particle mass, and the box size, as: 


^n x ,ny,n z g 


~~I ' (% 2 + n y 2 + n z 2 ) 


(Eq.5.15) 


Consequently, there exist different sets of quantum numbers with identical energy levels, for example 
El,1,2 — £*1,2,1 — £*2,1,1 • This phenomenon, called degenerate energy levels, is, for example, also found 
for the p-orbitals of an electron orbiting the hydrogen nucleus, as we will see later. 
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Finally, if the potential walls of our box are not of infinite energetic height, there exists a finite probability 
for the particle to exist outside the box or pass through the wall even if its kinetic energy is still much 
lower than the potential barrier. Note also that in this case the wave function no longer is a simple 
standing wave (simple sine- or cosine-function, see fig. 5.11), but approaches asymptotically 0 beyond 
the walls of the box (see fig. 5.12). This phenomenon is called tunneling, another effect besides energy 
quantization and Heisenbergs uncertainty principle not encountered in classical physics. 



Figure 5.12: particle in a box with finite potential barrier. Note that approaching the barrier the energy levels 
are conning closer. Above the barrier, quantization of energy is no longer found but any energy value is 
physically possible (energetic continuum) 

Formally, the simple model of a particle in a box already leads to some features also found in the more 
complex hydrogen atom. In both cases, the spatial location of a moving particle is limited by a potential 
energy barrier, leading to such general features as energy quantization, or energetically degenerated 
states if the box is 3-dimensional. In addition, the shape of the wave functions already reminds one of 
electron orbitals: the totally symmetric wave function flfl (x, y, z), for example, represents an ls-orbital 
in shape, whereas the degenerate axially antimetric wave functions V 2 t i t i(x, y, z), 2j i(x,y, z) or 
X V 1 12 (x,y,z) correspond in shape to 2p x -, 2p y - or 2p z -orbitals, respectively. The more complex wave 
function ^ 2 , 2,1 ( x > V > z ) finally looks similar to a 3d xy -orbital. 
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(ii) Our second example is the so-called harmonic oscillator, important to determine the energy levels 
of the oscillations of chemical bonds within molecules. Like in a simple spring model, where the force 
pulling back the spring is described by Hookes law (F = —kx, with k the spring constant), the potential 
energy of this system is given as: 

V = - / Fdx = \ kx 2 (Eq.5.16) 


The Schrodinger equation yields the following solutions for the discrete energy levels: 


E v = h ■ v 0 ■ (v + i), v = 0,1,2,... 


(Eq.5.17) 


with basic oscillation frequency v 0 depending on the spring constant and the moving mass as 


v 0 = 


(Eq.5.18) 


For a simple molecule consisting of two atoms only, k corresponds to the stiffness (or force constant) of 
the chemical bond, whereas the mass is the so-called reduced mass defined as (see fig. 5.14): 


m = 


m i 'EM 2 

7711+771 2 


(Eq.5.19) 


/ 

/ 


/ 

/ 

/ 



4 ► 


/ 

/ 

/ 

Figure 5.13: Harmonic oscillator (left) and relation to simple molecular vibrations of a diatomic molecule (right) 

For HC1, for example, you get a reduced mass of 0.97 g/mol. This means that the heavier Cl-atom is nearly 
not moving at all, formally representing the solid wall in figure 5.13, against which the much lighter 
H-atom is vibrating. Energy levels and wave functions of the harmonic oscillator are shown in figure 5.14. 
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Figure 5.14: energy Eigenvalues (a)), wave functions (b)), and squared amplitudes (c)) of the harmonic oscillator 
(c) dotted lines = classical probabilities) (from: Gerd Wedler und Hans-Joachim Freund, Lehrbuch der Physikalische 
Chemie, p. 548, 6.Auflage, Weinheim 2012. Copyright Wiley-VCH Verlag GmbH & Co. KGaA. Reproduced with permission.) 
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In analogy to the particle in a box, the ground state has a non-zero energy Note that the energy levels 
here are equi-distant, which is plausible since the potential barrier is getting wider with increasing 
energy To apply this concept to real molecular vibrations, we have to take into account that the harmonic 
potential at high energies formally allows for negative bond lengths (interatomic distance x < 0), and also 
ignores dissociation of the molecule or bond breaking at very high energies. A more realistic potential 
is therefore the Morse-potential considering these effects (see fig. 5.15): 

V = V 0 ■ (1 - exp(-a ■ x)) 2 (Eq.5.20) 



In Eq. 5.20, Vq provides a measure for the bond strength (or depth of the potential well = dissociation 
energy), and the parameter determines the reciprocal bond length. In contrast to the harmonic oscillator 
potential, the Morse potential is getting wider with increasing energy. As a consequence, the energy 
levels of the inharmonic oscillator are not any longer equidistant but the energy spacing gets smaller 
the higher the quantum number. Finally, if the bond is broken at E > V 0 , the energy spacing becomes 
zero, or there exists no longer any energy quantization, but we have freely moving atoms. 

(iii) Our next example is the so-called stiff rotator, which corresponds to a particle moving on a spherical 
orbit. In this respect, the model is similar to the particle in a box, i.e. the moving particle has no potential 
energy but is limited to angular positions between 0° and 360°. Consequently, the energy levels again 
scale with the quantum numbers squared, E m ~m 2 . The main difference to the particle in a box is that 
m = 0, 1, 2, i.e. here the ground state is Eq = 0! The corresponding wave function at all rotational 
angles has a constant amplitude. Therefore, whereas energy and momentum are well-defined (=0!), the 
positional probability is, in agreement with the Heisenberg-relation, totally undefined. 
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Figure 5.16: Sketch of the wave function of the ground state for the particle on a ring 

Alternatively, one may consider the wavelength of the ground state wave function to be infinitely 
large, leading according to the de-Broglie equation to a momentum of zero, and therefore also to zero 
kinetic energy. 

(iv) We conclude this section about simple quantum-mechanical models by briefly discussing the solutions 
of the Schrodinger-equation for the hydrogen atom. The potential energy within the Hamiltonian here 
is given by the Coulomb-attraction between positively charged nucleus and negatively charged moving 
electron. Further, due to its very high mass and corresponding momentum of inertia we may consider 
the nucleus, in respect to the fast moving electron, as stationary. Note that, in contrast to our previous 
examples, here we have a negative potential whose variation with electron-nucleus-distance is given as 

(Eq.5.21) 

Consequently, in contrast to the particle in a box we also expect negative energy levels, asymptotically 
increasing towards zero with increasing quantum number. In the Bohr model presented at the beginning 
of this chapter we already have seen that 

E n ~-^ 2 , n= 1,2,3,... (Eq.5.22) 

From what we have learned so far about quantum mechanics and standing waves, it is very plausible that 
the energy spacing for the hydrogen atom strongly decreases with increasing quantum number, since 
the potential barrier in this case is widening much stronger with increasing energy than the parabolic 
potential barrier of the harmonic oscillator (where we find a constant energy spacing or equidistant 
energy levels, see fig. 5.14). 
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E~ —rr(see Eq.5.3) 

n z 


Figure 5.17: Energy levels of the hydrogen atom depending on quantum numbers, red = hyperbolic Coulomb potential keeping 
the electron close to the positive core, s, p and d correspond to quantum state of the angular momentum of the electron. Note 
that for hydrogen, the energy eigenvalues only depend on the main quantum number n! 
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Going into more detail, there is not only the main quantum number n defining the energy of the 
electronic state of the hydrogen atom (corresponds to K-, L- or M-electronic shell, respectively), but 
also the quantum number of rotational momentum defining the shape of the orbital (corresponds to s-, 
p- or d-orbital), as well as the corresponding magnetic quantum number describing the orientation of 
non-isotropic electron orbitals in space (p x , p y or p z ). For the hydrogen atom, these energetic states are 
degenerate in the sense that the energy level only depends on the main quantum number w, i.e. 3s, 3p 
and 3d-orbitals all have identical energy levels. For more complex multi-electron atoms, this energetic 
degeneration is not found any longer. 


5.4 A brief introduction to optical spectroscopy 


e 2 


AAAAa 


E, 


e 2 





Figure 5.18: Schematics of absorption (left) and emission (right) spectroscopy as energetic transitions between energy 
eigenvalues E 1 and E 2 , involving either absorption or emission of light of energy A E = E 2 -E 1 =hv = h c ~. 


Optical spectroscopy deals with the absorption or emission of electromagnetic radiation in connection 
with a transition between energetic states of molecules. Depending on wavelength and therefore energy 
of the radiation, one may address transitions between electronic states (UV/Vis), molecular vibrations 
and rotations (IR and microwaves), or magnetic states of the nucleus (radio waves), which will be 
explained in more detail below. The spectrum in terms of probability of light absorption in dependence 
of the wavelength of the incident light then directly corresponds to the energy eigenvalues and therefore 
is characteristic for a given molecule. Consequently, the spectrum may be used as a finger print in 
qualitative chemical analytics. For quantitative analytics, one can determine the absorption of light 
at given wavelength, and the concentration of a substance in solution is then given by the law of 
Lambert-Beer: 

log = —e ■ c ■ d (Eq.5.23) 

with I 0 the intensity of the incident light, the intensity of the light passing an optical cuvette of thickness 
d, e, the molar decadic extinction coefficient specific for a given molecule and wavelength, and the 
molar concentration of this molecule in solution. 
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As shown in fig. 5.18, the energy of the electromagnetic radiation absorbed to trigger the transition 
between different energetic states of the molecule is given as: 

A E = E 2 — E 1 = h'V = h m j = h'C'V (Eq.5.24) 

with v the wave number of the incident radiation, often used as the important parameter in optical 
spectra. In contrast to the wavelength A, v is directly proportional to the energy of the transition. Wave 
lengths, energies and frequencies of electromagnetic waves are summarized in figure 5.19. 

wavelength /nm 10 200 400 750 


Xray 

far UV 

near UV 


IR 







wave number/cm' 1 10 6 5-10 4 2,5-10 4 1,3-10 4 

frequency/s' 1 3-10 16 15-10 14 7,5-10 14 4-10 14 

energy/kJ/1 Mol photons 11970 598,7 299,4 159,1 

Figure 5.19: wave lengths and energies of the electromagnetic spectrum 


In table 5.1, some important spectroscopic methods including energy regime and detected molecular 
properties are summarized: 


X 

V 

v/crrr 1 

AE/ kJ/mol 

molecular transition 

method 

100 m-1 m 

3 MHz-300 MHz 

io- 4 -o.oi 

10-6-1 o- 4 

nuclear spin 

NMR 

1 cm-100 pm 

30 GHz-3-10' 2 Hz 

1-100 

0,01-1 

molecule rotation 

microwave 

100 pm-1 pm 

3-10 12 Hz-3-10 14 Hz 

100-10 4 

1-100 

molecule vibration 

IR, Raman 

1 pm-10 nm 

3-10' 4 Hz-3-1 O ' 6 Hz 

10 4 -10 6 

100-10 4 

outer (binding) electrons 

UV/Vis, fluorescence 


Table 5.1: some important spectroscopic methods 


Next, let us consider some important principles of optical absorption spectroscopy: 


(i) For molecules, each electronic state corresponds to a given strength and length of the chemical bond, 
and therefore each electronic energy level is further divided energetically into vibrational states (= change 
of bond lengths and/or bond angles due to molecular vibrations) which themselves are divided into 
rotational states (= rotation of the whole molecule at given (average) bond lengths and bond angles): 
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Figure 5.20: Energy levels of molecules: each electronic state (black) contains a set of 
vibronic states (red), which itself contains a set of rotational states (blue). 
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To excite only the molecular vibrations without simultaneous electronic excitation, incident photons 
of lower energy or longer wave length are needed (infra red (IR)). Further, to excite only molecular 
rotations even longer wave lengths are necessary (micro waves). On the other hand, a molecular electronic 
transition excited by UV/Vis-absorption usually consist of simultaneous vibrational and rotational 
excitations, wherefore molecular UV/Vis-spectra, in contrast to the spectra of atoms, do not show sharp 
lines but characteristic broad absorption bands, if the absorption is plotted versus the wavelength of the 
incident light. 

(ii) The Franck-Condon-principle: during the excitation of electronic states of a molecule, typically caused 
by absorption of UV/Vis-light, the positions of the atomic cores remain unchanged because of the short 
time scale or high frequency of the incident light (~ 10 -15 s). For most molecules, the excited electronic 
state in respect to that of the ground state is not only corresponding to higher energy but also to larger 
interatomic distances or bond lengths (see fig. 5.21), since electronic excitation typically causes a decrease 
in strength and stiffness of the chemical bond. The probability of the electronic transition triggered 
by light absorption then is proportional to the overlap of the ground vibrational wave function of the 
electronic ground state, and the vibrational wave functions of the electronically excited state. Note again 
that, as described in the preceding paragraph (i), not only the electronic state but also the vibrational 
state of the molecule changes upon absorption of UV/Vis-light! Due to the rapid radiation-less decay of 
excited molecular vibrations, the emitted fluorescent light is usually red-shifted (= longer wave length 
or lower transition energy) in respect to the wavelength of absorption (see fig. 5.21). 



bond length): the transition has the largest probability at maximum overlap of the vibrational wave 
functions of the electronic ground and the excited electronic state. This, in combination with a rapid 
radiation-less decay of the molecules vibrations (dotted black arrows), explains why the emitted light is 
shifted towards longer wave length (in respect to the absorbed light). 
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(iii) Molecular energetic transitions can only be triggered by the absorption of light if not only the 
wavelength is fitting to the energy difference between ground and excited state, but also if the electric 
dipole momentum of the molecule is changing with the transition. This can be understood on a classical 
physical basis if we consider light as an electromagnetic wave, and light absorption as an energetic 
interaction of the electric field vector of the incident light with the electric dipole momentum of the 
molecular transition. For example, the symmetric oscillation of a C0 2 -molecule cannot be detected by 
absorption of IR-light, since the dipole momentum does not change for this mode of molecular vibration. 
On the other hand, the asymmetric vibration of the C0 2 -molecule can be detected by IR-absorption 
since here the vibrational motion leads to a change of the dipole momentum. 


As an illustrative example how to extract molecular characteristics from an optical spectrum, let us 
consider the rotation-vibration spectrum of HC1 in the gaseous state at low pressure. Only in this case 
the energetic levels are well defined and allow the separate detection of individual rotational levels 
within the vibrational excitation, whereas at higher pressure or even in the condensed phase peak 
broadening, mostly due to intermolecular interactions and collisions, leads to a comparatively blurred 
broad absorption band with no spectral fine structure. The allowed transitions and the resulting spectrum 
are shown in figure 5.22. 



A J = -l AJ = +1 



Figure 5.22: Energy transitions (right) and corresponding rotational-vibrational spectrum of HCI (left). Dotted line corresponds to the 
hypothetical transition aj = o not found in the HCI-spectrum. 
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The selection rules of allowed energetic transitions for the spectrum shown in fig. 5.23 are Av = +1, 
i.e. vibrational excitation by one energy level, and AJ = ±1, i.e. rotational excitation with a change 
in rotational quantum number by ±1. These lead to a spectrum containing two series of peaks, one 
corresponding to at AJ = +1 higher wave number, and the other corresponding to AJ = —1 at lower 
wave number. The gap in between, formally corresponding to AJ = 0, reflects a purely vibrational 
transition, and the corresponding wave number allows to directly extract the bond stiffness (see Eq.5.17). 

On the other hand, the energy spacing between the two peaks closest to this gap allows us to extract the 
bond length, since it is directly related to the rotational energy or momentum of inertia. Note that the 
selection rule is AJ = +1 based on the conservation of the rotational momentum, a principle also known 
in classical physics: the absorbed photon has a rotational momentum of 1, and therefore the loss of this 
momentum upon light absorption has to be compensated by a corresponding change in the rotational 
state of the molecule. (Note that the rotational momentum of the photon is related to a so-called spin, 
and also electrons and some atomic nuclei have a spin (the single electron, for example, has spin Vi). 
Some very important spectroscopic methods, for example nuclear magnetic resonance (NMR, see also 
tab. 5.1) which is also used in medicine for imaging, are related to pure spin transitions. The interested 
reader is referred here to Physical Chemistry textbooks for details). 
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Finally, one should note the interesting intensity distribution within a peak series: at room temperature, 
typically the 3 rd peak shows the highest absorption. This is caused by the fact that, according to the 
Boltzmann distribution, the 3 rd rotational state has the highest probability Finally, one should note that 
the absorption peaks in the spectrum actually are double-peaks, which is due to the fact that Cl-atoms 
are found in form of two isotopes of slightly different atomic masses. 

The mass of vibrating atoms at given bond strength may have a strong influence on the energy of the 
vibrational states (see Eq.5.17), or the peak position in an IR-absorption spectrum or, alternatively, a 
Raman scattering spectrum. For the later, the vibrational energy transition of the molecule corresponds 
to an energetic shift of the scattered light in respect to the incident wavelength. For a more detailed 
discussion of Raman spectroscopy, the interested reader should consult one of the textbooks given at 
the end of this compendium. 

The effect of atomic mass on the vibrational energy is best illustrated by the IR- or Raman-spectra of 
organic molecules containing either hydrogen or deuterium. The doubling in mass has a correspondingly 
strong isotope effect in the spectrum, i.e. a shift in transition energy by a factor of V2 (see Eq.5.17, and 
figure 5.23). 
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Figure 5.23: Raman-spectra of CHCI 3 (top) and CDCI 3 (bottom) showing the isotope effect (see peaks at highest wave 
number as indicated by arrows). Both spectra were measured with an iRaman-spectrometer by W. Schartl at university 
Mainz, Germany, August 2013. 
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Example 5.2: 

In the high-resolution IR spectrum of HC1 taken for a gaseous sample at low pressure, we find 
a band gap at frequency 8.66 ■ 10 -13 Hz, and an inner spacing between the two series of 
rotational transitions of 0.12 ■ 10 _13 Hz (see also fig. 5.22). Calculate the bond strength and 
bond length of the molecule. 


Solution: 

(i) bond strength 

For the HC1 molecule, the needed reduced mass is 


Mh ' Mg 
M h + M cl 


1-35 
1 + 35 


g/Mol = 0.98 g/Mol 


Note that this is almost just the mass of the hydrogen. Substituting the midpoint frequency into 
the expression containing the bond force constant gives: 


2n -v = 6.28 ■ 8.66 ■ 10 13 s -1 = 

0.98-HT 3 , 

k = (6.28 ■ 8.66 ■ 10 13 ) 2 ■ „ N/m = 481 N/m 

6.02 ■ 10 23 


k 

ft 


(ii) bond length 

The separation between the two series of spectral peaks is twice the rotational energy change 
from quantum number J = 0 to J = 1. This energy is related to the momentum of inertia, i.e. 



h- 0.12 ■ 10 13 s -1 


with I = g - r 2 


Therefore, the bond length is given as: 


r 2 = 


<£) 


f 6.626-10 -34 \ 

V 6.28 J 


ju-h-0.12-10 13 s~ 1 


!' I !! n 2 r 6 - 626 ' 10 ~ 34,0 - 12 ' 1013 

6 . 02 - 10 2! - :5 


m z = 1.72 ■ 10 20 77i 2 , or r — 0.13 nm 


We conclude this chapter (and our booklet) by reviewing the spectroscopic transitions between electronic 
states of molecules in more detail. The following scheme called Jablonski-diagram summarizes the 
possible transitions (see fig. 5.24): 
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Figure 5.24: Jablonski diagram showing electronic excitation and different pathways of energy relaxation 
(with and without (dotted) radiation absorption/emission). Circles indicate the electronic configuration of the 
respective electronic states, the arrows indicate the electron spin direction (+V 2 (up) or - V 2 (down)). 

See text for details. 

Upon UV/Vis-irradiation and absorption of light (A), the molecule shows an electronic transition from 
the ground state S Q to an excited electronic state S x (or S 2 ) by lifting one electron (or two electrons) to a 
higher energetic state while conserving the spin orientation. Note that simultaneously the excitation of 
molecular vibrations is usually found (see also fig. 5.21). Next, the excited energy is partially converted 
via a rapid radiation-less process called internal conversion (time-scale: l(f 12 seconds) (IC) into heat, 
before radiation is emitted from the vibrational ground state of S x to an excited vibrational state of S Q (F). 
This emission called fluorescence, on a time-scale of 10 9 - 10' 6 seconds, is therefore shifted, in respect 
to the wavelength of the incident light, towards longer wavelength (red-shift). 

Some molecules, especially those containing heavier atoms like S or P, show a competing process, a 
radiation-less transition from S x to an excited triplet state T x , where the orientation of the spin of the 
excited electron is reversed. This process called intersystem crossing (ISC) is usually forbidden and 
therefore very slow, since it violates the physical principle of the conservation of the rotational momentum. 
However, in case of heavier atoms the change in the electron spin orientation may be compensated by a 
corresponding change in the electrons angular momentum, a phenomenon called “spin-orbit-coupling” 
strongly speeding up the ISC to a time scale so it can compete with fluorescence. From the triplet state, 
the system decays via emission of radiation on a very slow time scale (1-100 s), and this emission is 
called phosphorescence (Ph). 
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